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FOREWORD 


The  research  described  in  the  report  involves  presentation  and  evaluation 
of  results  from  some  150  triaxial  compression  tests  run  on  specimens  of 
unfrozen  and  frozen  sand  at  varying  densities,  levels  of  confinement,  rates  of 
axial  strain  and  temperatures.  The  authors  had  plarmed  to  prepare  a  very  brief 
report  that  summarized  the  principal  findings  and  conclusions  from  three 
papers:  one  on  the  behavior  of  frozen  sand  at  small  strains  (from  0.001%  to 
about  1%  axial  strain);  a  second  on  the  behavior  of  frozen  sand  at  large  strains 
(where  fracturing  of  the  ice  matrix  and  the  frictional  resistance  of  the  sand 
skeleton  play  major  roles);  and  a  third  on  quantitative  analysis  of  the  peak 
strength  of  frozen  sand  using  a  dilatancy-hardening  model.  Whereas  the  first 
paper  was  written  last  spring  and  has  been  accepted  for  publication,  it  became 
obvious  that  considerably  more  analysis  of  the  data  was  required,  in  order  to 
submit  the  other  two  papers.  Consequently,  the  first  two  authors  devoted 
several  himdred  hours  since  last  July  to  further  evaluation  of  the  data  and  on 
preparation  of  this  report.  But  in  spite  of  their  effort,  the  report  is  less  polished 
than  desired,  which  explains  the  many  hand-drafted  figures  in  the  appendices, 

FoT 

and  why  the  text  on  large  strain  behavior  sometimes  reads  like  a  mystery  novel. 
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However,  we  hope  the  final  section  provides  a  coherent  summary  and  set  of 
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conclusions. 
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1.  INTRODUCTION 


1.1  BACKGROUND 

Frozen  soils  exist  in  Arctic  regions  as  perennially  frozen  soil  (permafrost)  that 
underlies  nearly  20%  of  the  land  surface  of  the  earth.  Of  all  geomaterials,  the 
mechanical  properties  of  frozen  soil  are  perhaps  the  most  variable  and  difficult  to 
understand  and  model.  Consequently,  the  placement  of  structures  on  frozen  soil 
poses  unique  design  problems  because  of  uncertainties  regarding  the  time  and 
temperature  dependence  of  the  frozen  soil  response.  Similar  problems  also  arise 
when  artificial  soil  freezing  is  used  as  a  method  of  ground  improvement  for 
underpinning  foundations  and  for  strengthening  soil  around  excavations  (tunnels, 
shafts,  etc.). 

Even  the  simplest  form  of  frozen  soil  (sand  with  ice-filled  pores)  entails  a 
highly  complex  interaction  between  the  skeleton  of  solid  particles  and  the  pore 
matrix  (composed  of  ice  and  unfrozen  water)  that  changes  continuously  with  time 
as  a  function  of  temperature  and  stress/ strain  level.  Consequently,  the 
development  of  reliable  constitutive  relationships  to  model  the  stress-strain-time- 
temperature  behavior  of  frozen  soils  lags  far  behind  modeling  capabilities  for  other 
geomaterials  such  as  soil  and  ice.  Long-term  progress  in  modeling  requires  research 
in  at  least  three  principal  areas: 

1)  experimental  studies  to  completely  define  the  strength-deformation 
properties  of  representative  frozen  natural  soils  as  a  function  of  the 
major  variables  thought  to  affect  their  behavior  at  the  macrolevel; 
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2)  coupled  experimental  and  theoretical  studies  to  test  hypotheses 
regarding  the  physical  mechanisms  controlling  the  behavior  of 
frozen  soils  at  various  stages  of  deformation;  and 

3)  theoretical  studies  which  absorb  the  insight  from  research  areas  1) 
and  2)  and  lead  to  the  development  of  physically  based  constitutive 
relationships  for  frozen  soils. 

Regarding  the  first  area,  the  major  variables  thought  to  affect  the  behavior  of 
frozen  soils  can  be  divided  into  three  categories.  Their  description  and  the  scope  of 
experimental  data  that  existed  prior  to  this  research  are  summarized  below. 

•  The  soil  type  that  will  be  frozen.  In  general,  this  means  the  particle 
size,  shape  and  gradation  of  granular  soils  and  the  plasticity 
characteristics  of  cohesive  soils.  Most  experimental  data  are 
restricted  to  uniform  sands,  such  as  commercially  available  Ottawa 
sand. 

•  The  initial  state  during  freezing.  This  refers  to  the  initial  density 
(granular)  or  stress  history  (cohesive)  of  the  specimen,  the  nature  of 
the  pore  fluid,  the  magnitude  of  the  applied  effective  confining 
stress,  and  the  freezing  temperature.  Most  existing  data  are  for  sand 
specimens  having  one  density,  and  all  specimens  were  frozen 
under  zero- confining  stress. 

•  The  conditions  during  shear.  The  major  variables  include  the 
confining  pressure,  the  rate  of  straining  or  the  creep  stress,  the 
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temperature,  and  the  mode  of  shearing.  Essentially  all  existing  tests 
have  been  restricted  to  monotonic  loading  in  triaxial  compression 
with  external  measurement  of  axial  strain  (which  precludes  reliable 
data  at  small  strains)  and  only  one  program  varied  confining 
pressure,  strain  rate  and  temperature. 

This  review  of  prior  research  indicates  that  no  test  program  has  produced 
stress-strain  data  from  very  small  to  very  large  strains  for  any  frozen  soil  as  a 
function  of  the  four  principal  variables  that  affect  its  response,  namely  specimen 
density,  confining  pressure,  strain  rate  and  temperature. 

1.2  RESEARCH  OBJECTIVES  AND  SCOPE 

The  long-term  objective  of  MIT's  research  is  to  develop  physically  based 
constitutive  relations  to  describe  the  stress-strain-time-temperature  behavior  of 
frozen  soils,  which  in  turn  first  requires  a  better  understanding  of  their 
fundamental  behavior.  Consequently,  the  near-term  objective  is  the  identification 
of  the  principal  physical  mechanisms  that  control  the  stress-strain  behavior  of 
frozen  soils  and  a  quantitative  (or  at  least  qualitative)  assessment  of  their  relative 
importance.  The  approach  used  by  MIT  to  achieve  this  near-term  objective  has 
three  components: 

1)  conduct  a  comprehensive  experimental  program  to  precisely 

measure  the  behavior  of  a  frozen  natural  sand  over  a  wide  range  of 
testing  conditions; 


2)  from  a  knowledge  of  the  general  behavior  of  polycrystalline  ice  and 
of  the  measured  stress-strain  behavior  of  the  same  sand  in  an 
unfrozen  state,  attempt  to  deduce  the  relative  importance  of 

•  the  ice  matrix  per  se 

•  the  frictional  resistance  of  the  sand  skeleton  per  se,  and 

•  the  interaction  between  the  ice  matrix  and  the  sand  skeleton; 

3)  where  possible,  employ  models  for  composite  materials  to  quantify 
the  physical  mechanisms  responsible  for  the  macro-behavior  of 
frozen  sands. 

Detailed  results  from  the  preceding  ARO  contract  (DAAL03-89-K0023)  for  the 
period  January,  1989,  through  December,  1991,  are  contained  in  Andersen  (1991)  and 
Andersen  et  al.  (1992).  The  principal  experimental  component  of  that  research  used 
a  high-pressure,  low-temperature  triaxial  compression  testing  system  with  on- 
specimen  axial  strain  measurements  and  lubricated  end  platens  to  characterize  the 
stress-strain-volume  change  characteristics  of  frozen  Manchester  Fine  Sand  (MFS)  at 
a  temperature  of  T  =  9.5°C  from  very  small  (0.001%)  to  very  large  (25%)  axial  strains. 

The  program  varied  the  relative  density  of  the  sand  from  D^  =  20  to  100%,  the 
confining  pressure  from  <s^  =  0.1  to  10  MPa,  and  the  strain  rate  from  e  =  3x10^/ s 
(slow)  to  SxlO’^^/s  (fast).  A  similar  triaxial  system  was  used  to  perform  an  initial 
series  of  undrained  and  drained  shear  tests  on  unfrozen  MFS.  The  research  also 
used  a  model  developed  for  composite  materials  (concrete)  to  evaluate  the  Young's 
modulus  of  frozen  MFS. 
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The  experimental  research  conducted  under  this  ARO  contract  extended  the 
prior  data  on  frozen  MFS  to  lower  temperatures  (nominal  values  of  T  =  -15°C,  -20°C 
and  -25°C),  completed  the  shear  testing  program  on  unfrozen  MFS,  and  conducted 
preliminary  tests  on  MFS  wherein  specimens  of  dense  sand  were  frozen  at  T  =  -9.5°C 

after  being  consolidated  to  stresses  ranging  from  o'c  =  0.1  to  10  MPa.  This  research 

also  evaluated  a  dilatancy-hardening  model  developed  by  Ladanyi  (1985)  for 
predicting  the  strength  of  frozen  sand. 

Unique  aspects  of  the  triaxial  compression  shear  testing  program  on  MFS 
include  the  following  features: 

•  the  first  set  of  data  on  any  frozen  soil  to  fully  characterize  behavior 
as  a  function  of  sand  density,  confining  stress,  strain  rate  and 
temperature; 

•  the  first  set  of  data  to  reliably  measure  axial  stress  vs.  axial  strain 
behavior  from  very  small  (0.001%)  to  very  large  (25%)  axial  strains, 
along  with  reasonably  accurate  measurements  of  volumetric 
strains; 

•  the  first  set  of  data  wherein  specimens  were  consolidated  to  varying 
pressures  prior  to  freezing,  although  the  results  are  rather 
surprising  and  should  be  checked; 

•  the  first  program  to  measure  the  shear  behavior  of  a  soil  in  both  its 
frozen  and  its  unfrozen  states. 
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Consequently,  even  if  some  of  the  conclusions  drawn  by  the  authors  from 
their  analyses  of  the  data  prove  to  be  questionable,  other  researchers  can  still  benefit 
from  their  own  interpretation  of  a  unique  set  of  experimental  data. 

1.3  ORGANIZATION  OF  THIS  REPORT 

The  objective  of  this  report  is  to  provide  a  detailed  summary  of  the  most 
important  experimental  results  and  of  the  authors'  conclusions  regarding  the 
relative  importance  of  the  physical  mechanisms  controlling  the  strength- 
deformation  behavior  of  frozen  sands.  The  technical  summary  includes  research 
conducted  under  both  ARO  contracts,  i.e.,  support  during  calendar  years  1989  to  1991 
and  from  June,  1992,  through  May,  1994. 

Section  2  presents  an  overview  of  the  scope  of  the  experimental  program  on 
Manchester  Fine  Sand  (MFS)  and  of  the  experimental  procedures  since  the  details 
are  contained  in  two  doctoral  theses  (Andersen,  1991  and  Swan,  1994)  and  the  Final 
Report  for  the  preceding  contract  (Andersen  et  ah,  1992).  Section  3  describes  the 

behavior  of  frozen  MFS  at  small  strains  (Eg  <  1%),  where  sand  density  and  confining 
pressure  are  seen  to  play  a  minor  role  in  the  observed  behavior.  Section  4  then 
describes  the  large  strain  behavior,  wherein  the  fracturing  of  the  ice  matrix  (as 
inferred  from  dilatant  volumetric  strains)  and  the  frictional  resistance  of  the  sand 
skeleton  now  play  a  major  role  in  the  observed  behavior.  Section  5  summarizes 
results  from  the  consolidate-freeze  tests  on  MFS,  i.e.,  influence  of  effective  confining 
stress  during  freezing  on  subsequent  behavior.  Section  5  also  applies  Ladanyi's 


dilatancy-hardening  model  to  frozen  MFS,  which  incorporates  results  from  triaxial 
compression  shear  tests  on  unfrozen  sand  specimens. 


1.4  PUBLICATIONS  AND  REPORTS 

Three  publications  and  reports  on  frozen  sand  were  sent  to  ARO  in  1994:  on 
January  31,  a  Technical  Progress  Report  was  submitted;  in  August,  a  copy  of  the 
doctoral  thesis  by  Swan  (1994)  was  submitted;  and  in  May,  we  sent  a  copy  of  the 
paper  titled,  "Small-Strain  Behavior  of  Frozen  Sand  in  Triaxial  Compression,"  by 
Andersen  et  al.  (1994)  that  was  submitted  to  the  Canadian  Geotechnical  Journal. 
Although  accepted  for  publication,  the  Andersen  et  al.  paper  was  expanded  and 
improved  in  response  to  detailed  comments  by  two  reviewers;  ARO  was  sent  a  copy 
of  the  revised  paper.  Two  other  papers  are  planned  for  submittal  to  the  to  the 
Canadian  Geotechnical  Journal  in  1995;  their  status  is  as  follows: 

•  "Large-Strain  Behavior  of  Frozen  Sand  in  Triaxial  Compression"  by 
C.C.  Ladd,  G.R.  Andersen,  C.W.  Swan  and  J.T.  Germaine.  Rough 
drafts  of  text,  tables  and  figures  were  prepared  by  Ladd  that  extend 
the  scope  of  analyses  contained  in  the  theses  by  Andersen  (1991)  and 
Swan  (1994).  Section  4  of  this  report  presents  the  principal  results 
and  conclusions. 

•  "Application  of  a  Dilatancy-Hardening  Model  to  the  Strength  of 
Frozen  Sand  in  Triaxial  Compression"  (tentative  title)  by  C.W. 

Swan,  C.C.  Ladd,  J.T.  Germaine  and  G.R.  Andersen.  The  selection  of 
parameters  and  results  of  preliminary  analyses  are  contained  in 


Swan  (1994);  these  are  summarized  in  Section  5  of  this  report,  along 
with  results  of  additional  analyses. 

1.5  STAFF  AND  ADVANCED  DEGREES 

Dr.  Charles  C.  Ladd,  Edmund  K.  Turner  Professor,  and  Dr.  John  T.  Germaine, 
Principal  Research  Associate,  served  as  co-principal  investigators  for  the  two-year 
project.  Both  are  members  of  the  Henry  L.  Pierce  Laboratory  in  MIT's  Department  of 
Civil  and  Environmental  Engineering.  In  simplified  terms.  Dr.  Germaine  was 
mainly  responsible  for  supervision  of  the  equipment  development  and 
experimental  work,  while  Professor  Ladd  administered  the  contract  and  was  mainly 
responsible  for  supervising  the  analyses  and  reporting  the  results. 

Christopher  W.  Swan  was  a  full-time  graduate  Research  Assistant  on  the 
project;  in  May,  1994,  he  received  the  Doctor  of  Science  Degree  in  Civil  and 
Environmental  Engineering.  He  is  now  an  Assistant  Professor  in  the  Department  of 
Civil  and  Environmental  Engineering  at  Tufts  University  in  Medford, 
Massachusetts. 

The  project  involved  three  undergraduate  students  under  MIT's 
Undergraduate  Research  Opportunities  Program  (UROP).  They  were  Victor 
Salvador  during  the  summer  of  1992,  Munera  Abdur-Rahim  during  the  summer  of 
1993,  and  Adrian  Maung  for  one  year  starting  September,  1992. 
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2.  EXPERIMENTAL  PROGRAM  AND  PROCEDURES 


2.1  DESCRIPTION  OF  MANCHESTER  HNE  SAND 

A  large  supply  of  Manchester  Fine  Sand  (MFS)  was  obtained  from  the  banks 
of  the  Merrimak  River,  Manchester,  New  Hampshire,  and  processed  to  reduce  the 
non-plastic  fines  content  (portion  passing  No.  200  sieve)  from  20%  to  7%.  Pertinent 
characteristics  of  the  MFS  abstracted  from  Andersen  et  al.  (1992,  1994)  and  Swan 
(1994)  follow: 

•  composition — mainly  subangular  quartz  and  feldspar  with  some 
mica  flakes  (Martin  et  al,  1981) 

•  gradation — mean  particle  size,  dso  =  0.18  mm;  99%  passing  No.  45 
sieve  (0.35  mm);  7.2%  passing  No.  200  sieve  (0.074  mm);  uniformity 

coefficient,  C^  =  —  =  2.35 

dio 

•  specific  gravity — Gs  =  2.688  +  0.003 

3 

•  minimum  dry  density — =  1408  kg/m  ; 
maximum  void  ratio,  e^ax  ==  0.909 

3 

•  maximum  dry  density — =  1701  kg/m  ; 

minimum  void  ratio,  e^jn  =  0.580 

•  particle  crushing  after  undrained  shearing  at  10  MPa  corifining 
pressure — negligible  for  frozen  sand;  moderate  for  unfrozen  sand 
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for  intermediate  size  particles,  e.g.,  percent  finer  than  0.125  mm 
increases  from  about  25%  to  35% 

•  unfrozen  water  content — ^less  than  0.14%  at  T  =  9.5°C 

2.2  SCOPE  OF  EXPERIMENTAL  PROGRAM 

The  experimental  program  involved  triaxial  compression  tests  performed  on 
tJ\re0  types  of  specimens:  1)  specimens  that  were  frozen  in  a  mold  before  being 
placed  in  the  triaxial  cell,  i.e.,  the  conventional  method  for  testing  frozen  soil;  2) 
specimens  that  were  frozen  after  being  consolidated  to  varying  effective  stress  levels 
in  the  triaxial  cell,  which  are  called  consolidate-freeze  tests;  and  3)  specimens  that 
were  not  frozen. 

Table  2.1  summarizes  the  testing  matrix  for  the  parametric  study  in  terms  of 
the  number  of  conventional  frozen  sand  tests  at  each  condition.  The  four  variables 
were 

•  the  relative  density  of  the  sand  specimens,  which  varied  from  = 

20  to  100%,  where  D,  =  (e^ax  "  e)/ (Smax  •  emin)/ 

•  the  confining  pressure  during  shearing,  which  varied  from  =  0.1 
to  10  MPa; 

•  the  nominal  testing  temperature,  which  varied  in  5°C  increments 
from  T  =  -10  to  -25°C; 

•  the  nominal  axial  strain  rate,  which  varied  from  e  =  3x10  / s  (slow) 
to  e  =  5x10^/ s  (fast 
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Andersen  (1991)  and  Andersen  et  al  (1992)  presented  results  from  the  tests  at 
T  =  -10°C,  which  included  the  most  detailed  study  of  the  influence  of  sand  density 

.  -5 

and  confining  pressure  at  the  moderate  strain  rate  of  e  =3.5  x  10  /s.  The  current 
contract  performed  the  tests  at  the  lower  temperatures  and  typically  included  loose 
and  dense  specimens  at  low  and  high  confining  pressures  (Oc  =  0.1  and  10  MPa, 
respectively). 

Ten  consolidate-freeze  tests  were  run  on  dense  MFS  consolidated  to  effective 
stresses  of  =  0.1,  2, 5  and  10  MPa,  frozen  at  T  =  -10°C  and  then  sheared  at  the 

moderate  strain  rate.  Section  5.3  describes  the  general  scope  of  the  49  undrained  and 
drained  triaxial  compression  tests  run  on  unfrozen  MFS  for  qualitative  comparison 
with  results  for  frozen  sand  and  to  obtain  parameters  required  for  application  of 
Ladanyi's  dilatancy-hardening  model. 

2.3  COLD  ROOM  FACILITY  AND  TRIAXIAL  TESTING  SYSTEM 
2.3.1  Cold  Room  Facility 

Pertinent  aspects  of  this  facility  include  the  following  features 
(abstracted  from  Andersen,  1991;  Andersen  cf  ah,  1992;  and  Swan,  1994): 

•  a  growth  room  (at  T  =  0°C)  for  specimen  freezing;  a  staging  room  (at 

T  =  -4  to  -5°C)  for  preparing  triaxial  specimens;  and  a  testing  room 
maintained  a  few  degrees  colder  than  the  desired  testing 
temperature;  and 


12 


•  an  environmental  chamber  to  maintain  the  air  temperature 
surrounding  the  triaxial  cell  to  within  ±  0.17°C. 

2.3.2  Triaxial  Testing  System 

Figure  2.1  shows  a  schematic  of  a  Wykeham  Farrance  high-pressure 
triaxial  cell  that  was  modified  for  conventional  frozen  soil  testing  to  accommodate 
lubricated  end  platens,  on-specimen  axial  strain  measurements,  an  internal  load  cell 
to  measure  the  axial  force,  and  thermistors  to  record  temperatures  in  the  silicone  oil 
cell  fluid  near  the  top  and  bottom  of  the  test  specimens  (nominal  dimensions  of  38 
mm  diameter  by  75  mm  long).  Confining  pressures  were  generated  using  a  closed- 
loop  pressure/ volume  control  system  consisting  of  a  DC  servomotor  attached  to  a 
hydraulic  piston  acting  on  a  reservoir  of  fluid.  For  tests  on  unfrozen  sand,  a  similar 
system  provided  back  pressure  saturation  and  volume  changes  during  consolidation 
and  drained  shear.  For  the  consolidate-freeze  tests,  the  top  cap  was  modified  to 
allow  drainage  of  water  as  the  specimen  was  frozen  by  circulating  a  freezing  fluid  at  - 
8°C  through  the  base  pedestal. 

For  conventional  tests  on  frozen  sand,  the  axial  strain  rate  was  applied 
by  a  screw-driven  variable-speed  45kN  Wykeham  Farrance  (T-57)  loading  frame.  As 
a  result  of  the  use  of  lubricated  end-platens  and  compliance  of  the  triaxial  system, 
the  actual  strain  rate  experienced  by  the  specimen  increased  during  shear. 

Compared  to  the  maximum  rate,  which  was  usually  achieved  after  about  one 

percent  axial  strain,  the  initial  secant  value  (up  to  ~  0.03%)  was  about  one  tenth  of 
the  maximum  rate,  and  the  tangent  value  at  the  upper  yield  stress  (defined  in 
Section  2.5  and  occurring  at  Eg  =  0.6  ±  0.3%)  was  typically  10  ±  5%  lower  than  the 
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maximum  rate.  These  actual  rates  are  used  for  presentation  of  the  data.  The 
consolidate-freeze  tests  and  most  of  the  tests  on  unfrozen  sand  used  hydraulically- 
driven  loading  frames  that  also  required  several  tenths  percent  strain  before 
reaching  the  maximum  (nominal)  strain  rate. 

Three  experimental  aspects  of  the  program  of  conventional  tests  on 
frozen  sand  are  especially  noteworthy: 

•  The  use  of  enlarged,  lubricated  (via  ice  caps  and  silicone  grease)  end 
platens  to  facilitate  lateral  straining  of  the  specimen.  Although 
constrained  by  a  pin  at  each  end,  data  at  large  strains  were 
discounted  for  some  tests,  because  of  very  poor  stability  due  to 
excessive  lateral  sliding  (as  observed  after  shearing). 

•  Measurement  of  on-specimen  axial  strains  via  specially  designed 
yokes  holding  two  AC-type  linear  variable  differential  transformers 
on  opposite  sides  of  the  specimen  (see  Fig.  2.1).  These  can  resolve 
axial  strains  to  0.001%  and  eliminate  the  initially  S-shaped  stress- 
strain  curves  commonly  reported  in  the  literature  based  on 
displacement  of  the  two  end  platens  or  the  loading  piston.  Initial 
stress-strain  data  were  discounted  when  the  two  ACDT's  showed 
poor  agreement.  Note:  displacement  of  the  loading  piston  was 

used  at  >  2.5%. 

•  Measurement  of  volumetric  strains.  These  were  computed  from 
motions  of  the  cell  pressure /volume  controller  piston  and  the  axial 
loading  piston;  then  they  were  corrected  to  account  for  cell  fluid 
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leakage,  compressibility  of  the  cell  fluid  and  triaxial  chamber,  and 
flexure  of  the  base  of  the  triaxial  cell  (Andersen,  1991).  This 

technique  has  a  maximum  volumetric  strain  error  of  ±  0.2%,  which 
is  small  compared  to  the  computed  volumetric  strains  at  large  axial 
strains  and  low-to-moderate  confining  pressures. 

2.4  TESTING  PROCEDURES 

(Abstracted  from  Andersen,  1991;  Andersen  et  al,  1992,  1994;  Swan,  1994) 

2.4.1  Conventional  Tests  on  Frozen  Sand 

The  molds  used  for  preparation  and  freezing  had  been  used  by  Martin 
et  al.  (1981)  and  were  obtained  from  CRREL.  All  but  seven  specimens  were  prepared 
using  multiple  sieve  pluviation  (MSP),  which  produced  specimens  of  very  uniform 
density  (as  verified  by  radiography).  The  other  seven  used  wet  tamping,  as  did 
Martin  et  al.  (1981)  and  many  other  experimental  programs.  As  demonstrated  by 
Baker  and  Konrad  (1985),  this  technique  causes  non-uniform  densities  and 
significant  scatter  in  measured  behavior.  All  specimens  were  frozen  from  the  top 
(at  -15°C)  down  in  about  four  hours,  with  resultant  degrees  of  ice  saturation  equal  to 

99.2%  ±  1.4%  SD,  and  stored  at  -25°C  prior  to  being  carefully  trimmed  for  testing, 

usually  with  ice  end  caps  of  1-  to  3-mm  thickness.  All  specimens  had  a  prophylactic 
membrane  and  were  maintained  under  the  cell  pressure  prior  to  shearing  for  at 
least  12  hours  or  the  expected  shearing  duration,  in  order  to  allow  temperature 
stabilization  and  measurement  of  the  cell  fluid  leakage  rate.  Calculation  of  the 
specimen  dry  density  was  based  on  the  volume  of  the  trimmed  specimen  and 
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weights  measured  before  and  after  shearing.  The  reported  values  of  relative  density 
have  an  estimated  standard  deviation  of  ±  2%. 

Measurement  of  temperature  within  the  triaxial  cell  produced  the 
following  results  (±  1  std.  dev.). 

Nominal  Temp.  (°C)  Mean  Temperature  Gradient 

_ Temp.  (°C) _ TO _ 

^10  -9.55  ±  0.05  0.35  ±  0.05 

-15  -15.60  ±0.10  0.40  ±0.10 

-20  -20.50  ±  0.20  0.55  ±  0.20 

-25  -25.35  ±  0.20  0.30  ±  0.15 

The  temperature  gradient  is  thought  to  be  caused  by  heat  generated  from  the  on- 

specimen  axial  strain  device,  since  the  top  of  the  specimen  was  always  warmer  than 

the  bottom.  Note:  MIT  would  like  to  construct  a  slightly  larger  cell  with  an  internal 

mixer  to  eliminate  the  temperature  gradient  and  enable  testing  at  T  =  -5°C. 

2.4.2  Tests  on  Unfrozen  Sand 

Most  tests  were  either  isotropically  or  anisotropically  consolidated  after 
back  pressure  saturation  (usually  with  a  computer-automated  system)  and  then 
sheared  undrained  in  triaxial  compression  with  pore  pressure  measurements  at  a 

strain  rate  of  about  lO’^/s  (i.e.,  the  moderate  rate  for  frozen  sand  testing)  and  are 
designated  as  CIUC  and  CAUC  tests,  respectively.  Some  tests  were  sheared  with 
drainage  at  constant  mean  effective  stress  (CIDC  tests).  As  with  the  frozen  tests,  the 
program  employed  multiple  sieve  pluviation  (with  the  mold  sitting  on  the  base 
pedestal),  lubricated  end  platens  (usually  greased  latex  rubber)  and  on-specimen 
axial  strain  measurements. 
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2.4.3  Consolidate-Freeze  Tests 

The  consolidate-freeze  test  procedures  combined  those  of  the  unfrozen 
and  conventional  frozen  testing  programs,  but  also  required  development  of  a 
technique  for  freezing  the  specimen  in  the  triaxial  cell  while  under  a  consolidation 

stress  that  ranged  from  =  0.1  to  10  MPa.  Note;  since  the  program  wanted 
specimens  having  the  same  pre-shear  relative  density  (D^  =  95%),  the  initial  relative 
density  was  progressively  reduced  with  increasing  consolidation  stress.  After 
several  attempts,  the  freezing  process  used  the  following  procedure: 

•  saturate  and  consolidate  specimens  at  T  =  +3°C; 

•  circulate  freezing  fluid  at  T  =  -8°C  through  a  redesigned  base 
pedestal  to  cause  freezing  from  the  bottom  to  the  top  of  the 
specimen  (see  Sheet  D5  in  Appendix  D); 

•  during  freezing,  which  took  about  four  hours,  measure  the  volume 
of  water  flow  out  of  the  top  cap  via  a  drainage  line  connected  to  the 
back  pressure /volume  change  controller; 

•  also  monitor  the  temperature  via  three  thermistors  located  in  the 
cell  fluid  adjacent  to  the  bottom,  middle  and  top  of  the  specimen; 

•  note  that  the  freezing  process  ends  when  there  is  no  measured 
outflow  of  water  and/or  the  top  thermistor  registers  zero  °C;  and 

•  reduce  the  temperature  in  the  testing  room  and  the  environmental 
chamber  to  T  =  -13  and  -10°C,  respectively,  and  equilibrate  for  at 
least  12  hours  prior  to  shearing. 


17 


Although  Swan  (1994)  and  Section  5.4  report  results  from  ten 
consolidate-freeze  tests,  the  data  should  be  regarded  as  preliminary  for  two  reasons: 
first,  three  tests  developed  oil  leaks  into  the  top  drainage  line  (hence,  questionable 
estimates  of  the  volume  of  expelled  water);  and  second,  the  volume  of  expelled 
water  for  all  specimens  was  typically  less  than  85%  of  the  theoretical  volume 
increase  due  to  freezing  of  the  pore  water  (taken  as  a  9%  volume  increase). 

2.5  OVERVIEW  OF  STRESS-STRAIN  BEHAVIOR  AND  DEFINITION  OF 
STRESS-STRAIN  PARAMETERS 
2.5.1  Overview  of  Stress-Strain  Behavior 

Figure  2.2  plots  deviator  stress  (Q  =  Oi  -  o  3)  versus  engineering  axial 

strain  (Eg)  in  order  to  illustrate  the  range  in  types  of  stress-strain  curves  measured 
during  the  conventional  testing  program  on  frozen  MFS.  Two  stress-strain  curves 
are  plotted  at  each  strain  rate  and  temperature.  One  corresponds  to  specimens 
having  a  low  relative  density  (Dj,  =  35%)  and  low  confining  pressure  (o^,  =  0.1  MPa) 

and  hence  presumably  reflecting  the  minimum  frictional  resistance  of  the  sand 
skeleton.  The  other  corresponds  to  specimens  having  a  high  relative  density 
(Dr  =  95%)  and  high  confining  pressure  (a^  =  10  MPa)  and  presumably  reflecting  the 
maximum  frictional  resistance  of  the  sand  skeleton.  Figure  2.2a  shows  the  effect  of 
increasing  the  strain  rate  (e)  by  two  orders  of  magnitude  from  slow  to  fast  for 
specimeris  sheared  at  T  =  -10°C;  Figure  2.2b  shows  the  effect  of  decreasing  the 
temperature  by  10°C  for  specimens  sheared  at  the  moderate  strain  rate. 
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All  of  the  stress-strain  curves  in  Figure  2.2  exhibit  a  very  distinctive 
yield  point  (the  knee  of  the  curve)  at  levels  of  axial  strain  ranging  from  0.3  to  1.0%, 
as  shown  by  the  circles  in  the  figure.  This  point  represents  the  onset  of  highly  non¬ 
linear  behavior  and  the  development  of  very  significant  plastic  deformations;  it  is 
termed  the  upper  yield  point. 

A  qualitative  assessment  of  the  data  in  Figure  2.2  shows  that  the 
behavior  in  the  small-strain  region  is  largely  unaffected  by  changes  in  relative 
density  and  confining  pressure,  both  of  which  affect  the  frictional  resistance  of  the 
sand  skeleton.  The  results  in  Figure  2.2  also  clearly  demonstrate  that  strain  rate  and 
temperature,  which  affect  ice  behavior,  have  a  tremendous  influence  on  the 
magnitude  of  the  deviator  stress  at  the  upper  yield  point,  termed  the  upper  yield 
stress  (Quy)-  Section  3  will  evaluate  the  behavior  of  frozen  sand  in  the  small-strain 

region  where  the  measured  volumetric  strains  were  essentially  zero. 

At  larger  strains  beyond  the  upper  yield  region,  the  data  in  Figure  2.2 
show  that  changes  in  relative  density  and  confining  pressure  obviously  have  a 
tremendous  influence  on  stress-strain  behavior.  These  two  variables  also  affect  the 
frictional  resistance  of  the  sand  skeleton.  Moreover,  at  large  strains,  all  of  the  test 
specimens  exhibited  varying  degrees  of  dilation  (volumetric  expansion)  that 
depended  mainly  on  sand  density  and  confining  pressure.  Section  4  will  evaluate 
the  behavior  of  frozen  sand  in  the  large-strain  region. 
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2.5.2  Definition  of  Parameters 


Small-Strain  Behavior 

Figure  2.3  presents  actual  stress-strain  curves  drawn  to  two  different 
strain  scales  in  terms  of  the  deviator  stress  (Q  =  Oi  -  03)  and  the  axial  strain  (eg).  The 
parameters  that  have  been  selected  for  the  evaluation  of  the  small-strain  behavior 
are  the  Young’s  modulus,  the  yield  offset  at  0.01%  strain,  and  the  upper  yield  stress. 
The  graphical  construction  techniques  used  to  obtain  these  parameters  are 
siunmarized  on  the  figure  and  explained  as  follows: 

•  The  Young's  modulus  (E).  The  initial  slope  of  the  stress-strain 
curve  determined  visually  from  the  average  of  the  ACDT’s  (Fig. 

2.3a)  at  levels  of  axial  strain  less  than  0.01%. 

•  The  yield  offset  stress  at  0.01%  strain  (Qyo).  The  intersection  point  of 
the  stress-strain  curve  and  a  line  with  a  slope  of  the  Young's 
modulus  translated  by  0.01%  strain  (Fig.  2.3a). 

•  The  upper  yield  stress  (Quy)-  The  first  point  on  the  stress-strain 

curve  where  the  slope  (dQ/dfia)  either  becomes  zero  or  reaches  its 
minimum  positive  value  before  significant  strain  hardening 
(Fig.  2.3b). 

Tables  5  to  8  in  Andersen  et  al.  (1994)  and  Sheets  A1  to  A5  in  Appendix  A  present 
these  three  parameters  and  the  strain  at  the  upper  yield  stress  (Cuy)  measured  from 
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the  stress-strain  curves  for  all  tests  in  Table  2.1.  However,  values  of  E,  Qyo  and  e^y 
from  tests  with  poor  ACDT  agreement  are  excluded  from  the  subsequent  analysis  in 
Section  3. 

Large-Strain  Behavior 

The  deviator  stress  versus  axial  strain  curves  in  Figure  2.2  have 
different  characteristics  regarding  the  degree  of  post-upper  yield  strain  softening  or 

strain  hardening,  i.e.,  whether  Q  decreases  or  increases  after  reaching  Qyy;  and 
whether  or  not  strain  softening  occurs  after  the  specimen  reaches  its  peak  strength 
(Qp).  These  different  characteristics,  which  depend  mainly  on  sand  density  and 
confining  pressure,  are  illustrated  in  the  upper  portion  of  Figure  2.4  for  tests 
conducted  at  low  and  high  confining  pressures  (Oc  =  0.1  and  10  MPa).  Note  that  the 
deviator  stress  is  normalized  by  the  upper  yield  stress,  since  Q^y  does  not  vary  with 
Df  and  Figure  2.4  also  summarizes  characteristic  volumetric  behavior  measured 

at  low  and  high  confining  pressures. 

Figure  2.4  labels  four  curve  types,  and  Table  2.2  describes  their 
characteristics  and  the  typical  testing  conditions  leading  to  each  type.  In  way  of 
summary.  Type  C  curves  predominate  at  low  confining  pressures  =  0.1  MPa)  and 

are  characterized  by  post-yield  strain  hardening  to  give  a  peak  strength  (Qp)  at 
moderate  strains,  followed  by  significant  strain  softening  and  significant  dilation. 
For  loose  sand  sheared  at  higher  strain  rates  and/or  lower  temperatures,  this 
behavior  switches  to  Type  A  curves,  i.e.,  Qp  =  Q^y  followed  by  significant  strain 
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softening  and  volumetric  expansion.  In  contrast.  Type  D  curves,  predominate  at 
high  confining  pressures  =  10  MPa)  and  are  characterized  by  significant  post-yield 

strain  hardening  to  give  a  peak  strength  at  large  strains,  all  with  minimal  dilation. 
For  loose  sand  sheared  at  higher  strain  rates  and/or  lower  temperatures,  this 
behavior  switches  to  Type  B  curves,  i.e.,  initial  strain  softening  followed  by  strain 

hardening  to  produce  Qp  «  Quy/  again  with  minimal  dilation. 

The  above  qualitative  trends  are  quantified  in  Section  4  by  analyzing 
the  influence  of  the  four  testing  variables  (D^,  Gf.,  e  and  T)  on  the  following 
parameters: 

•  the  amount  of  post-upper  yield  strengthening  (AQ  =  Qp  -  Quy)  ^nd 

the  axial  strain  at  peak  strength  (Cp); 

•  the  volumetric  behavior  in  terms  of  the  maximum  rate  of  dilation 
(dev/degmax  =  maximum  slope  of  the  Cy  versus  Eg  curve)  and  the 

amount  of  dilation  at  Eg  =  20%  (Ev2o)/  ^nd 

•  the  amount  and  rate  of  post-peak  strain  softening,  both  of  which 
will  be  related  to  the  volumetric  behavior. 

Section  4  also  presents  a  conceptual  explanation  for  the  above  behavior  in  terms  of 
strengthening  mechanisms  related  to  the  frictional  resistance  of  the  sand  skeleton- 
ice  matrix  and  weakening  mechanisms  related  to  damage  (fracture)  of  the  ice  matrix. 
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Table  2,1  Scope  of  Conventional  Triaxial  Tests  on  Frozen  Manchester  Fine  Sand 
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Effect  of  Relative  Density  and  Confuting  Pressure  on  Stress- 
Strain  Curves  for  Frozen  MFS:  a)  At  Varying  Strain  Rates;  b) 
At  Varying  Temperature 
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Fig.  2.4  Types  of  Stress-Strain  Curves  Describing  Large  Strain  Behavior 

(Note:  The  magnitude  of  Q  and/or  the  rates  of  strain  hardening/ 

softening  and  dilation  vary  for  each  curve  type) 


3.  SMALL-STRAIN  BEHAVIOR  OF  FROZEN  SAND 


3.1  BACKGROUND 

The  small-strain  region  refers  to  behavior  measured  up  to  the  upper  yield 
region  (see  Fig.  2.2),  which  typically  occurs  at  axial  strain  levels  of  less  than  one 
percent.  Section  3  summarizes  data  on  Manchester  Fine  Sand  (MFS)  in  terms  of 
three  parameters  (E  =  Young's  modulus;  Qyo  =  yield  offset  stress;  and  Quy  =  upper 
yield  stress)  and  compares  these  results  with  information  reported  in  the  literature 
for  other  frozen  sands.  However,  a  review  of  the  literature  by  Andersen  et  al.  (1994) 
shows  the  following: 

•  Of  the  ten  programs  that  report  modulus  data  (tangent,  secant,  shear,  etc.), 
only  two  used  small-strain  measurement  techniques  capable  of  obtaining 
reliable  values  of  Young's  modulus,  i.e.,  within  the  linear  range  of  stress- 
strain  behavior  corresponding  to  axial  strains  less  than  0.01%. 

•  Of  the  15  programs  that  present  stress-strain  curves  with  an  upper  yield 
stress  (commonly  referred  to  as  the  initial,  first  or  lower  yield  or  the  first 
peak),  none  have  investigated  the  combined  effects  of  sand  density, 
confining  pressure,  strain  rate  and  temperature.  Moreover,  only  two 
studied  the  effect  of  sand  density,  and  only  five  evaluated  the  effect  of 
confining  pressure  at  varying  strain  rate  or  temperature. 

Sheets  Al  to  A5  of  Appendix  A  tabulate  the  specimen  properties,  testing  conditions 
and  measured  small  strain  parameters  for  all  conventional  tests  on  frozen  MFS. 
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3.2  YOUNG'S  MODULUS 


3.2.1  Data  on  Manchester  Fine  Sand 

Figure  3.1  plots  Young's  modulus  (E)  versus  relative  density  from  all 
tests  having  good  or  fair  ACDT  agreement  and  shows  the  nominal  temperature  for 
each  test.  Mean  and  standard  deviation  values  of  modulus  in  GPa  at  each 
temperature  are  26.6  ±  4.4  at -10°C;  23.9  ±  2.5  at -15°C;  28.1  ±4.0  at -20°C,  and  27.1  ± 
1.3  at  -25°C.  This  variation  is  not  statistically  significant.  Andersen  et  al.  (1994) 
conclude  that  modulus  does  not  vary  with  strain  rate,  but  does  show  a  slight 
decrease  with  higher  confinement.  The  collective  data  for  all  tests  indicate  that  E 

decreases  by  about  10  to  15%  as  Oc  increases  from  0.1  to  10  MPa.  This  is  considered  to 
be  a  second-order  effect. 

Linear  regression  on  the  data  in  Figure  3.1  predicts  a  slight  increase  in  E 

2 

with  increasing  relative  density,  but  with  very  low  correlation  coefficient  (r  =  0.06). 

Nevertheless,  the  authors  conclude  that  the  trend  is  real,  based  on  the  analysis 
presented  in  Section  3.2.3.  But  from  a  practical  viewpoint,  E  for  this  frozen  sand  is 
independent  of  the  four  testing  variables  (£,0^,  D^,  and  T)  and  equals  26.5  +  4  GPa. 

3.2.2  Comparisons  with  Prior  Data 

The  Young's  modulus  of  saturated,  frozen  sand  has  been  measured 
reliably  in  two  programs.  Kaplar  (1963)  used  the  resonant  beam  technique  to  test 
two  sands  over  a  temperature  range  of  T  =  -1.1  to  -27.8°C.  He  recorded  E  =  34.5  to 
42.4  GPa  for  Peabody  Gravelly  Sand  and  E  =  28.3  to  34.5  GPa  for  McNamara  Concrete 
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Sand.  For  both  soils,  E  increased  with  decreasing  temperature,  but  tended  to  level 
off  below  about  -10°C.  Baker  and  Kurfurst  (1985)  used  two  techniques  for  testing  16- 
100  Ottawa  Sand.  They  recorded  E  =  23  ±  9  GPa  from  on-specimen  extensometer 

tests  at  T  =  -10°C  and  E  =  27  ±  3  GPa  from  acoustic  wave  propagation  tests  at  T  =  -3.2 
and  -10°C  (no  influence  of  temperature).  In  both  test  series,  E  increased  with 
increasing  sand  density,  which  varied  between  1550  and  1780  kg/m  .  Note:  these 
authors  also  showed  a  substantial  drop  in  E  with  slight  increases  in  the  salinity  of 
the  pore  fluid,  especially  at  -3.2°C. 

The  Young’s  modulus  of  frozen  MFS  agrees  well  with  the  Kaplar  (1963) 
data  and  very  well  with  the  Baker  and  Kurfurst  (1985)  data.  Hence,  the  actual  sand 
type  has  a  relatively  minor  effect  on  Young's  modulus.  In  contrast,  values  of 
"initial  modulus"  reported  in  the  literature  from  numerous  programs  that  did  not 
use  on-specimen  strain  devices  or  acoustic  wave  propagation  techniques  are  much 
lower,  e.g.,  typically  in  the  range  of  only  0.1  to  8  GPa.  These  modulus 
determinations  were  made  at  strain  levels  much  higher  than  the  Eg  <  0.01%  used  for 
frozen  MFS.  These  collective  data  also  indicate  a  tendency  for  modulus  to  increase 
with  decreasing  confining  pressure,  increasing  strain  rate  and  decreasing 
temperature.  As  shown  in  Section  3.3,  these  trends  are  consistent  with  those 
observed  for  the  yield  offset  stress  (Qyo)  of  frozen  MFS;  therefore,  they  undoubtedly 
include  some  nonlinear  plastic  deformations. 
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3.2.3  Prediction  of  Young's  Modulus 

The  Young's  modulus  of  frozen  MFS  is  about  three  times  larger  than 
that  reported  for  polycrystalline  ice  (i.e.,  9  GPa  by  Sinha  1989)  and  one  to  two  orders 
of  magnitude  larger  than  measured  in  consolidated-undrained  triaxial  compression 
tests  on  unfrozen  MFS  having  consolidation  stresses  ranging  from  0.1  to  10  MPa 
(i.e.,  0.2  to  4  GPa  by  Swan  1994).  A  logical  explanation  for  why  E  of  frozen  sand  is 
larger  than  either  of  its  components  is  that  frozen  sand  acts  as  a  composite  material 
consisting  of  a  stiff  aggregate  (sand  particles)  within  a  soft  matrix  (ice).  That  is,  it  is 
the  stiffness  of  the  sand  particles,  and  not  that  of  the  sand  skeleton,  that  dominates 
the  overall  stiffness  of  the  frozen  sand  system  at  small  strains. 

Andersen  (1991)  and  Andersen  et  al  (1992,  1994)  applied  Counto's 
(1964)  isostrain  composite  material  model,  which  was  developed  for  concrete,  to 
frozen  MFS.  Flis  model  can  be  visualized  as  a  prism  of  a  solid  aggregate  material 
(silicate  mineral  representing  the  sand  phase)  embedded  in  a  cube  of  matrix  material 
(poly crystalline  ice  representing  the  ice  phase),  with  perfect  bonding  between  the 
sand  particles  and  the  pore  ice.  The  resulting  equation  for  the  composite  Young's 
modulus  is 
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where  E^.  is  the  composite  modulus;  Ej  and  Eg  are  the  Young's  moduli  for  the  ice 


matrix  and  the  sand  grains,  respectively;  and  Vg  is  the  volume  fraction  of  the  sand 


grains. 
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Using  a  Young's  modulus  for  sand  particles  of  90  GPa  (from  Table  12.5  of 
Lambe  and  Whitman  1969,  assuming  an  average  value  for  Quartzite)  and  9  GPa  for 
polycrystalline  ice  from  Sinha  (1989),  the  estimated  composite  Young  s  modulus  for 

frozen  MFS  ranges  from  24.0  to  29.6  GPa  for  volume  fractions  ranging  from  0.52  (D^ 
=  0%)  to  0.63  (Dr  =  100%).  Figure  3.1  compares  the  predicted  and  measured  Yoimg’s 

modulus  for  frozen  MFS  having  relative  densities  varying  from  20  to  100%. 

The  excellent  agreement  between  the  predicted  and  measured  Young's 
modulus  in  Figure  3.1  and  the  agreement  with  prior  reliable  measurements 
indicates  that  the  Counto  isostrain  model  can  offer  insight  as  to  why  the  Young's 
modulus  of  frozen  sand  is  so  much  greater  than  either  of  its  two  components.  It 
also  helps  to  explain  why  the  type  of  sand  has  little  effect  on  Young  s  modulus  of  a 
frozen  cohesionless  soil  (Section  3.2.2).  The  silicate/ice  interface  of  frozen  MFS  in 
the  range  of  temperatures  and  strain  rates  studied  herein  appears  capable  of 
transmitting  the  shear  stresses  necessary  to  allow  the  silicate  particles  to  reinforce 
the  ice  matrix  in  composite  action,  thus  causing  the  frozen  sand  to  behave  as  a 
composite  material  at  very  small  strains.  The  mineralogy  of  the  silicate  particles 
(governing  their  elastic  properties)  is  much  more  important  to  the  initial  stress- 
strain  response  than  is  the  structure  of  the  sand  skeleton.  Additionally,  the  dry 
density  of  the  sand  is  important  insofar  as  it  affects  the  volume  fraction  of  the 
silicate  minerals  in  the  frozen  composite  material. 
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3.3  YIELD  OFFSET  STRESS 

The  yield  offset  stress  is  considered  in  order  to  include  a  parameter  lying 
between  the  truly  elastic  region  (Young's  modulus)  and  the  onset  of  highly  non¬ 
linear  and  plastic  behavior  at  the  upper  yield  stress.  Figure  3.2  plots  the  yield  offset 

stress  (Qyo)  versus  temperature  as  a  function  of  strain  rate  for  the  two  extremes  of 
confining  pressure  (o^.  =  0.1  and  10  MPa).  Analysis  of  the  extensive  data  at  T  =  -10  C 
shows  no  effect  of  sand  density.  For  the  other  three  variables,  the  results  in  Figure 
3.2  and  similar  analysis  for  Qyo  versus  strain  rate  indicate  the  following; 

•  Qyo  is  strongly  dependent  on  temperature  (increasing  with  decreasing 
temperature),  with  dQyo/dT  (MPa/°C)  decreasing  consistently  from  about  - 
0.02  at  the  slow  strain  rate  to  about  -0.25  at  the  fast  strain  rate. 

•  Qyo  is  also  strongly  dependent  on  initial  secant  strain  rate.  Using  a  power 
law  relationship  of  the  form  (Glenn  1955,  Sayles  1973) 


1 


where  B  is  a  temperature-dependent  constant,  the  power  law  coefficient 
(n)  has  values  of  13.2  ±  5.2  (and  r^  =  0.49  ±  0.19)  with  no  consistent 
variation  as  a  function  of  temperature  and  confining  pressure. 

•  Qyo  decreases  with  increasing  confining  pressure.  Although  this  pressure 
sensitivity  cannot  be  defined  as  a  function  of  strain  rate  and  temperature. 
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the  collective  data  suggests  that  the  dQyo/ da^  =  -0.1  ±  0.05  and  perhaps 
tending  towards  zero  at  the  lowest  temperature. 

3.4  UPPER  YIELD  STRESS 

3.4.1  Data  on  Manchester  Fine  Sand 

Definition  of  the  upper  yield  stress  (Quy)  as  a  fimction  of  the  testing 
variables  is  considered  to  be  one  of  the  most  important  contributions  of  this 
research.  Figure  3.3  illustrates  the  fact  that  relative  density  does  not  affect  Quy  for 
specimens  sheared  at  varying  strain  rates  and  T  =  -10°C.  For  tests  sheared  at  the 
moderate  strain  rate  and  a  temperature  of  -10°C,  the  data  in  Figure  3.3  and  the 
tabulated  results  in  Appendix  A  show  a  fairly  consistent  decrease  in  Quy  with 

increasing  confining  pressure,  i.e.,  dQuy/da,.  =  -0.06  with  r  =  0.2.  However,  the 
pressure  sensitivity  of  the  upper  yield  stress  is  considered  a  second  order  effect  since 
the  collective  data  at  the  other  strain  rates  and  lower  temperatures  do  not  show  the 

same  consistent  change  in  Quy  with  confining  pressure.  Therefore,  Figures  3.4  and 

3.5  can  be  used  to  define  Quy  as  a  function  of  strain  rate  and  temperature, 
respectively,  independent  of  sand  density  and  confining  pressure.  The  strain  rate 
used  in  Figure  3.4  is  the  tangent  (instantaneous)  strain  rate  at  the  upper  yield  stress, 

which  typically  was  10  ±  5%  less  than  the  nominal  (maximum)  strain  rate. 

For  the  data  in  Figure  3.4,  a  power  law  relationship  for  the  strain  rate 
dependency  can  be  used  of  the  form 
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Q.,=a(£)"  . P.31 

where  A  is  a  constant  for  a  given  temperature.  For  the  data  in  Figure  3.5,  the 
temperature  dependency  can  be  expressed  using  a  linear  relationship  of  the  form 

O  =C  +  DT  . [3.4] 

>cuy 

where  C  and  D  are  constants  for  a  given  strain  rate. 

Sheet  A6  presents  the  results  of  linear  regression  analyses  on  the  data, 

which  show  the  following  trends: 

•  For  the  data  in  Figure  3.4,  n  increases  in  a  very  consistent  fashion 

from  about  4.6  to  6.6  as  the  temperature  decreases  from  -10  to  -25°C. 

That  is,  the  upper  yield  stress  becomes  less  strain  rate  dependent  at 
lower  temperatures. 

•  For  the  data  in  Figure  3.5,  D  =  dQ^y/dT  varies  in  a  very  consistent 

fashion  from  -0.58  at  the  slow  strain  rate  to  -0.97  at  the  fast  strain 
rate.  That  is,  the  upper  yield  stress  becomes  more  temperature 
dependent  with  increasing  strain  rate. 

Appendix  A  lists  the  values  of  axial  strain  (Euy)  at  the  upper  yield  stress, 
which  range  from  0.25  to  1.1%  over  all  testing  conditions.  For  shearing  at  the  slow 
and  moderate  strain  rates,  mean  values  of  e^y  and  Qyy  at  each  e  and  T  give  the 

relationship  (r  =  0.85) 
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£uy(%)  =  0.12  +  0.04  Quy(MPa) 
where  the  mean  Quy  ranged  from  5.1  to  20.7  MPa. 


[3.5] 


3.4.2  Comparisons  with  Prior  Data 

Baker  and  Konrad  (1985)  represent  the  only  other  program  to  measure 

upper  yield  stress  versus  sand  density.  Their  results  show  no  change  in  Quy  for 

pluviated  30-100  Ottawa  sand  with  ranging  from  30  to  100%  and  tested  in 
unconfined  compression  with  on-specimen  axial  strain  measurements  at 
c  =  1.7  X  10"^/sec  and  T  =  -10°C.  This  behavior  is  consistent  with  the  frozen  MFS 

data  in  Figure  3.3.  Their  mean  value  of  Quy  also  is  within  about  15%  of  that 

predicted  for  MFS  at  the  same  strain  rate  and  temperature. 

The  primary  focus  of  other  programs  reported  in  the  literature  was  on 
peak  strength;  they  contained  relatively  little  discussion  of  upper  yield  behavior. 
The  fact  that  these  programs  did  not  employ  on-specimen  axial  strain 
measurements  and  usually  tested  "compacted"  specimens  also  prevented  reliable 
detection  of  the  upper  yield  stress  in  a  consistent  fashion.  These  limitations  should 
be  born  in  mind  when  comparing  prior  data  with  results  from  this  research. 

Figure  3.6  plots  the  upper  yield  stress  as  a  function  of  confining 
pressure  for  frozen  Ottawa  sands  tested  by  Sayles  (1973),  Chamberlain  et  al.  (1972) 
and  Parameswaran  and  Jones  (1981).  The  mesh  sizes  for  the  three  sands  were  20-30, 
100-200  and  30-100,  respectively.  The  Sayles  tests  were  conducted  at  -3.85°C;  the 
others  were  conducted  at  T  =  -10°C.  The  current  data  at  the  moderate  strain  rate  for 
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all  relative  densities  and  a  temperature  of  -10°C  are  also  plotted.  Although  the 
Parameswaran  and  Jones  (1981)  test  program  produced  a  higher  upper  yield  stress 

than  expected  for  Manchester  Fine  Sand  at  the  same  strain  rate  (i.e.,  Qyy  =  10  MPa 
from  Figure  3.4),  the  modest  pressure  sensitivity  is  consistent  with  the  collective 
data  on  MFS  which  showed  dQuy/dOc  =  0  ±  0.1  from  nine  sets  of  data  (albeit  that  the 

Parameswaran  and  Jones  data  exhibit  a  consistent  increase  with  confinement).  Note 
that  the  confining  pressure  in  the  Parameswaran  and  Jones  tests  was  not  constant, 
but  increased  during  each  test  (the  confining  pressure  used  in  Figure  3.6  is  that 
reported  at  the  upper  yield  stress). 

Two  of  the  three  tests  by  Chamberlain  et  al.  (1972)  produced  values  of 

significantly  lower  than  expected  for  Manchester  Fine  Sand  at  the  same  strain  rate 

(i.e.,  18  MPa  for  Figure  3.4).  Both  the  Sayles  (1973)  data  at  lower  confining  pressures 

and  the  Chamberlain  et  al.  (1972)  data  have  much  larger  pressure  sensitivities  than 

observed  for  MFS.  Perhaps  coarse-grained  sands  and  higher  strain  rates  cause  earlier 

fracturing  of  the  ice  matrix,  leading  to  a  pressure  sensitivity  similar  to  that  reported 

by  Jones  (1978, 1982)  for  polycrystalline  ice.  At  higher  confining  pressures,  the  Sayles 

data  indicated  less  pressure  sensitivity  in  a  manner  consistent  with  the 

Parameswaran  and  Jones  (1981)  results  and  those  in  the  current  testing  program. 

-7 

Parameswaran  and  Roy  (1982)  present  strain  rate  data  (e  =  5  x  10  to 
6  X  10'^/ sec)  on  upper  yield  stress  for  a  compacted,  saturated  Ottawa  30-100  sand  in 

unconfined  compression  at  -30°C.  They  reported  a  power  law  coefficient  of  10.6, 
which  means  a  lower  rate  sensitivity  than  for  frozen  MFS.  However,  a  direct 
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comparison  with  the  power  law  coefficients  of  the  present  data  may  not  be  valid, 
since  they  did  not  have  the  true  strain  rate  (and  their  irutial  stress-strain  curves 
were  S-shaped)  and  they  used  compacted  rather  than  pluviated  specimens. 

3.4.3  Discussion  of  Mechanisms  and  Analogy  with  Peak  Ice  Strength 

Various  researchers  have  presented  short  discussions  of  possible 
mechanisms  that  may  operate  within  the  upper  yield  region  during  the  shearing  of 
frozen  soils.  These  are  focused  on  drawing  an  analogy  between  the  upper  yield 
behavior  of  frozen  soil  and  the  peak  strength  behavior  of  polycrystalline  ice,  e.g., 
Chamberlain  et  al  (1972),  Bragg  and  Andersland  (1980),  Parameswaran  and  Jones 
(1981),  and  Sayles  (1989),  among  others.  Brief  excerpts  from  these  discussions  follow: 

Chamberlain  et  al  (1972)  "The  distinct  yield  points  in  the  stress- 
strain  curves  for  the  OWS  'Ottawa  sand’  samples  tested  in  the  range  of 
Gc  =  0.5  -  5.0  kip/in^  are  most  probably  related  to  breaking  of  the  ice-sand  bonds  and 
shearing  of  ice  crystals." 

Bragg  and  Andersland  (1980)  "The  initial  yield  occurs  as  the  pore  ice 
reaches  its  final  yield  stress  and  as  interparticle  friction  and  dilatancy  begin  to 
contribute  to  the  shear  strength.... With  increasing  strain  rates  the  ice  matrix  strength 
increases  giving  a  higher  initial  yield  stress." 

Parameswaran  and  Jones  (1981)  "This  indicates  that  deformation 
behavior  in  the  early  elastic  region  is  governed  mainly  by  the  ice  matrix.. ..The  larger 
value  of  the  yield  strength  for  frozen  sand  as  compared  to  ice  is  due  to  the 
mobilization  of  soil  grain  friction,  in  addition  to  cohesion  of  ice.  As  in  ice,  an  initial 
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increase  in  yield  stress  with  pressure  can  be  explained  by  the  closure  of  voids  in  the 
sand-ice  specimens,  thereby  allowing  the  ice  matrix  to  deform  plastically." 

Sayles  (1989)  "It  should  be  noticed  that  yield  (initial  sharp  bending  of 
the  curve)  occurs  at  less  than  1%  strain,  even  though  the  maximum  stress  is  not 
reached.  There  is  evidence  to  suggest  that  initial  cracking  of  the  ice  matrix  occurs  at 
this  point....That  the  initial  yield  represents  the  domination  of  the  ice  matrix 
cohesive  component... is  reinforced  by  the  shape  of  the  Mohr  s  strength 
envelope....The  envelope  for  the  initial  yield  (i.e.,  first  peak)...is  curved  at  lower 
values  of  confining  pressures,  indicating  that  some  internal  friction  is  involved  at 
these  lower  levels  of  confining  pressures.  This  'first  peak'  envelope  approaches  a 
constant  value  at  the  higher  confining  pressures,  suggesting  that  ice  cohesion 
provides  nearly  all  the  initial  yield  resistance  at  these  confining  pressures  and  small 
strains." 

Thus,  there  is  general  consensus  that  the  upper  yield  behavior  of 

frozen  sand  is  dominated  by  the  ice  matrix. 

Figure  3.7  compares  the  upper  yield  stress  as  a  function  of  strain  rate  for 
frozen  MFS  (from  Fig.  3.4)  and  the  peak  strength  of  polycrystalline  ice  at 
temperatures  ranging  from  -7  to  -12°C,  confining  pressures  up  to  10  MPa  and  grain 
sizes  of  0.7  and  1.0  mm.  These  ice  data  were  reported  by  Flawkes  and  Mellor  (1972) 
and  Jones  (1982).  The  power  law  coefficients  (n)  for  ice  range  from  4.3  to  5.4  and 
correspond  well  to  those  of  the  Q^y  of  frozen  MFS.  This  suggests  that  the  behavior 

of  the  frozen  sand  in  the  upper  yield  region  is  qualitatively  similar  to  that  of 
polycrystalline  ice.  However,  the  magnitude  of  the  upper  yield  stress  at  comparable 


temperatures  is  approximately  double  that  of  the  strength  of  polycrystalline  ice  in 
unconfined  compression.  Apparently,  the  presence  of  the  sand  particles  causes  a 
strengthening  effect,  even  though  the  frictional  characteristics  of  the  sand  skeleton 
seem  to  be  unimportant  as  evidenced  by  the  insensitivity  to  changes  in  D,.  and  0^ 
(Fig.  3.3).  It  is  also  important  to  note  that  the  ice  peak  strength  data  of  Jones  (1982) 
show  a  confining  pressure  sensitivity  unlike  the  upper  yield  stress  of  frozen 
Manchester  Fine  Sand. 

The  strong  dependence  of  the  upper  yield  behavior  on  the  applied 
strain  rate  and  temperature,  its  independence  of  the  relative  density  and  its  very 
slight  dependence  on  the  confining  pressure  (considered  a  second-order  effect)  all 
indicate  that  the  structure  of  the  sand  skeleton  (i.e.,  its  frictional  resistance  and 
tendency  to  dilate  or  contract)  probably  does  not  play  a  significant  role  in  controlling 
the  behavior  of  frozen  sand  in  the  upper  yield  region.  On  the  other  hand,  the  close 
agreement  between  the  power  law  coefficients  for  the  peak  strength  of 
polycrystalline  ice  and  the  upper  yield  stress  of  frozen  sand  over  a  wide  range  in 
strain  rates  and  temperatures  suggests  that  the  principal  physical  mechanisms 
occurring  within  the  upper  yield  region  are  probably  similar  to  those  controlling  the 

behavior  of  polycrystalline  ice  at  its  peak  strength. 

Unfortunately,  the  ice  data  in  Figure  3.7  are  at  strain  rates  beyond  the 
ductile  (i.e.,  flow  or  creep)  range  wherein  the  effects  of  strain  rate  and  temperature 
can  be  modeled  by  a  combination  of  the  Norton-type  power  law  and  the  Arrhenius 
temperature  equation  according  to  Glen(1955),  i.e.,  the  constant  power  law 
coefficient  and  activation  energy.  Consequently,  an  understanding  of  the  physical 
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mechanisms  controlling  the  upper  yield  behavior  of  frozen  sand  (in  the 
deformational  regime  investigated  in  this  research)  will  first  require  an  improved 
knowledge  of  the  physical  mechanisms  controlling  the  behavior  of  ice  within  the 
ductile-to-brittle  transition  region. 

It  is  interesting  to  note  that  the  initial  shapes  of  the  stress-strain  curves 
presented  in  Figure  2.2  are  qualitatively  similar  to  those  reported  by  Baker  et  ah 
(1982),  Chamberlain  et  ah  (1972),  and  Parameswaran  and  Jones  (1981).  Additionally, 
they  are  remarkably  similar  to  the  stress-strain  curves  presented  by  Mellor  and  Cole 
(1982)  for  polycrystalline  ice. 

The  authors  concur  with  prior  researchers  that  the  upper  yield 
behavior  is  attributable  to  the  pore  ice,  but  have  seen  no  evidence  of  a  frictional 
contribution  as  suggested  by  Sayles  (1973,  1989)  and  Parameswaran  and  Jones  (1981) 
for  frozen  Manchester  Fine  Sand  in  the  ranges  of  temperature  and  strain  rate 
investigated  in  this  program. 

3.5  SUMMARY  AND  CONCLUSIONS 

The  measured  small-strain  response  of  frozen  Manchester  Fine  Sand  (MFS)  is 
reported  in  terms  of  the  Young's  modulus,  the  yield  offset  stress  at  0.01%  strain,  and 
the  upper  yield  stress  (first  point  on  the  stress-strain  curve  where  the  slope  either 
becomes  zero  or  reaches  its  minimum  positive  value  before  significant  strain 
hardening).  Figures  2.2  and  2.3  present  stress-strain  curves  that  show  the  upper 
yield  region. 
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Figure  3.1  demonstrates  that  the  Young's  modulus  (E)  is  nearly  independent 
of  testing  variables.  There  is  a  slight  (10-15%)  decrease  in  E  with  increasing  from 
0.1  to  10  MPa  (considered  a  second-order  effect),  and  a  slight  increase  (15%)  with  D, 
increasing  from  20  to  100%.  The  measured  modulus  equals  26.5  ±  4.0  GPa  over  all 
the  testing  variables.  This  value  agrees  well  with  reliable  measurements  on  three 
other  frozen  sands. 

The  mineralogy  of  the  sand  grains,  the  sand  density,  and  the  ice/ silicate  bond 
strength  are  considered  to  be  the  most  significant  factors  that  contribute  to  the  initial 
stiffness  properties  of  frozen  MFS  in  the  ranges  of  temperature  and  strain  rate 
investigated  in  this  program.  A  composite  materials  model  developed  by  Counto 
(1964)  for  concrete  has  been  extended  to  frozen  sands.  The  excellent  agreement 
between  predicted  and  measured  data  in  Figure  3.1  indicates  that  composite  action 
between  the  pore  ice  and  silicate  particles  is  the  mechanism  responsible  for  a 
Young's  modulus  that  is  several  times  greater  than  that  of  polycrystalline  ice  and  up 
to  two  orders  of  magnitude  greater  than  that  of  the  sand  skeleton. 

An  upper  yield  stress  (Quy)  was  identified  from  stress-strain  curves  such  as 
illustrated  in  Figure  2.3  for  a  range  in  testing  conditions.  The  measured  volumetric 
strains  up  to  the  upper  yield  stress  were  essentially  zero.  The  value  of  the  upper 
yield  stress  does  not  depend  on  relative  density  and  is  only  slightly  affected  by  the 
magnitude  of  the  confining  pressure,  as  illustrated  in  Figure  3.3.  This  indicates  that 
the  frictional  characteristics  of  the  sand  skeleton  play  a  minor  role,  since  these  two 
variables  have  a  strong  influence  on  the  stress-strain  behavior  of  unfrozen  sand. 
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However,  the  upper  yield  stress  is  very  dependent  on  the  strain  rate  and  the 
temperature,  as  demonstrated  in  Figures  3.4  and  3.5,  respectively.  Although  the 
upper  yield  stress  of  frozen  sand  is  larger  than  the  peak  strength  of  polycrystalline 
ice  (Fig.  3.7),  both  materials  have  similar  power  law  coefficients,  i.e.,  in  the  range  of 

5.5  ±  1.0.  This  suggests  that  the  principal  physical  mechanisms  occurring  in  the 

upper  yield  region  are  probably  similar  to  those  controlling  the  behavior  of 
polycrystalline  ice.  This  is  in  general  consensus  with  the  opinions  of  other 
researchers.  However,  the  actual  nature  of  the  mechanisms  controlling  the  upper 
yield  response  are  still  poorly  understood. 

The  yield  offset  stress  at  0.01%  axial  strain  (Qyo)  is  presented  in  order  to  bridge 
the  gap  between  the  initial  modulus  and  the  upper  yield  behavior.  The  relative 
density  has  no  influence,  and  there  is  a  modest  decrease  in  Qy^,  with  increasing 
confining  pressure.  It  is  strongly  influenced  by  strain  rate  and  temperature,  but  to  a 
lesser  degree  than  the  upper  yield  stress  (Fig.  3.2).  This  indicates  that  the  behavior  of 
frozen  sand  in  the  small-strain  region  becomes  more  strain  rate  and  temperature 
sensitive  as  the  level  of  strain  increases. 

These  data  present  the  most  comprehensive  treatment  of  the  small-strain 
behavior  of  a  frozen  sand  available  in  the  literature;  we  hope  they  will  assist  in  the 
development  of  improved  modeling  capabilities  for  describing  the  small-strain 
behavior  of  frozen  sands. 
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4.  LARGE-STRAIN  BEHAVIOR  OF  FROZEN  SAND 


4.1  BACKGROUND 

Table  B  (Sheets  B1  to  B4  in  Appendix  B)  summarizes  experimental  programs 
on  saturated  (or  nearly  so)  frozen  sands  that  report  large-strain  behavior.  For  each 
program,  the  table  lists  the  reference  and  the  following  information: 

•  name  or  type  of  sand  and  gradation  characteristics,  the  method  of 
specimen  preparation  (C  =  compacted,  P  =  pluviation,  V  =  vibration  and 
MSP  =  multi-sieve  pluviation)  and  the  degree  of  ice  satmation; 

•  whether  or  not  the  program  used  an  on-specimen  axial  strain  device  and 
measured  volumetric  strains; 

•  the  range  of  testing  variables,  i.e.,  sand  density,  confining  pressure  (Oc), 

axial  strain  rate  (e)  and  temperature;  and 

•  the  principal  focus  of  the  programs. 

The  programs  listed  in  Table  B  have  the  following  attributes  and  limitations: 

•  about  three-quarters  tested  Ottawa  sands,  the  most  prevalent  being 
uniform  coarse  (20-30  sieve)  and  uniform  medium  (=  20-100  sieve); 

•  about  one-third  varied  confining  pressure  and  one-third  varied  strain  rate 
and  temperature  (with  =  0),  while  only  one  program  (I)  studied  these 

three  variables  in  some  detail; 

•  only  two  included  a  detailed  study  of  relative  density; 
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•  about  half  prepared  their  specimens  via  compaction  (rather  than 
pluviation),  which  leads  to  non-uniform  densities  and  scattered  results 
(Baker  and  Konrad  1985);  and 

•  four  made  measurements  of  volumetric  strains,  and  only  one  measured 
on-specimen  axial  strains. 

Therefore,  the  program  on  frozen  Manchester  Fine  Sand  (MFS)  is  umque  by 
measuring  volumetric  and  on-specimen  axial  strains  on  pluviated  specimens  with 
lubricated  ends;  the  specimens  had  a  wide  range  of  relative  densities  (Dj)  and  were 
sheared  at  varying  confining  pressures  (Oj.)/  strain  rates  (e),  and  temperatures  (T),  as 
summarized  in  Table  2.1.  The  only  known  limitation  involves  compliance  of  the 
loading  system  that  required  about  0.5  to  1%  axial  strain  before  achieving  the 
nominal  (maximum)  strain  rate.  However,  the  authors  believe  that  the  prior 
slower  strain  rate  should  not  have  a  significant  effect  on  the  measured  large-strain 
behavior,  since  Mellor  and  Cole  (1983)  show  that  ice  has  a  direct  correspondence 
between  constant  strain  rate  and  constant  stress  creep  tests.  On  the  other  hand,  there 
is  some  error  in  reporting  axial  stresses  normalized  by  the  upper  yield  stress  (i.e.,  use 

of  the  stress  ratio  =  Q/Quy),  since  Qyy  usually  occurred  at  an  instantaneous  strain  rate 
about  10  ±  5%  lower  than  the  maximum  strain  rate.  For  typical  power  law 

coefficients  of  n  =  10  ±  5,  this  corresponds  to  an  error  in  Q/Quy  ranging  from  0.5  to 
3%.  Given  the  inherent  scatter  in  the  data,  no  correction  was  made  for  this  small 


error. 


Sheets  B5  to  Bll  present  tabulated  data  from  the  conventional  triaxial 
compression  tests  rim  on  frozen  MFS. 

4.2  OVERVIEW  OF  STRESS-STRAIN  BEHAVIOR  AT  LARGE  STRAINS 

Section  3.4  concluded  that  the  upper  yield  stress  (Quy)  of  frozen  MFS  varies 
only  with  strain  rate  and  temperature,  and  that  this  behavior  is  consistent  with  the 
general  behavior  of  poly  crystalline  ice.  But  at  large  strains,  changes  in  sand  density 
and  confining  pressure  have  a  very  significant  effect  on  both  the  deviator  stress 

(Q  =  Oj  -  03)  and  the  volumetric  strain  (Ey).  Section  2.5.2  divided  the  large-strain 

behavior  into  four  basic  types  of  behavior.  These  are  designated  as  curve  types  A,  B, 
C  and  D,  as  described  in  Table  2.2  and  Figure  2.4.  Note  that  the  stress-strain  curves 

in  Figure  2.4  are  normalized  by  using  the  stress  ratio,  Q/ Quy. 

Figures  4.1  through  4.3  plot  actual  data  in  order  to  support  the  basis  for  the 
four  curve  types  and  to  illustrate  how  changes  in  density  and  confinement  affect  the 
large-strain  behavior.  These  changes  can  be  summarized  as  follows. 

•  Figure  4.1  shows  that  increasing  relative  density  (D,)  at  low  confinement 

(Oc  =  0.1  MPa)  initially  changes  the  post-upper  yield  behavior  from  strain 
softening  to  strain  hardening,  causes  an  increase  in  the  peak  strength  (Qp) 

that  is  followed  by  more  pronounced  post-peak  strain  softening,  and 
causes  an  increase  in  the  amount  of  volumetric  expansion  (dilation).  The 
curve  designation  changes  from  Type  A  or  B  at  low  density  to  Type  C  at 
higher  densities. 
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•  Figure  4.2  shows  that  increasing  confinement  during  shear  of  dense  sand 
(Dj.  =  90%)  causes  an  increase  in  the  degree  of  post-upper  yield  strain 

hardening,  which  produces  higher  peak  strengths  at  larger  axial  strains. 
This  behavior  is  accompanied  by  less  post-peak  strain  softening  and  a  large 
decrease  in  the  amount  of  dilation.  The  curve  designation  changes  from 
Type  C  to  Type  D. 

•  Figure  4.3  shows  that  increasing  relative  density  at  high  confinement 

(Oc  =  10  MPa)  causes  more  pronounced  Type  D  behavior  for  shearing  at  the 

moderate  strain  rate  and  T  =  -10°C;  it  causes  a  change  from  Type  B  to  Type 
C  behavior  for  shearing  at  the  fast  strain  rate  and  T  =  -15°C. 

Section  2.5.2  described  the  parameters  that  will  be  used  to  quantify  the 

influence  of  the  four  testing  variables  (Dj,  Cf.,  e  and  T)  on  large-strain  behavior. 

These  parameters  relate  to  the  amount  of  post-upper  yield  strengthening 

(AQ  =  Qp  -  Quy)  and  the  peak  strength  (Qp),  the  volumetric  behavior  in  terms  of  the 

maximum  rate  of  dilation  and  the  amount  of  dilation  at  20%  axial  strain,  and 
various  approaches  for  relating  the  amount  of  post-peak  strain  softening  to  the 
volumetric  behavior.  These  aspects  are  treated  in  Sections  4.3,  4.4  and  4.5, 
respectively.  Section  4.6  presents  a  conceptual  hypothesis  to  explain  the  large-strain 
behavior  in  terms  of  strengthening  mechanisms  related  to  the  frictional  resistance 
of  the  sand  skeleton-ice  matrix  and  weakening  mechamisms  related  to  damage 
(fracture)  of  the  ice  matrix. 
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4.3  POST-UPPER  YIELD  STRENGTHENING  AND  PEAK  STRENGTH 
4.3.1  Comparison  of  Peak  Strengths  with  Prior  Data 

Figure  4.4  compares  the  peak  strength  of  frozen  MFS  with  data  on 
other  frozen  sands  as  a  function  of  strain  rate  at  temperatures  varying  from  T  =  -10° 
to  -30°C.  All  prior  programs  that  investigated  the  effects  of  strain  rate  and 
temperature  we  performed  on  dense  to  medium-dense  sand,  and  most  used 

unconfined  compression  tests  (Oj.  =  0).  Figure  4.4a  compares  results  (from  linear 

regression)  at  low  confinement  (Oc  <0.1  MPa)  along  with  the  corresponding  values 

of  the  power  law  coefficient  (n).  The  location  of  the  MFS  and  its  values  of  n  are 
representative  of  frozen  sand  data  with  one  exception:  Bragg  and  Andersland  (1980) 

show  no  strength  increases  beyond  e  «  10*^/sec  (also  true  for  their  tests  at  T  =  -6  and 
-15°C).  Also  note  the  two-  to  three-fold  variation  in  Qp  at  the  same  strain  rate  for 

sands  tested  at  T  =  -10°C.  This  variation  probably  reflects  differences  in  relative 
density,  grain  size  distribution,  degree  of  ice  saturation,  and  the  method  of  specimen 
preparation  (compaction  versus  pluviation). 

Figure  4.4b  compares  peak  strength  versus  strain  rate  for  the  only  other 

program  with  similar  testing  at  high  confinement  (a^  =  10  MPa).  The  T=-10°C 
Shibata  et  al.  (1985)  data  on  dense,  uniform  fine  Toyoura  sand  (that  may  have  been 
pluviated)  lies  between  the  MFS  results  at  T  =  -10  and  -15°C.  Note:  the  capacity  of 
MIT's  load  frame  prevented  shearing  at  the  fast  strain  rate  at  lower  temperatures. 

Figure  3.4  showed  that  the  upper  yield  stress  of  frozen  MFS  became  less 
rate  sensitive  (i.e.,  n  increasing  from  4.6  to  6.6)  at  lower  temperatures.  The  collective 
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results  in  Figure  4.4  show  that  the  peak  strength  of  dense  sands  also  tends  to  become 
less  rate  sensitive  at  lower  temperatures.  However,  the  rate  sensitivity  of  the 
strength  of  dense  frozen  sands  is  also  much  lower  than  that  for  the  upper  yield 
stress,  i.e.,  higher  n  values. 

4.3.2  Effects  of  Relative  Density  and  Confining  Pressure 

Effect  of  Relative  Density 

Figure  4.5  plots  peak  strength  (Qp)  versus  relative  density  for  tests  on 
frozen  MFS  at  different  temperatures  at  the  moderate  strain  rate  and  low 
confinement  (a^  =  0.1  MPa).  All  tests  at  T  =  -10  and  -15°C  had  Type  C  curves  and 

produce  a  well-defined,  linear  increase  in  Qp  with  increasing  Dj.  However,  the 

looser  specimens  at  lower  temperatures  had  either  Type  A  curves  (Qp  =  Quy)  or  Type 

C  curves  with  Qp  only  sUghtly  larger  than  Q^y.  These  data  sets  produce  a  bilinear  Qp 

versus  D,  relationship  where  definition  of  the  slope  at  higher  densities  can  be 

tmcertain.  This  is  illustrated  by  the  results  at  T  =  -20°C,  where  the  tests  at  ==  45  to 

50%  had  lower  Q„v  values  than  the  test  at  =  36%.  The  solid  and  dashed  lines 

represent  two  possible  interpretations  of  this  data  set. 

In  an  attempt  to  minimize  this  problem,  which  basically  arises  from 

scatter  in  the  Q^y  data  (remember  that  Q^y  should  not  vary  with  D,  or  as  per 
Section  3.4),  the  influence  of  density  and  confinement  on  the  peak  strength  is 
evaluated  in  terms  of  the  strength  increase  relative  to  the  upper  yield  stress.  In 
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other  words,  the  results  are  also  analyzed  using  AQ  =  Qp  -  Quy/  where  Q^y  is  the 
value  measured  for  each  individual  test. 

Figure  4.6  plots  AQ  versus  for  tests  with  low  and  high  confinement 
at  all  strain  rates  and  temperatures.  Note  that  results  from  tests  with  Type  A  curves 
are  plotted  below  the  X  axis,  with  the  label  zero,  since  they  should  not  be  used  to 
define  the  slope  of  AQ  versus  relationships.  The  most  extensive  data  sets 
(moderate  e  at  T  =  -10°C  with  Oc  =  0.1  and  10  MPa)  show  a  linear  relationship,  i.e.,  a 
constant  value  of  dAQ/dD^.  Hence,  linear  regression  was  used  to  fit  data  sets  having 

three  or  more  points,  leading  to  the  lines  shown  in  Figure  4.6  and  having  slopes 
given  in  parentheses.  In  general,  faster  strain  rates  and  lower  temperatures  cause  a 

downward  shift  in  the  AQ  versus  Dj.  relationships. 

Figure  4.7  represents  an  attempt  to  quantify  the  influence  of  relative 
density  on  peak  strengths  as  a  function  of  confining  pressure,  strain  rate  and 

temperature.  This  analysis  applied  linear  regression  to  both  AQ  (giving  dAQ/ dDj.) 
and  Qp  data  (giving  dQp/dD^,  but  only  for  tests  with  curve  Types  C  and  D).  For 
simplicity,  both  ratios  will  be  refered  to  as  dQ/ dDj,  even  though  dQp/ dD,  usually 
had  slightly  higher  values  than  dAQ/ dDj  for  unexplained  reasons.  The  results  in 
Figure  4.7  show  two  consistent  trends: 
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•  The  moderate  e ,  -10°C  data  in  Figure  4.7a,  produce  a  linear  increase  in 
dQ/dDj  with  increasing  confinement,  with  a  slope  of  0.67  based  on  the 

collective  AQ  and  Qp  analyses. 

•  At  Oc  =  10  MPa,  the  slow  and  moderate  strain  rate  data  produce  similar 
values  of  dQ/dD, «  11.9  ±  0.6  MPa  independent  of  temperature. 

In  contrast,  the  values  of  dQ/ dD^  at  =  0.1  are  much  more  scattered  than  at  high 
confinement  as  a  function  of  e  and  T.  The  collective  results  give  dQ/ dD^  =  5.6  ±  1.0 
MPa.  For  slow  shearing,  the  mean  data  at  a,.  =  0.1  and  10  MPa  produce  a  slope  of 
0.67,  which  is  the  same  slope  as  that  for  the  moderate  rate  data  at  T  =  -10°C. 

In  summary,  for  test  conditions  leading  to  curve  Types  C  and  D, 

dQ/dDj.  increases  from  5.6  ±  1.0  MPa  at  low  confinement  to  11.9  ±  0.6  MPa  at  high 
corrfinement  for  shearing  at  the  slow  and  moderate  strain  rates  (fast  shearing 
produced  mostly  Type  A  curves  that  preclude  analysis). 

Baker  and  Konrad  (1985)  provide  the  only  prior  data  with  detailed 
measurements  of  peak  strength  versus  relative  der\sity.  Their  data  on  pluviated  16- 
100  Ottawa  sand  tested  in  unconfined  compression  at  e  =  1.7  x  10"^/ sec  and  T  =  -10°C 

give  dQp/dDj »  12.5  MPa  at  Dj  >  65%.  This  value  is  double  that  for  frozen  MFS  at 

low  confinement,  even  though  both  sands  have  similar  frozen  strengths  at  =  90% 
(18  versus  13  MPa  for  MFS).  This  suggests  the  need  for  further  research  on  the 
influence  of  density  on  the  peak  strength  of  frozen  sands. 
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Effect  of  Confining  Pressure 


Figure  4.8  plots  Qp  and  AQ  =  Qp  -  Quy  data  versus  confining  pressure  for 
shearing  at  the  moderate  strain  rate  and  T  =  -10°C.  Linear  regression  was  used  to 
obtain  values  of  Qp  and  AQ  at  Dj.  =  30,  60  and  90%.  As  for  the  analysis  of  relative 

density,  dQ/dOj.  refers  to  ratios  in  terms  of  both  Qp  and  AQ.  Figure  4.8  also  shows 
values  of  the  Mohr-Coulomb  friction  angle  (<}>),  that  were  calculated  using  the 
relationship 

sin  <1)  =  (dQ/da<.)/(dQ/dac  +  2)  . [4.1] 

The  results  in  Figure  4.8  show  the  following  trends: 

•  For  dense  sand  (D^  =  90%),  the  Qp  data  produce  a  linear  increase  in 

strength  with  confinement  (<])  =  14°),  whereas  the  AQ  data  produce  a  slight 
curvature  (<1)  decreases  from  18°  to  12°). 

•  A  decrease  in  relative  density  causes  the  Qp  and  AQ  versus  relationships 
to  become  both  more  curved  and  flatter.  That  is,  the  pressure  sensitivity 
decreases  at  lower  density  and  higher  confinement,  with  <|)  decreasing  to 
about  5°  or  less. 

The  results  in  Figure  4.8  represent  the  only  data  set  having  tests 
sheared  at  Oc  =  0.1,  2,  5  and  10  MPa.  Consequently,  analyses  to  determine  how  the 
pressure  sensitivity  varies  with  strain  rate  and  temperature  had  to  rely  on  values  of 
dQ/d  <Tc  based  only  on  strengths  determined  at  =  0.1  and  10  MPa,  which  can  be 
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somewhat  misleading  at  the  lower  densities  due  to  the  curvature  in  Qp  and  AQ 
versus  D,.  Values  of  Qp  and  AQ  for  =  30,  60  and  90%  were  obtained  from  linear 

regression  and  then  used  to  compute  dQ/ do^  (but  excluded  tests  having  Type  A 
curves).  Figure  4.9  presents  the  results  of  these  analyses,  which  are  based  on  mean 
values  of  dQp/dOc  and  dAQ/da^  and  shows  the  following  trends: 

•  dQ/dOj.  increases  with  increasing  density,  as  expected  from 
Figure  4.8; 

•  at  the  slow  e,  decreasing  temperature  causes  little  or  no  decrease  in 
pressure  sensitivity;  and 

•  at  the  moderate  e,  decreasing  temperature  causes  a  larger  decrease  in 
pressure  sensitivity,  especially  for  dense  sand. 

In  summary,  the  pressure  sensitivity  (i.e.,  the  value  of  dQ/ dOj.)  of 

frozen  MFS  is  affected  mainly  by  relative  density.  Dense  sand  (D,  =  90%)  exhibits  a 
more  or  less  linear  increase  in  Qp  and  AQ  =  Qp  ■  Quy  'vith  increasing  o^,  leading  to 
friction  angles  on  the  order  of  ([)  =  10°  to  15°  that  tend  to  decrease  at  faster  strain  rates 
and  lower  temperatures.  Lower  densities  lead  to  greater  curvature  in  the  strength 
versus  confining  pressure  relationship,  with  mean  friction  angles  on  the  order  of  (j) 

=  5°  or  less.  At  low  densities,  increases  in  strain  rate  and  decreases  in  temperature 
produce  Type  A  curves  (Qp  =  Quy)  that  have  zero  pressure  sensitivity  (dQ/d  Oc  =  0). 

Several  of  the  programs  listed  in  Table  B  studied  the  effect  of  confining 
pressure  on  peak  strength,  but  only  for  medium-dense  to  dense  sands.  Figure  4.10 
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summarizes  data  for  tests  conducted  at  T  =  -10  or  -12°C  and  having  less  than  20 


MPa.  Four  of  the  five  data  sets  have  dQp/dOc  ratios  ranging  from  0.56  to  0.9  (or  <t)  = 

13°  to  18°),  which  are  in  general  agreement  with  results  for  dense  MFS  (Fig.  4.9).  In 
contrast,  the  tests  by  Alkire  and  Andersland  (1973)  run  on  20-30  Ottawa  sand 

produced  a  much  higher  pressure  sensitivity  (<[)  =  28°).  Since  the  reported  ice 

saturation  was  close  to  100%,  this  behavior  appears  to  be  anomalous.  Shibata  et  al 
also  ran  tests  at  T  =  -30°C  on  dense  Toyoura  sand  at  three  strain  rates  similar  to  those 
used  for  the  current  program.  These  data  did  not  show  a  consistent  change  in 
pressure  sensitivity  with  temperature  or  with  strain  rate,  whereas  Figure  4.9  for 
dense  MFS  indicates  a  decrease  in  pressure  sensitivity  at  faster  strain  rates  and  lower 
temperatures.  Consequently,  the  pressure  sensitivity  of  frozen  sands  does  not 
appear  to  follow  the  same  trends. 

Conditions  Producing  Qp  =  Quy 

As  previously  discussed.  Figure  4.6  shows  a  general  tendency  for  a 
downward  shift  in  AQ  =  Qp  -  Qyy  versus  relative  density  with  increasing  strain  rate 
and  decreasing  temperature.  This  trend  can  be  quantified  by  determining  the  values 
of  Dj  that  produce  Type  A  curves  for  which  Qp  =  Quy  (and  for  Type  B  curves  with 
Q  »  Q  ).  Figure  4.11  does  this  for  low-confinement  tests  based  on  mean  results  of 
analyses  using  both  AQ  and  Qp  data  (i.e.,  as  done  for  Fig.  4.7).  Figure  4.11  suggests  a 
linear  increase  in  the  Dj.  producing  Qp  =  Quy  with  decreasing  temperature  and  that 
this  temperature  dependency  becomes  more  pronounced  at  higher  strain  rates. 
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strain  at  Peak  Strength 


Sheets  B14  and  B15  from  Swan  (1994)  show  the  following  trends  in  the 
axial  strain  at  the  peak  strength  (Ep)  for  Type  C  and  D  curves: 

•  at  low  confinement  (Oc  =  0.1  MPa),  Ep  typically  equals  about  5.5  ±  2%,  more 
or  less  independent  of  D,,  e  and  T  (all  Type  C  curves); 

•  increases  in  confinement  cause  a  large  increase  in  Ep  (e.g.,  see  Sheet  B15); 
and 

•  at  high  confinement  (Oc  =  10  MPa),  Ep  typically  ranges  between  15  and  25% 
and  tends  to  decrease  with  higher  densities,  higher  strain  rates  and  lower 
temperatures. 

4.3.3  Effects  of  Strain  Rate  and  Temperature 

Figure  4.12  plots  peak  strength  and  upper  yield  stress  versus  strain  rate 
for  frozen  MFS  at  two  extreme  conditions: 

(a)  low  density  and  low  confinement  (hence,  minimal  effect  of  frictional 
resistance  of  the  sand  skeleton);  and 

(b)  high  density  and  high  confinement  (hence,  maximum  frictional 
resistance). 

The  results  in  Figure  4.12a  show  relatively  little  difference  between  Qp  and  Qyy 
(except  for  slow  shearing  at  -10°C),  whereas  the  results  in  Fig.  4.12b  show  values  of 
Qp  much  larger  than  Q^y.  This  comparison  illustrates  the  extreme  importance  of 

comparing  these  two  stresses,  since  Section  3.4  concluded  that  Q^y  is  controlled  by 
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the  behavior  of  the  ice  matrix  and  therefore,  involves  very  little  sand  skeleton 
frictional  resistance.  It  also  demonstrates  that  inspection  of  Qp  data  alone  can  lead  to 
erroneous  conclusions  regarding  the  physical  mechanisms  controlling  frozen  sand 
behavior  at  large  strains. 

Figure  4.13  plots  AQ  =  Qp  -  Qyy  for  loose  (Dj  =  30%)  and  dense  (D,  =  90%) 

specimens  (interpolated  from  Fig.  4.6)  versus  strain  rate  at  varying  temperatures  at 
low  and  high  confinement.  Inspection  of  the  results  indicates  the  following  trends 
(using  the  slow  strain  rate  of  -10°C  as  the  reference): 

•  For  loose  sand  at  low  confinement,  the  minimal  strength  gain  rapidly 
vanishes  at  higher  strain  rates. 

•  for  loose  sand  at  high  confinement  and  for  dense  sand  at  low 
confinement,  the  moderate  strength  gain  disappears  more  rapidly  with 
faster  shearing  at  the  lower  temperatures. 

•  For  dense  sand  at  high  confinement,  the  large  strength  gain  is  temperature 
independent  at  the  slow  rate;  lower  temperatures  then  cause  a 
progressively  smaller  strength  gain  at  the  higher  strain  rates. 

4.3.4  Summary  and  Conclusions 

Section  4.3  presents  a  detailed  analysis  of  the  amount  of  strengthening 
that  occurs  after  reaching  the  upper  yield  stress  (Quy)-  Quy  is  taken  as  the  reference 
stress,  since  its  value  does  not  vary  with  sand  density  and  confining  pressure  and  is 
thought  to  be  dominated  by  the  same  physical  mechanisms  that  control  the 

behavior  of  polycrystalline  ice.  The  analysis  focuses  on  combinations  of  Dj.,  a^,  e 


64 


and  T  that  produce  curve  Types  C  and  D  (see  Fig.  2.4  and  Table  2.2)  for  which 
AQ  =  Qp  -  Quy  is  positive. 

1)  Relative  Density  (Figs.  4.5, 4.6  and  4.11) 

There  is  a  linear  increase  in  AQ  with  increasing  D,.  The  rate  of 

increase  (dQ/dD^)  is  roughly  independent  of  e  and  T  and  is  highly  dependent  on 
confinement;  e.g.,  dQ/dD^  =  5.6  ±  1.0  MPa  and  11.9  ±  0.6  MPa  at  =  0.1  and  10  MPa, 
respectively.  Analyses  are  made  to  obtain  values  of  D,  producing  Qp  =  Quy  (i.e., 

AQ  =  0  for  Type  A  curves)  for  tests  run  at  low  confinement  as  a  function  of  e  and  T. 
These  values  of  D,  range  from  about  zero  to  60%  for  slow  and  moderate  strain  rates. 

2)  Confining  Pressure  (Figs.  4.8  and  4.9) 

The  pressure  sensitivity  (dQ/ d  Oj.)  is  much  higher  for  dense  than  for 
loose  sand  and  tends  to  decrease  at  very  high  confining  pressures.  Mean  values  of 
dQ/dOc  (from  data  at  Oc  =  0.1  and  10  MPa)  for  slow  and  moderate  shearing  range 
from  about  0.2  to  0.7  (corresponding  to  friction  angles  of  5  to  15°)  and  vary  in  a 
complex  fashion  with  D,,  e  and  T. 

3)  Strain  Rate  and  Temperature  (Figs.  4.12  and  4.13) 

Figure  4.13  shows  the  post-upper  yield  strength  gain,  AQ  =  Qp  -  Quy, 
for  loose  and  dense  sand  at  low  and  high  confinement  as  a  function  of  strain  rate 
and  temperature.  The  magnitude  of  AQ  varies  in  a  complex  fashion  as  a  function  of 
Dj.,  Of.,  e  and  T.  However,  one  can  say  that  and  dominate  at  lower  strain  rates 
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and  higher  temperatures,  whereas  fast  shearing  and  low  temperatures  produce  little 
strengthening  except  for  dense  sand  at  high  confinement. 

Sheet  B13  evaluates  the  strength  gain  normalized  by  the  upper  yield 

stress,  (Qp  -  Quy)/Quy/  arid  plots  this  ratio  versus  temperature  for  loose  (D^  =  35%) 
and  dense  (D,  =  95%)  sand,  as  a  function  of  strain  rate  and  confining  pressure.  One 
observes  that  (Qp  -Quy)/Quy  undergoes 

•  a  significant  increase  with  increasing  density  (note  the  different  scales  for 
loose  and  dense  sand),  as  expected; 

•  a  significant  increase  with  increasing  confining  pressure,  as  expected;  and 

•  a  significant  increase  with  decreasing  strain  rate  and  increasing 
temperature.  In  other  words,  the  relative  degree  of  strengthening  is  much 

larger  at  low  values  of  Qyy  than  at  high  values  of  Qyy. 

4)  Relative  Importance  of  Testing  Variables 

Sheet  B12  illustrates  the  relative  importance  of  the  testing  variables, 
taking  medium-dense  sand  sheared  at  low  confinement,  the  moderate  strain  rate 
and  T  =  -15°C  as  the  reference  condition  (Type  C  curve).  Both  an  order  of  magnitude 
change  in  strain  rate  and  a  5°C  change  in  temperature  have  a  large  effect  on  the  peak 
strength,  due  mainly  to  changes  in  the  upper  yield  stress  (Quy)-  Changes  in  relative 

density  by  AD^  =  ±  30%  have  a  small  effect  (unless  increases  to  extremely  high 
values). 
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4.4  VOLUMETRIC  BEHAVIOR 


4.4.1  Background 

The  current  program  is  one  of  the  few  to  measure  volumetric  strains 
(Ev)  and  only  Shibata  et  al.  (1985)  report  Ey  data  at  varying  confining  pressures,  strain 
rates  and  temperatures  (see  Table  B).  Both  programs  show  essentially  zero 
volumetric  strain  at  low  axial  strains  (up  to  Q^y)  and  then  varying  amoimts  of 
volumetric  expansion  (dilation)  that  continues  to  increase  with  further  straining, 
e.g.,  see  Figures  4.1  through  4.3. 

An  ice-saturated  system  having  a  Poisson's  ratio  less  than  0.5  and 
sheared  by  increasing  the  axial  stress  at  constant  cell  pressure  can  undergo  only 
volumetric  compression.  Consequently,  volumetric  expansion  (dilation)  must 
reflect  fracturing  (cracking)  of  the  ice  matrix  and  loss  of  bonding  between  the  ice 
matrix  and  the  sand  particles.  In  other  words,  the  onset  of  dilation  suggests  a 
decrease  in  the  cohesive  strength  of  the  frozen  sand  system  caused  by  damage  of  the 
ice  matrix. 

The  presentation  uses  two  parameters  to  quantify  the  dilation  behavior 
at  large  strains: 

•  the  maximum  rate  of  dilation  (MRD)  =  the  maximum  slope  of  the 
versus  Eg  curve  =  (dEy/d  ea)inax/ 

•  the  volumetric  expansion  at  20%  axial  strain  =  e^jo- 
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4.4.2  Maximum  Rate  of  Dilation 

The  maximum  rate  of  dilation  (MRD)  is  affected  mainly  by  the  amoimt 
of  confinement,  as  illustrated  in  Figure  4.14a  for  loose  and  dense  specimens  sheared 
at  the  moderate  strain  rate  and  T  =  -10°C.  For  this  most  extensive  data  set,  MRD 

equals  0.5  ±  0.2  at  o,.  =  0.1  MPa  and  drops  to  less  than  0.1  at  Oc  =  10  MPa. 

Figures  4.14b  and  4.14c  show  the  influence  of  relative  density  at  =  0.1 
and  10  MPa,  respectively,  for  all  strain  rates  and  temperatures.  At  low  confinement, 
the  MRD  increases  with  sand  density  more  of  less  independent  of  strain  rate  and 
temperature.  At  high  confinement,  lower  temperatures  tend  to  increase  MRD,  and 
sand  density  plays  a  very  minor  role. 

4.4.3  Volumetric  Expansion  at  20%  Axial  Strains 

The  Eyjo  parameter  follows  the  same  basic  trends  as  MRD,  as  illustrated 

in  Figure  4.15.  At  low  confinement,  increases  with  Dj.  from  about  3  to  8%  and 
tends  to  increase  at  higher  strain  rates  and  lower  temperatures.  At  high 
confinement,  e^jo  =  1-05%  ±  0.85SD  independent  of  D^,  but  with  a  tendency  to 
increase  at  the  fast  strain  rates. 

It  is  interesting  to  note  that  loose  specimens  at  low  confinement 
experience  about  twice  as  much  dilation  at  20%  axial  strain  as  unfrozen  MFS  (Swan 
1994).  These  specimens  often  expanded  to  relative  densities  that  approach  zero 
percent.  Hence,  fracturing  of  the  ice  matrix  increases  the  dilation  of  frozen  sand. 
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4.5  POST-PEAK  STRAIN  SOFTENING 


4.5.1  Objective 

The  combination  of  lubricated  end  platens,  which  enhances  the  ability 
to  obtain  reliable  stress-strain  data  to  very  large  strains,  and  measurement  of 
volumetric  strains  offered  a  unique  opportunity  to  evaluate  the  post-peak  strain 
softening  behavior  of  a  frozen  sand.  The  analyses  looked  at  a  variety  of  parameters 
in  an  attempt  to  quantify  the  degree  of  strength  loss  as  a  function  of  the  four  testing 

variables  (D^,  e  and  T),  with  particular  emphasis  on  the  role  of  volumetric 
expansion  (dilation).  The  stress  parameters  used  in  the  analyses  included  various 
forms  of  the  rate  of  strain  softening,  dQ/ deg,  and  stress  ratios  such  as  Q20/ Qp 

Q2o/Quy/  where  Q20  equals  (Oi  -  03)  at  20%  axial  strain.  Parameters  selected  to 
represent  the  volumetric  behavior  included  the  maximum  rate  of  dilation, 

MRD  =  (d  Ey/d  ea)max/  volumetric  strain  at  20%  axial  strain,  e^jo- 

The  evaluation  focuses  on  tests  conducted  at  Oj.  =  0.1  MPa,  since  low 
confinement  produces  post-peak  strain  softening  of  two  types:  curve  Type  A  with 
Qp  =  Quy  occurring  at  small  strains  (ep=  Euy  <  1%);  and  curve  Type  C  with  Qp  >  Quy 

occurring  at  moderate  strains  (Ep  =  5.5  ±  2%). 

4.5.2  Results  of  Analyses 

Behavior  at  Low  Confinement 
Figure  4.16  plots  the  normalized  rate  of  strain  softening, 

NRSS  =  [(Qp  -  Q2o)/Quy]/(  £20  •  Ep)/  versus  relative  density  from  24  low-confinement 
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tests  having  Type  C  curves.  The  data  produce  a  well-defined  trend,  in  spite  of  large 
changes  in  both  the  strain  rate  and  the  temperature.  The  increase  in  NRSS  with 
relative  density  comes  from  the  fact  that  Qp/Quy  increases  with  Dj.  and  Q2o/Quy 
essentially  constant  (e.g.,  see  normalized  stress-strain  curves  in  Figure  4.1  and 
subsequent  discussion).  However,  the  apparent  independence  of  NRSS  on  strain 

rate  and  temperature  at  the  same  D,  was  rather  surprising.  This  independence 
results  from  two  basic  trends  that  cancel  each  other,  as  demonstrated  below: 

•  The  normalized  peak  strength  ratio,  Qp/Quy/  decreases  with  increasing 
strain  rate  and  decreasing  temperature,  as  illustrated  in  Figure  4.17. 

•  However,  the  normalized  strength  at  20%  axial  strain,  Q20/ Quy/ 
decreases  with  increasing  strain  rate  and  decreasing  temperature,  as 
illustrated  in  Figure  4.18. 

Note  that  Figure  4.18  also  includes  data  from  Type  A  curves  and  that  both  curve 
types  end  up  with  the  same  normalized  stress  at  20%  axial  strain  independent  of 
sand  density. 

Rate  of  Strain  Softening 

Figure  4.19  plots  the  rate  of  strain  softening,  RSS  =  -dQ/de^,  versus  the 
maximum  rate  of  dilation,  MRD  =  (d  Ey  /d  Eg)  max-  The  RSS  is  measured  at  the  axial 

strain  corresponding  to  the  MRD,  which  typically  occurred  at  Eg  =  15  ±  5%.  The 
figure  includes  data  for  all  curve  types  and  for  all  testing  conditions.  The  results 
show  the  following  trends: 
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•  Most  of  the  specimens  having  Type  D  curves  are  still  strain  hardening  at 
MRD;  hence,  they  have  negative  values  of  RSS.  All  these  data  come  from 

tests  with  high  confinement  (a^  =  10  MPa  and  a  few  at  2  and  5  MPa). 

•  The  Type  B  curves  are  mostly  for  loose  specimens  with  high  confinement 
and  tested  at  the  higher  strain  rates  and  lower  temperatures.  As  for  the 
Type  D  curves,  the  RSS  is  typically  less  than  10  MPa. 

•  The  Type  A  and  C  curves  produce  significant  scatter  in  the  data,  but  with  a 
clear  trend  of  an  increasing  rate  of  strain  softening  with  increasing 
maximum  rate  of  dilation  and  no  consistent  difference  between  the  two 
curve  types.  Linear  regression  on  these  data  suggests  somewhat  higher 
values  of  RSS  at  the  lower  temperatures.  Note:  normalization  of  RSS  to 
Qyy  does  not  reduce  the  scatter  by  a  significant  amount. 

As  previously  emphasized  in  Section  4.4.1,  the  volumetric  expansion 
(dilation)  of  frozen  sand  must  reflect  fracturing  of  the  ice  matrix  such  that  the 
specimen  behaves  more  like  a  drained  shear  test  on  a  particulate  material. 

However,  for  drained  triaxial  compression  tests  on  moderate  to  dense  specimens  of 
unfrozen  MFS,  the  peak  strength  occurs  when  the  sand  is  experiencing  its 
maximum  rate  of  dilation  (Swan  1994).  Thus,  frozen  sand  behaves  opposite  to 
unfrozen  sand,  since  the  maximum  rate  of  dilation  for  the  fractured  sand-ice  system 
corresponds  to  a  post-peak  decrease  in  resistance. 
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strength  Loss 


Figure  4.20  plots  the  stress  ratio  Q2o/Qp  versus  the  volumetric  strain  at 

20%  axial  strain,  The  trends  are  similar  to  those  discussed  for  Figure  4.19.  That 
is.  Type  A  and  C  curves  show  larger  strength  losses  (a  lower  stress  ratio)  with  higher 
amoimts  of  dilation,  again  with  significant  scatter  in  the  data. 

Normalized  Rate  of  Strain  Softening 
Figure  4.16  shows  a  very  consistent  relationship  between  the 

normalized  rate  of  strain  softening,  NRSS  =  [(Qp  -  Q2o)/Quy]/(  ^20  "  relative 

density  for  all  tests  at  Oj,  =  0.1  MPa  having  Type  C  curves.  Although  not  plotted,  tests 
having  Type  A  curves  did  not  fit  this  trend.  For  example,  seven  Type  A  curves 
from  tests  having  D,  *  35  to  40%  produced  NRSS  values  ranging  from  0.9  to  2.65. 

However,  correlating  NRSS  with  the  amount  of  dilation  at  20%  axial  strain 
gives  a  very  consistent  relationship  for  both  curve  types.  This  is  demonstrated  in 
Figure  4.21,  which  probably  represents  the  most  significant  finding  from  the  analysis 
of  strain  softening  behavior. 

4.5.3  Summary  and  Conclusions 

The  use  of  lubricated  end  platens  and  measurement  of  volumetric 
strains  have  made  possible  a  unique  evaluation  of  the  post-peak  strain  softening 
behavior  of  frozen  MFS.  For  tests  at  low  confinement  that  produce  curve  Types  A 
and  C  (i.e.,  post-peak  strain  softening),  there  appears  to  be  a  unique  relationship 
between  the  rate  of  strain  softening  (after  normalization  to  the  upper  yield  stress) 
and  the  amount  of  volumetric  expansion  (dilation).  Figure  4.21  presents  this 
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relationship  in  terms  of  the  normalized  rate  of  strain  softening, 

NRSS  =  [(Qp  -  Q2o)/Quyl/(  ^20 "  Cp)/  versus  the  amount  of  dilation  at  20%  axial  strain, 

£^20-  The  collective  data,  which  cover  a  wide  range  of  relative  densities  (Dj  =  35  to 

95%),  have  relatively  little  scatter  about  the  mean  trend  and  do  not  show  a 
consistent  deviation  as  a  function  of  strain  rate  or  temperature.  However,  the 
individual  stress  components  of  the  NRSS  do  vary  significantly  with  the  testing 
variables,  as  summarized  below  (for  curve  Types  A  and  C): 

•  Qp/Quy  increases  linearly  with  relative  density;  these  lines  move 
downward  and  become  flatter  with  increasing  strain  rate  and  decreasing 
temperature.  (See  Fig.  4.17  for  Type  C  curves.) 

•  Q2o/Quy  is  independent  of  relative  density,  but  also  has  a  smaller 

magnitude  with  increasing  strain  rate  and  decreasing  temperature.  (See 
Fig.  4.18  for  both  Type  A  and  C  curves.) 

Analysis  of  test  results  at  all  confining  pressures  (i.e.,  up  to  =  10  MPa) 

produced  the  results  presented  in  Figures  4.19  and  4.21.  For  all  curve  types  (i.e., 
including  Type  B  and  D  with  little  or  no  post-peak  strain  softening,  the  rate  of  strain 

softening  increases  with  the  maximum  rate  of  dilation,  and  the  stress  ratio  Q20/ Qp 
decreases  with  the  amount  of  dilation.  Although  both  plots  have  significant  scatter, 
they  provide  additional  evidence  that  the  degree  of  post-peak  strain  softening  or 
strength  loss  is  strongly  related  to  the  rate  or  amount  of  dilation. 
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4.6  CONCEPTUAL  MODEL  OF  THE  PHYSICAL  MECHANISMS  GOVERNING 
THE  LARGE-STRAIN  BEHAVIOR  OF  FROZEN  SAND 
4.6.1  Introduction 

Section  3.4  showed  that  the  upper  yield  stress  (Quy)  of  frozen  MFS  is 
independent  of  relative  density  (Dj)  and  confining  pressure  (Oc)  and  varies  with 

strain  rate  (e)  and  temperature  (T)  in  a  fashion  similar  to  that  of  polycrystalline  ice. 
The  volumetric  strain  (e^)  at  Q^y  is  also  essentially  equal  to  zero.  Section  4.6.2  gives 
a  brief  overview  of  ice  behavior  and  suggests  possible  mechanisms  to  help  explain 
why  Quy  is  greater  than  the  peak  strength  of  ice. 

In  contrast,  the  post-upper  yield  behavior  of  frozen  sand  is  highly 
dependent  on  sand  density  and  confinement,  in  addition  to  strain  rate  and 
temperature.  The  four  curve  Types,  A  through  D,  are  used  to  illustrate  this 
dependency  on  Dj.,  £  and  T  (Table  2.2  and  Fig.  2.4).  Section  4.6.3  presents  a 
conceptual  model  to  explain  the  strain  hardening/strain  softening  and  volumetric 
behavior  associated  with  these  four  curve  types  in  terms  of 

•  weakening  mechanisms  associated  with  damage  to  the  ice  matrix  due  to 
cracking  that  is  reflected  by  dilation  (expansion)  of  the  frozen  sand;  and 

•  strengthening  mechanisms  associated  with  the  frictional  resistance  of  the 
frozen  sand  system. 

Section  4.6.4  then  examines  the  large-strain  behavior  of  frozen  MFS  in  terms  of 
these  weakening  and  strengthening  mechanisms. 
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4.6.2  Ice  Behavior  and  Upper  Yield  Strength  of  Frozen  Sand 

Figure  4.22  illustrates  typical  stress-strain  curves  for  polycrystalline  ice 
for  the  range  of  testing  conditions  used  for  frozen  MFS.  (The  Mellor  and  Cole  1982 
data  at  T  =  -5°C  should  be  representative  of  ice  at  T  =  -10°C  at  the  slow  and  fast  strain 
rates.)  According  to  Mellor  and  Cole  (1982)  and  Murrell  et  al.  (1989),  the  onset  of 
internal  cracking  (crack  nucleation)  begins  at  the  first  yield  point  (at  axial  strains  on 
the  order  of  0.05  to  0.5%).  The  continued  formation  and  growth  of  cracks  result  in 
the  peak  strength  occurring  at  about  0.5  to  1%,  followed  by  strain  softening  with 
dQ/ dEg  becoming  less  negative  with  increasing  strain.  The  amount  of  post-peak 

strain  softening  increases  at  higher  strain  rates  and  lower  temperatures  and  is 
reduce  by  added  confinement.  However,  in  all  cases  of  interest,  there  is  a  very 
significant  post-peak  loss  in  the  strength  of  polycrystalline  ice  having  a  grain  size  of 
about  1  mm. 

Figure  3.7  shows  that  the  upper  yield  stress  (Quy)  of  frozen  MFS  is 
approximately  double  the  peak  strength  of  ice  in  unconfined  compression  at 
T  =  -10°C.  If  Quy  is  primarily  controlled  by  the  cohesive  strength  of  the  ice  matrix,  as 
hypothesized  in  Section  3.4,  then  the  presence  of  the  sand  must  cause  significant 
strengthening.  The  authors  agree  with  Ting  et  al.  (1983),  who  attribute  this 
enhanced  strength  of  the  intact  ice  matrix  to 

•  a  difference  in  ice  structure.  Although  the  grain  size  of  ice  in  frozen  sand 
is  probably  at  least  an  order  of  magnitude  smaller  than  typically  tested,  no 
data  exist  to  prove  that  a  smaller  grain  size  will  necessarily  produce  a 
higher  strength  of  intact  granular  ice. 
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•  an  increase  in  the  strain  rate  due  to  the  presence  of  the  sand  grains. 

•  changes  in  the  state  of  stress  and  the  added  deformational  restraints  due  to 
the  presence  of  the  sand  grains. 

Since  the  peak  strength  of  normal  granular  ice  increases  with  confinement  at  the 
higher  strain  rates  and  lower  temperatures  (due  to  more  prevalent  cracking 
activity),  one  might  expect  the  upper  yield  stress  of  frozen  sand  to  increase  with 
confining  pressure.  However,  this  was  not  observed  for  frozen  MFS,  although  it  did 
occur  for  the  fast  strain  rate  tests  of  Chamberlain  et  al.  (1972)  as  per  Figure  3.6. 

Perhaps  the  Q^y  of  frozen  MFS  at  high  strain  rates  and  low  temperatures  occurs 
when  the  ice  matrix  is  closer  to  its  first  yield  stress,  which  may  be  less  pressvure 
sensitive  than  the  peak  strength.  In  any  case,  the  authors  have  assumed  that  the  ice 
matrix  in  frozen  sand  follows  the  same  general  trends  as  illustrated  in  Figure  4.22. 

4.6.3  Post-Upper  Yield  Weakening  and  Strengthening  Mechanisms 

After  frozen  sand  reaches  its  yield  stress,  it  is  hypothesized  that  the 
cohesive  strength  of  the  ice  matrix  decreases  because  of  the  continued  formation 
and  growth  of  cracks.  This  cracking  activity  will  decrease  the  compressive  and 
tensile  strength  of  the  ice  and  will  also  decrease  the  bonding  between  the  ice  matrix 
and  the  sand  particles,  thus  permitting  volumetric  expansion  (dilation)  of  the  ice- 
sand  composite  material.  Hence,  continued  cracking  of  the  ice  matrix  represents  a 
weakening  mechanism  that,  by  itself,  will  cause  post-upper  yield  strain  softening  of 
the  frozen  sand  composite  such  as  observed  in  Type  A  curves.  It  is  fmther 
hypothesized  that  increasing  rates  of  dilation  reflect  greater  cracking  activity  and 
more  pronounced  weakening. 
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On  the  other  hand,  continued  straining  of  frozen  sand  beyond  Quy  also 

can  produce  a  strengthening  mechanism  due  to  the  frictional  resistance  of  the  ice- 
sand  composite  that  progressively  behaves  more  like  a  particulate  system  due  to 
increasing  damage  of  the  ice  matrix  and  increasing  interaction  between  the  particles 
of  the  sand  skeleton.  Although  the  ice  matrix  prevents  the  sand  skeleton  from 
exhibiting  the  same  behavior  as  observed  for  unfrozen  sand,  it  still  may  be  valid  to 
use  this  behavior  for  general  guidance.  Based  on  data  summarized  in  Section  5, 

(e.g..  Figs.  5.2  to  5.6),  one  observes  the  following  behavioral  trends: 

1)  For  shear  at  constant  volume  (perhaps  applicable  at  low  to  moderate 
strains  for  frozen  sand,  i.e.,  before  significant  dilation): 

•  loose  sand  at  low  consolidation  stresses  (<Tc  =0.1  MPa)  requires  very  large 

strains  to  reach  a  very  low  peak  strength  (say  Qp  «  1  MPa  at  Ep  =  25%); 

•  dense  sand  at  low  consolidation  stresses  requires  very  large  strains  to  reach 
a  moderate  peak  strength  (say  Qp  =  5  MPa  at  =  25%);  and 

•  dense  sand  at  high  consolidation  stresses  (a’^  =  10  MPa)  requires  relatively 

small  strains  to  reach  a  moderate  peak  strength  (say  Qp  =  5  MPa  at  Eg  =  1  to 

2%),  followed  by  continued  straining  at  constant  Q. 

It  needs  to  be  emphasized  that  initial  application  of  confining  pressure  to  frozen 
sand  does  not  cause  a  corresponding  increase  in  the  consolidation  stress  (effective 

stress)  acting  on  the  soil  skeleton,  since  all  or  most  of  is  probably  carried  by  the 

incompressible  ice  matrix.  It  also  should  be  noted  that  Ladanyi's  (1985)  dilatancy- 
hardening  theory  discussed  in  Section  5  assumes  constant  volume  shear  of  the  sand 
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skeleton  and  intact  properties  for  the  ice  matrix.  For  example,  the  tensile  strength  of 
ice  is  assumed  to  provide  significant  added  confinement  and,  hence,  strength  during 
constant  volume  shear  of  dense  sand  sheared  at  low  confinement. 

2)  For  shear  at  constant  effective  confining  stress,  i.e.,  so-called  drained  shear 
(perhaps  applicable  at  moderate  to  large  strain  for  frozen  sand,  i.e.,  during  significant 
dilation): 

•  both  loose  and  dense  sand  at  low  confinement  (a'c  =  0.1  MPa)  have  a 
negligible  peak  strength  (Qp  <  0.2  MPa),  which  occurs  at  moderate  strains 
(say  Eg  =  5%);  and 

•  with  increasing  confinement,  there  is  a  proportional  increase  in  strength 
and  strain  at  failure  (say  Qp  =  10-20  MPa  at  =  10-20%  for  a'c  =  5-10  MPa  ). 
In  other  words,  drained  shear  at  high  confinement  (a'^.  =  10  MPa)  produces 

very  significant  strengths,  compared  to  ice  and  frozen  sand.  However, 
these  high  strengths  also  require  significant  volumetric  compression 
during  shear. 

In  addition  to  the  strengthening  mechanism  caused  by  the  frictional 
resistance  of  the  sand  skeleton  (as  described  above),  one  can  envision  strengthening 
due  to  interference  between  the  sand  particles  and  the  damaged  ice  matrix,  and 
perhaps  some  frictional  resistance  of  the  fractured  ice  matrix.  In  any  case,  these 
strengthening  mechanisms  are  attributed  to  various  forms  of  frictional  resistance.  It 
is  further  hypothesized  that  these  strengthening  mechanisms  for  frozen  sand  will 
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•  increase  with  increasing  relative  density,  with  this  effect  being  more 
important  for  undrained  than  for  drained  shear;  and 

•  increase  with  increasing  confining  pressure,  with  this  effect  being  much 
more  important  for  drained  than  for  undrained  shear. 

Furthermore,  based  on  the  observed  volumetric  behavior  of  frozen  MFS,  one  might 
expect  the  ice-sand  system  to  behave  in  an  undrained  (constant  volume)  manner  at 
low  to  moderate  strains.  For  this  condition,  the  effective  stresses  causing  frictional 
resistance  are  generated  within  the  ice-sand  system.  At  larger  strains,  as  dUatancy 
becomes  more  pronounced,  the  ice-sand  system  would  be  expected  to  approach  a 
drained  shear  condition  wherein  the  effective  stresses  causing  frictional  resistance 
are  generated  by  the  applied  confining  pressure. 

4.6.4  Post-Upper  Yield  Behavior  of  Frozen  MFS 

Behavior  at  Low  Confinement  (Cc  =  0.1  MPa) 

Type  C  curves  represent  the  most  prevalent  type  of  behavior.  After 
reaching  the  upper  yield  stress,  the  frozen  sand  exhibits  strain  hardening,  since  the 
strengthening  mechanisms  are  stronger  than  the  weakening  mechanisms.  The 
degree  of  strengthening  and  the  resultant  peak  strength  (Qp)  increase  with 

increasing  relative  density  (D^),  as  illustrated  in  Figures  4.1,  4.6  and  4.17,  probably 

similar  to  imdrained  shear.  With  continued  straining,  dilatancy  reduces  the 
internal  frictional  resistance,  and  the  weakening  mechanisms  now  predominate  to 
cause  strain  softening.  Moreover,  the  normalized  rate  of  strain  softening  is 
uniquely  related  to  the  amount  of  dilation,  as  per  Figure  4.21. 
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Faster  strain  strains  and  lower  temperatures  produce  Type  A  curves, 
e.g.,  test  number  36  in  Figure  2.2a  and  test  number  97  in  Figure  4.1b,  especially  at  the 
lower  sand  densities.  The  weakening  mechanisms  dominate,  due  to  both  a  more 
rapid  rate  of  ice  cracking  and  the  relatively  low  post-yield  frictional  resistance.  It  is 
important  to  note  that  the  normalized  rate  of  strain  softening  versus  volumetric 
expansion  relationship  for  Type  A  curves  is  the  same  as  for  Type  B  curves  (Fig.  4.21). 
However,  Type  A  curves  exhibit  a  more  rapid  initial  strength  loss  immediately  after 

yielding  at  Qp  =  Quy  than  occurs  after  the  peak  strength  for  Type  C  curves.  But  both 
curve  types  show  an  almost  linear  decrease  in  strength  and  increase  in  dilation  at 
the  larger  strain  levels  (e.g..  Fig.  4.1). 

Behavior  at  High  Confinement  (Cc  =  10  MPa  ) 

Type  D  curves  represent  the  most  prevalent  type  of  behavior.  After 
reaching  the  upper  yield  stress,  the  frozen  sand  exhibits  continued  strain  hardening 
until  reaching  the  peak  strength  at  large  strains.  The  amount  of  strength  increase 
with  increasing  relative  density  is  much  larger  than  for  tests  at  low  confinement 
(Fig.  4.6),  since  high  confinement  increases  the  frictional  resistance  (more 
strengthening)  and  reduces  the  degree  of  ice  cracking  (less  weakening). 

Faster  strain  rates  and  lower  temperatures  produce  Type  B  curves,  e.g.,  test 
number  129  in  Figure  4.3b,  especially  at  the  lower  sand  densities.  The  weakening 
mechanisms  dominate  immediately  after  yielding  for  the  same  basic  reasons  as  for 
Type  A  curves.  But  thereafter,  the  increased  confinement  produces  sufficient 
frictional  resistance  to  offset  damage  of  the  ice  matrix. 
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Fig.  4.3  Normalized  Stress-Strain  Data  for  Frozen  MFS  Showing  the 
Effect  of  Relative  Density  at  High  Confinement  (cr^=10  MPa): 
(a)  Moderate  Strain  Rate  and  T=-10‘’C:  (b)  Fast  Strain  Rate 
and  T=-15°C  [x=peak  stress  ratio] 
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Fig.  4.4  b  Comparison  of  Peak  Strength  vs.  Strain  Rate  of  Dense  MFS 
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Fig.  4.7  Effect  of  Confining  Pressure  on  Values  of  dQ/ dD^  at 
Varying  Strain  Rates:  (a)  Moderate  Strain  Rate; 

(b)  Slow  and  Fast  Strain  Rates 
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Effect  of  Relative  Density  on  Strength  Gain  with  Confinement 
vs.  Temperature  for  Frozen  MFS  at  Slow  and  Moderate  Strain 
Rates  (Mean  values  from  linear  regression  analyses  on  Qp 
andAQ  data  at  a  =0.1  and  10  MPa) 
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Rate  of  Strain  Softening  vs  Maximum  Rate  of  Dilation  of  Frozen  MFS 
as  Function  of  Curve  Type  for  All  Testing  Conditions 
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Normalized  Rate  of  Strain  Softening  vs.  Volumetric  Strain  at  Varying 
Relative  Density,  Strain  Rate  and  Temperature  for  Frozen  MFS  (Type 
and  C  Curves  at  a  =0.1  MPa) 


Fig.  4.22  Typical  Stress-Strain  Curves  for  Polycrystalline  Ice 


103 


5.  QUANTITATIVE  ANALYSIS  OF  THE  STRENGTH  OF  FROZEN  SAND 

5.1  INTRODUCTION 

Ladanyi  et  al  developed  a  dilatancy-hardening  model  to  predict  the  strength 
of  frozen  dense  sand.  Basically,  it  adds  the  shear  strength  of  the  pore  ice,  which  is 
assumed  equal  to  that  measured  in  triaxial  compression  tests  on  polycrystalline  ice, 
to  the  undrained  shear  strength  of  the  sand  skeleton,  which  is  assumed  equal  to  that 
measured  in  triaxial  compression  tests  on  unfrozen  sand  after  accounting  for  the 
potential  development  of  higher  effective  stresses  provided  by  the  higher  tensile 
strength  of  the  pore  ice  relative  to  unfrozen  pore  water.  Section  5.2  describes  this 
model  and  prior  tests  to  check  its  validity.  Section  5.3  summarizes  the  triaxial 
compression  behavior  of  unfrozen  Manchester  Fine  Sand  (MFS)  and  presents  the 
unfrozen  sand  parameters  that  are  required  for  application  of  the  dilatancy- 
hardening  model.  Section  5.4  presents  results  from  consolidate-freeze  triaxial 
compression  tests  on  dense  frozen  MFS,  in  order  to  illustrate  how  the  magnitude  of 
the  effective  stress  acting  on  the  sand  skeleton  prior  to  freezing  affects  frozen  sand 
behavior.  These  data  also  are  relevant  to  application  of  Ladanyi's  model.  Finally, 
Section  5.5  compares  strengths  predicted  by  the  dilatancy-hardening  model  to  those 
measured  from  both  conventional  and  consolidate-freeze  triaxial  compression  tests 
on  frozen  MFS. 

5.2  DILATANCY-HARDENING  MODEL 
5.2.1  Assumptions  and  Development 

Ladanyi  (1985)  developed  the  dilatancy-hardening  model  to  predict  the 
strength  of  frozen  dense  sand,  and  Ladanyi  and  Morel  (1990)  presented  experimental 
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results  to  check  the  validity  of  the  model.  The  basic  concept  assumes  that  the 
effective  stress  path  followed  by  the  soil  skeleton  in  a  frozen  sand  will  be  the  same  as 
that  for  unfrozen  sand  provided  that  1)  both  systems  are  subjected  to  the  same 
strain  path  (axial  loading  at  constant  confining  pressure  and  constant  volume);  and 

2)  the  sand  skeleton  in  both  systems  starts  from  the  same  state  (density  and  effective 
stress).  Unfrozen  dense  sand  sheared  undrained  at  low  confinement  tends  to  dilate 
(expand  in  volume).  This  causes  development  of  negative  pore  pressures  that 
increase  the  effective  stresses  acting  on  the  sand  which  results  in  increased  frictional 
resistance.  The  dilatancy-hardening  model  makes  the  following  assumptions 
(Ladanyi  and  Morel  1990): 

1)  All  of  the  pore  water  in  the  frozen  sand  is  considered  to  be  frozen. 

2)  The  behavior  of  the  sand  skeleton  is  a  function  of  the  state  parameter  'P  as 
introduced  by  Been  and  Jefferies  (1985)  and  described  in  Section  5.3. 

3)  No  consolidation  occurs  prior  to  shear;  i.e.,  the  tests  are  considered  as 
"unconsolidated."  The  shear  starts  from  a  known  density  and  confining 
pressure. 

4)  So  long  as  the  pore  ice  is  continuous  and  unbroken  during  shear,  the  test 
will  behave  as  "undrained."  When  the  pore  ice  fails,  i.e.,  breaks  up,  the 
test  behaves  as  "drained."  This  breakup  of  the  ice  has  been  associated  with 
axial  strains  of  1  to  2%  in  the  literature. 

5)  As  long  as  the  test  behaves  undrained  (i.e.,  before  the  pore  ice  becomes 
broken),  the  dilatancy-hardening  principles  established  for  unfrozen  sand 
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are  also  applicable  to  the  frozen  sand  skeleton  provided  that  proper 
account  is  made  for  the  tensile  strength  of  the  pore  ice  matrix. 

Before  further  description  of  the  dilatancy-hardening  model,  the 

definition  of  the  critical  confining  pressure,  is  needed.  Seed  and  Lee  (1967)  and 
others  have  observed  that  if  a  dense  sand  is  sheared  in  an  undrained  mode  with  no 
back  pressure  in  a  state  where  it  exhibits  a  tendency  to  dilate,  i.e.,  it  has  a  large  -'F 

parameter,  the  pore  water  pressure  decreases,  and  the  effective  stresses  increase  on 
the  sand  skeleton  until  one  of  two  conditions  occurs.  Either  the  pore  water  pressure 
reaches  the  tensile  strength  of  water,  i.e.,  it  decreases  sufficiently  to  cause  cavitation, 
or  the  effective  stresses  increase  until  the  "state  point"  reaches  the  Critical  Void 
Ratio  Line  (CVRL)  or  alternatively,  the  Steady  State  Line  (SSL).  The  transition 
confining  pressure  between  these  two  conditions  depends  on  the  "critical  confining 
pressure,"  Seed  and  Lee  (1967)  define  the  "critical  confining  pressure"  as  that 

confining  pressure  for  a  given  preshear  void  ratio  which  results  in  no  net  volume 
change  at  the  peak  strength  in  a  drained  triaxial  compression  test  on  unfrozen  sand 
(or  no  net  change  in  pore  pressure  at  failure  for  an  undrained  test).  As  used  by 
Ladanyi  et  al,  a'„  is  the  initial  confining  pressure  that  will  result  in  no  further 

change  in  pore  pressure  at  failure  for  undrained  shear  of  unfrozen  sand  at  a  given 
initial  density  (or  void  ratio).  Note:  Ladanyi  (1989)  acknowledges  that  the  Seed  and 
Lee  definition  "might  have  been  better"  and  that  definition  has  been  adopted  by 
MIT. 
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Although  Ladanyi  et  al.  presented  stresses  in  the  form  Q  =  (oj  -  03)  and 
®oct  =  ^3)/  application  of  the  model  will  be  illustrated  using  the  MIT 

format:  q  =  0.5  (a^  -  03),  p  =  0.5(ai  +  03)  and  p'  =  0.5(a'i  +  a'3).  This  is  done  in  Figure 
5.1,  which  has  been  drawn  to  scale,  based  on  specific  values  for  the  critical  confining 
pressure  (a'c)  and  friction  angle  (<!>')  of  the  sand  skeleton  and  the  shear  strength  (qj) 

and  tensile  strength  (Tj)  of  the  ice  matrix. 

The  undrained  shear  strength  of  the  frozen  sand  skeleton  is  equal  to 
qfs  =  P’f  «'  =  P'f  sin  <!)’/  where  p'f  is  the  value  of  p’  at  failure  and  (])’  is  the  effective 

stress  friction  angle  (assumed  equal  to  that  for  unfrozen  sand).  The  total  stress 
(failure)  envelope  for  the  frozen  sand  skeleton,  TSE  (fs),  has  the  same  initial  slope  as 
the  effective  stress  envelope,  ESE  (us  and  fs),  but  translated  to  the  left  by  an  amount 
equal  to  the  tensile  strength  of  the  pore  ice,  Tj.  For  this  initial  region,  the  sand 
strength  is  controlled  by  the  magnitude  of  the  negative  "pore  pressure"  that  can 
develop  before  cavitation  of  the  pore  fluid. 

The  transition  point  between  the  cavitating  and  non-cavitating 
regimes  occurs  at  a  confining  pressure  (total  stress)  of 

Cc  =  o'er  -  T  . [5-1] 

where  T  is  the  tensile  strength  of  the  pore  fluid.  For  pore  water,  T  =  T^^,  which  is 
about  0.1  MPa;  for  pore  ice,  T  =  Tj,  where  Tj  =  2  MPa  for  this  example.  For  a 
confining  pressure  greater  than  o'er  -  T,  i.e.,  for  an  initial  p  (=  a^)  greater  than  pj  in 
Figure  5.1,  the  total  stress  envelope,  TSE(fs),  assuming  an  undrained  test  condition. 
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becomes  horizontal  with  a  value  of  q  =  qfg.  The  predicted  strength  of  the  soil 
skeleton  in  this  region  is  therefore  solely  controlled  by  the  magnitude  of  the  critical 
consolidation  stress  and  is  calculated  as 

qf,  =  0.5'a’„-(Rf-l)  . [5.2] 

where  Rf  is  equal  to  (a'i/a'3)f  =  (1  +  sin  4)')/(l  -  sin  <()’).  For  a  confining  pressure  less 
than  -  Tj,  the  strength  of  the  sand  is  governed  by  the  applied  confining  pressure 
(03  =  Oc)  and  the  tensile  strength  of  the  pore  ice,  so  that 

qfe  =  0.5(a3  +  Ti)(Rf-l)  . [5.3] 

Note;  Although  the  predicted  value  of  qfs  in  the  cavitating  regime  is 
controlled  by  Tj  and  the  total  confining  stress  (03  =  0^.),  the  undrained  shear  data 
presented  in  Section  5.3  will  show  that  the  amount  of  strain  required  to  reach  q^s  is 
very  dependent  upon  the  magnitude  of  the  preshear  effective  stress,  a  '3  =  o  'c. 

Since  pore  ice  can  support  both  ter\sile  and  shear  stresses,  Ladanyi  and 
Morel  (1990)  propose  that  the  shear  strength  of  the  ice  should  be  added  directly  to 
the  shear  strength  of  the  sand  skeleton,  which  also  has  been  increased  by  the  tensile 
strength  of  the  pore  ice.  Using  their  methodology,  the  strength  of  the  frozen  sand  is 
represented  by  a  total  stress  line,  TSL  (FS  =  fs  +  ice),  where  the  additional  ice  strength 
is  added  to  the  sand  skeleton  strength  directly  above  the  effective  stress  envelope 
(ESE).  For  example,  for  the  non-cavitating  case  starting  at  pj,  the  sand  skeleton  ends 
up  at  point  B  on  the  ESE.  The  shear  strength  of  the  ice,  qj,  is  then  added  to  obtain 
point  C,  which  defines  the  break  in  the  TSL  (FS)  envelope.  In  contrast,  the  MIT 
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methodology  would  use  the  TSE  (FS)  which  represents  the  actual  end  points  of 
predicted  total  stress  paths  (i.e.,  points  q^,  q2  and  qa  in  Figure  5.1). 

The  total  shear  strength  of  the  frozen  system  can  now  be  described  for 
both  the  cavitating  and  non-cavitating  case.  In  the  non-cavitating  region,  the 
predicted  shear  strength  of  the  frozen  sand  is  given  by 

qFs  =  0.5a'„(Rf-l)  +  qi  . [5.4] 

where  a'cr  is  the  critical  confining  pressure  for  the  initial  void  ratio  (density)  of  the 
sand  skeleton.  This  equation  predicts  that  the  strength  of  the  frozen  sand  is 
independent  of  the  actual  confining  pressure  for  =  03  >  p^. 

For  the  cavitating  case,  i.e.,  03  <  o’er  -  T,,  the  shear  strength  of  the 
frozen  sand  (qps)  is  predicted  to  be 

qps  =  qfs  +  qi  =  0-5  ■  (<^3  +  Tj)  •  (Rf  - 1)  +  qj . [5.5] 

where  03  is  the  confining  (cell)  pressure,  Tj  is  the  tensile  strength  of  ice,  and  qj  is  the 
shear  strength  of  the  pore  ice.  Hence,  the  shear  strength  in  this  region  increases 
linearly  with  increasing  confining  pressure  with  the  same  slope  as  the  ESE. 

5.2.2  Prior  Experimental  Validation 

Ladanyi  and  Morel  (1990)  present  the  results  of  drained  triaxial 
compression  tests  on  unfrozen  sand  and  conventional  triaxial  compression  tests  on 
frozen  specimens  of  20-30  Ottawa  sand  to  check  the  dilatancy-hardening  model.  The 

tests  were  performed  on  vibrated,  dense,  saturated  specimens  (D,  =  90  ±  6%)  at  a 
mean  strain  rate  at  failure  of  2.7  x  10'^/ sec.  Confining  pressures  (or  effective 
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consolidation  stresses)  prior  to  shear  ranged  from  zero  to  0.3  MPa.  The  temperature 
for  frozen  tests  was  -5°C.  For  the  frozen  tests,  a  well-pronounced  peak  occurred  at  or 
a  little  above  2%  axial  strain  (based  on  measurement  of  axial  deformation). 

Results  from  23  frozen  tests  showed  an  extremely  high  pressure 

sensitivity  of  dQp/dOc  =  2.75  (<t)  =  35°),  compared  to  data  summarized  in  Figures  4.8, 
4.9,  and  4.10;  an  extremely  high  rate  of  dilation  at  failure  of  de^/ dEg  =  3.4  ±  0.15, 
compared  to  the  MFS  data  in  Figure  4.14;  and  a  mean  peak  strength  of  Qp  «  8.8  MPa. 
Two  tests  on  ice  with  =  0.2  MPa  and  e  =  2.8  x  lO'^/sec  gave  Qi  =  4.9  MPa  at  2%  axial 
strain.  Using  ({>'  =  32°  (from  drained  triaxial  compression  tests  on  unfrozen  sand), 

=  0.15  MPa  and  Tj  =  1.75  MPa  in  Eq.  5.3  gives  =  2  qfg  =  4.0  MPa.  Hence,  the  model 
predicts  Qps  =  4.0  +  4.9  =  8.9  MPa,  which  is  in  excellent  agreement  with  the  measured 
peak  strength. 

5.3  TRIAXIAL  COMPRESSION  BEHAVIOR  OF  UNFROZEN  MFS 
5.3.1  Scope  of  Test  Program 

Sections  2.3  and  2.4  summarize  the  experimental  procedures  used  by 
Swan  (1994)  to  conduct  a  total  of  49  useful  triaxial  compression  tests  on  MFS 
specimens  prepared  by  multiple  sieve  pluviation  (i.e.,  same  deposition  techmque  as 
used  for  the  tests  on  frozen  MFS).  Note  that  this  is  the  first  program  to  evaluate  the 
behavior  of  unfrozen  sand  for  direct  comparison  with  the  behavior  of  frozen  soil. 

After  careful  backpressure  saturation  at  a  consolidation  stress  of  =  0.1  MPa,  the 
following  three  types  of  triaxial  compression  tests  were  run: 
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1)  isotropic  consolidation  to  varying  stresses  followed  by  undrained  axial 
loading  at  constant  cell  pressure,  denoted  as  CIUC  tests; 

2)  anisotropic  consolidation  (K^  =  a'hc/cj'vc  =  to  varying  stresses 

followed  by  undrained  axial  loading  at  constant  cell  pressure,  denoted  as 
CAUC  tests;  and 

3)  isotropic  consolidation  to  varying  stresses  followed  by  drained  shear  at 
coi\stant  p'  =  0.5(a'i  +  a'3),  denoted  as  CIDC  tests.  These  tests  were 

computer  controlled  to  maintain  A<y'3  =  -Acr'i  during  shear. 

The  Series  B  test  program  included  24  CIUC  tests  with  ranging  from 

44  to  104%  and  a'c  =  0.1,  2,  5  and  10  MPa  (i.e.,  same  confining  pressures  as  used  for 
frozen  MFS  sheared  at  the  moderate  strain  rate  and  T  =  -10°C).  The  Series  C  test 
program  included  19  CIUC  tests,  3  CAUC  tests  and  13  CIDC  tests  with  D,  varying 

from  30  to  115%  and  o'c  varying  from  0.1  to  12.5  MPa.  Sheet  C6  in  Appendix  C 
summarizes  typical  compression  curves  for  isotropic  consolidation  up  to  a'c  =  10 
MPa  .  Note  that  consolidation  to  this  high  stress  causes  secondary  compression 
(drained  creep)  and  often  produced  relative  densities  greater  than  100%  (the 
minimum  void  ratio  from  vibration  at  low  confinement  is  0.580). 

5.3.2  Stress-Strain-Effective  Stress  Behavior 

Influence  of  Consolidation  Stress 
Figures  5.2  and  5.3  illustrate  the  very  significant  influence  of 
consolidation  stress  on  the  undrained  shear  behavior  of  dense  MFS  (preshear 
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Dr  =  95%),  as  measured  in  one  CAUC  and  five  CIUC  tests.  For  the  CIUC  tests,  one 

observes  that  increases  in  cr’c  cause  the  following: 

•  a  large  increase  in  the  pseudo  upper  yield  stress  (first  obvious  change  in 
the  slope  of  the  Q  versus  curve  at  roughly  1%  axial  strain).  Sheet  C7 

shows  a  linear  increase  in  Quy  with  a'c- 

•  a  large  decrease  in  the  strain  required  to  reach  (or  nearly  reach)  the  peak 
undrained  strength.  However,  the  magnitude  of  Qp  remains  essentially 
constant. 

•  a  change  from  highly  dilatant  behavior  at  low  stresses  (i.e.,  development 
of  negative  excess  pore  pressures  and  an  effective  stress  path  that  climbs 
up  the  effective  stress  failure  envelope  as  in  test  C-34)  to  highly  contractive 
behavior  at  high  stresses  (i.e.,  development  of  positive  excess  pore 
pressure  and  an  effective  stress  path  that  curves  to  the  left  before  reaching 
the  effective  stress  envelope  as  for  tests  C-22  and  C-25).  Note  that  the  CIUC 
tests  used  a  very  high  back  pressure  in  order  to  prevent  cavitation  in 
highly  dilatant  specimens,  such  as  test  C-34  which  reached  Ug  =  -2  MPa. 

Figure  5.4  shows  that  consolidation  stress  also  has  a  very  significant 
influence  on  the  drained  shear  behavior  of  dense  MFS.  As  for  undrained  shear, 
increases  in  o'c  cause  a  proportional  increase  in  the  initial  "stiffness"  and  a  change 
from  dilatant  (volumetric  expansion)  to  contractive  (volumetric  compression) 
behavior.  However,  increases  in  o'c  also  cause  a  proportional  increase  in  the  peak 
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strength  and  the  strain  required  to  reach  Qp,  whereas  undrained  shear  produces  the 
same  strength  at  large  strains.  Note:  although  the  CIDC  tests  in  Figure  5.4  were  run 
at  constant  p',  drained  shear  with  increasing  axial  stress  and  constant  confining 

pressure  (i.e.,  a'3  =  a'c)  will  show  similar  basic  trends. 

Influence  of  Relative  Density 

Figures  5.5  and  5.6  illustrate  the  influence  of  relative  density  on  the 
undrained  shear  behavior  of  MFS  as  measured  in  CIUC  tests  having  low  and 
moderate  consolidation  stresses,  a'c  =  0.1  and  2  MPa,  respectively.  Sheets  C8  and  C9 
show  the  corresponding  effective  stress  paths.  One  observes  that  increases  in  D, 

have  a  minor  effect  on  the  initial  stiffness  (i.e.,  at  £3  less  than  1%),  but  produce  much 
higher  strengths  at  large  strains.  However,  the  reverse  occurs  for  drained  shear  as 
illustrated  in  Sheet  CIO,  i.e.,  increasing  causes  an  increase  in  initial  stiffness  but 

has  much  less  effect  at  large  strains. 

Implications  for  Dilatancy-Hardening  Model 

Section  3.4  showed  that  changes  in  sand  density  (D^)  and  confining 

pressure  (0^)  did  not  affect  the  upper  yield  stress  (Quy)  of  frozen  MFS,  which  appears 
to  represent  the  onset  of  cracking  of  the  ice  matrix  and  which  occurs  at  less  than  1% 
axial  strain.  Section  4  showed  that  sand  density  and  confining  pressure  have  a  large 
effect  on  the  post-upper  yield  behavior.  If  the  strain-strain  behavior  of  the  sand 
skeleton  is  similar  to  that  measured  for  unfrozen  MFS,  one  would  expect  increases 
in  sand  density  and  confinement  to  exhibit  the  following  pattern: 


113 


Relative  Importance  of  Variable  on 
Frictional  Resistance  of  Sand  Skeleton 
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Effective  Stress  Failure  Envelope 

Figure  5.7  plots  values  of  q  and  p'  at  the  maximum  obliquity  condition 
(i.e.,  when  R  =  a’i/a'3  reaches  it  maximum  value)  from  the  CIUC  and  CAUC  tests 
run  on  unfrozen  MFS.  This  condition  occurs  at  about  1%  axial  strain  at  low 
confinement  (o'^  =  0.1  MPa)  and  increases  to  Eg  =  10%  at  a'c  =  10  MPa.  Linear 

regression  on  these  data  gives  an  effective  stress  friction  angle  of  <|)'  =  34.9°  that  is 

used  in  Section  5.5  for  application  of  Ladanyi's  dilatancy-hardening  model. 

5.3.3  Steady  State  Analysis  and  Critical  Confining  Pressure 

Application  of  Ladanyi's  model  requires  determination  of  the  critical 

confining  pressure  (cr'cr)  in  order  to  define  the  transition  point  between  cavitating 
and  noncavitating  regimes  for  shearing  of  frozen  sand  (e.g.,  Eq.  5.1).  As  stated  in 
Section  5.2.1,  MIT  has  defined  as  that  consolidation  stress  for  a  given  preshear 
density  (void  ratio)  that  produces  zero  excess  pore  pressure  at  failure  during 
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undrained  triaxial  compression.  This  condition  is  equivalent  to  setting  Skempton’s 
A  parameter  equal  to  zero,  i.e.,  Af  =  (Au  -  A  03)/ (Aa^  -  Ao^)  =  0,  where  Au  is  the 
measured  change  in  pore  pressure  at  failure.  Note;  for  the  MIT  CIUC  tests, 

Aa3  =  0  and  Aoj  «  AQ. 

Section  5.5  applies  Ladanyi's  model  to  frozen  MFS  for  a  wide  range  of 
densities,  which  in  turn  requires  determination  of  the  critical  confining  pressure  as 
a  function  of  void  ratio.  Since  Af  varies  with  both  density  and  confinement,  the 

analysis  used  the  state  parameter  v  proposed  by  Been  and  Jefferies  (1985),  in  order  to 
account  for  the  influence  of  these  two  variables.  As  illustrated  in  Figure  5.8,  \|r  is 

defined  as  the  initial  (preshear)  void  ratio  of  a  test  specimen  minus  the  void  ratio 
on  the  steady  state  line  (SSL)  at  the  same  effective  mean  stress.  Been  and  Jefferies 
(1985)  selected  the  first  invariant  of  the  stress  tensor  to  define  the  effective  stress, 

Ii  =  ^'oct  =  +  ^'2  ‘^'3)- 

The  steady  state  line  (SSL)  is  defined  as  the  shearing  condition 
representing  continued  straining  at  constant  void  ratio,  pore  pressure  (for 
undrained  shear  tests)  and  shear  shear  and,  hence,  at  a  constant  value  of  Igg-  Sheets 
Cl  and  C2  tabulate  test  parameters  for  CIUC  and  CAUC  tests  at  or  near  a  steady  state 
condition;  Sheet  Cll  plots  the  values  of  e^s  and  I^s  used  to  define  the  steady  state 
line  for  unfrozen  MFS.  This  steady  state  line  is  also  shown  by  the  dashed  line  in 
Figure  5.8.  Two  points  should  be  emphasized.  First,  the  steady  state  condition 

occurs  at  large  strains  that  often  preclude  continued  strairung  at  exactly  constant  Igg 


115 


due  to  limitations  of  the  triaxial  test;  hence,  selection  of  the  steady  state  condition 
involved  judgment  for  most  of  the  tests.  Second,  the  tests  used  to  select  the  steady 

state  points  plotted  in  Sheet  Cll  had  state  parameters  varying  from  y  =  -0.45 

(preshear  state  far  below  the  SSL  and  hence  highly  dilatant  behavior)  to  y  =  +0.09 

(preshear  state  above  the  SSL  and  hence  contractive  behavior).  Hence,  in  spite  of 
potential  errors  in  selecting  steady  state  points,  the  selected  SSL  should  be  reasonably 
accurate.  For  example.  Sheet  C12  shows  a  very  consistent  relationship  between  the 

peak  undrained  shear  strength  ratio  and  y  for  CIUC  tests. 

Sheets  C3  and  C4  tabulate  various  shear  parameters  (i.e.,  (j)’  at 
maximum  obliquity,  Af  and  the  peak  undrained  strength  ratio)  as  a  function  of  the 
state  parameter  y  for  22  CIUC  tests,  and  Figure  5.9  plots  Af  versus  y.  This  figure 
shows  a  well-defined  relationship  that  gives  Af  =  0  at  y  =  -0.061.  Consequently,  the 
critical  confining  pressure  line  for  unfrozen  MFS  is  located  below  the  SSL  by  a 
difference  in  void  ratio  equal  to  Ae  =  -0.061,  as  shown  by  the  solid  line  in  Figure  5.8. 

This  relationship  predicts  a'^r  =  0.27  and  2.52  MPa  at  =  20  and  95%,  respectively. 

5.4  RESULTS  FROM  CONSOLIDATE-FREEZE  TESTS  ON  MFS 
5.4.1  Objectives  and  Scope 

Figure  5.2  showed  that  increases  in  the  preshear  consolidation  stress 
(o'c)  cause  an  enormous  increase  in  the  initial  stiffness  (initial  slope  of  the  Q  versus 
Ea  curve)  of  consolidated-undrained  triaxial  compression  tests  run  on  unfrozen 
dense  MFS.  Consequently,  one  could  expect  that  the  magnitude  of  the  effective 
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stress  acting  on  the  sand  skeleton  prior  to  (and  during)  freezing  could  affect  the 
strength-deformation  behavior  of  frozen  sand. 

A  review  of  the  literature  in  Section  3.6  of  Andersen  et  al.  (1992)  and  in 
Section  2.3.4  of  Swan  (1994)  indicates  that  no  tests  have  been  run  to  compare  the 
behavior  of  frozen  sand  measured  in  conventional  triaxial  compression  tests  with 
that  measured  in  tests  wherein  the  soil  was  frozen  after  application  of  confinement. 
However,  the  results  in  Figure  5.10  suggest  that  the  magnitude  of  the  pre-freezing 
effective  stress  can  be  an  important  variable.  The  figure  compares  stress-strain 
behavior  from  two  conventional  tests  run  on  frozen  dense  MFS  in  which  one 
specimen  (No.  55)  was  prepared  by  multiple-sieve  pluviation,  and  the  second 
specimen  (No.  26)  was  prepared  by  wet  tamping.  The  second  specimen  exhibits  a 
much  higher  rate  of  post-upper  yield  strain  hardening,  a  smaller  strain  at  the  peak 
strength,  a  higher  rate  of  post-peak  strain  softening,  and  a  more  rapid  development 
of  dilation.  (Note:  similar  comparisons  at  higher  confining  pressures  show  similar 
trends  but  little  change  in  the  peak  strength.)  One  possible  explanation  for  this 
behavior  is  that  the  high  compaction  effort  used  to  prepare  the  wet  tamped 
specimen  locked  in  significant  effective  stresses  that  were  preserved  during  freezing. 

Section  2.4.3  describes  the  experimental  procedures  used  for  the 
program  of  ten  consolidate-freeze  tests  run  on  dense  MFS  at  the  moderate  strain  rate 

(£  =  3.5  X  lO'^/sec)  and  T  =  -10°C.  This  preliminary  program  had  three  objectives: 

1)  to  develop  an  experimental  technique  that  produces  homogeneous  frozen 
specimens  (as  accomplished  for  the  conventional  tests)  after  consolidation 

to  stresses  ranging  from  o'c  =  0.1  to  10  MPa; 
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2)  to  produce  data  for  comparison  with  the  results  from  conventional  tests  in 
order  to  obtain  a  better  understanding  of  the  role  of  sand  skeleton  frictional 
resistance  on  the  strength-deformation  behavior  of  frozen  sand;  and 

3)  to  determine  if  the  magnitude  of  the  pre-freezing  effective  stress  is  an 
important  variable  when  applying  Ladanyi's  dilatancy-hardening  model 
to  predict  the  strength  of  frozen  sands. 

5.4.2  Experimental  Results  and  Comparison  with  Conventional  Frozen  Tests 

Consolidation  and  Freezing 

Sheets  D1  and  D4  summarize  data  measured  during  the  consolidation 
phase  of  the  consolidate-freeze  (CF)  tests,  which  resulted  in  pre-freezing  relative 
densities  of  D,  =  93.5  ±  6.3%  SD.  Note  that  volumes  changes  in  three  of  the  tests  (CF 
06,  07  and  09)  were  estimated  due  to  leakage  of  cell  fluid  into  the  internal  drainage 
line  that  provided  drainage  from  the  top  cap  (see  sheet  D5). 

As  described  in  Section  2.4.3,  freezing  was  accomplished  by  circulating 
ethylene  glycol  through  a  redesigned  base  pedestal  (Sheet  D5),  while  monitoring  the 
temperature  in  the  cell  fluid  at  three  locations  next  to  the  specimen  and  the  volume 
of  water  flow  from  the  internal  drainage  line  connected  to  the  top  cap.  Sheet  D6 
plots  results  during  freezing  of  Test  CF02,  and  Sheet  D2  summarizes  data  for  all  ten 
tests.  Sheets  D2  and  D7  show  that  the  measured  outflow  of  water  during  freezing 
was  always  less  than  the  calculated  change  in  volume  due  to  expansion  of  the  pore 
water  as  it  turns  into  ice  (based  on  a  9%  volume  increase),  resulting  in  efficiencies 
(ratio  of  measured  to  calculated  volumes)  ranging  from  87%  to  only  20%.  This 
discrepancy  is  problematic,  since  an  efficiency  of  50%  corresponds  to  an  overall 
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volumetric  expansion  of  the  specimen  of  about  1.7%.  Although  the  assumption  of  a 
9%  pore  water  expansion  during  freezing  may  be  too  high,  this  aspect  of  the  freezing 
technique  requires  further  study. 

Small  Strain  Results 

Sheet  D3  tabulates  values  of  Young's  modulus  (E)  and  yield  offset  stress 
(Qyo)  from  six  of  the  ten  tests;  Sheet  D8  compares  these  values  with  data  from 
conventional  tests.  Although  the  mean  E  from  the  CF  tests  is  about  10%  lower,  the 
overall  agreement  is  considered  acceptable,  since  Section  3  concluded  that  E  and  Qyo 

should  not  be  affected  by  the  frictional  resistance  of  the  sand  skeleton. 

The  upper  yield  region  represents  one  of  the  most  unexpected  aspects 
of  the  consolidate-freeze  tests,  because  only  two  of  the  ten  tests  exhibited  a  distinct 
yield  point,  i.e.,  a  plateau  in  the  slope  of  the  Q  versus  curves.  Sheet  D3  shows 

values  of  Quy  of  8.6  and  7.7  MPa  for  CFOl  and  03  with  =  0.27  and  0.17  MPa, 

respectively,  which  agrees  with  the  Q^y  =  8.3  ±  0.6  MPa  for  all  conventional  tests 

sheared  at  the  same  strain  rate  and  temperatures.  The  lack  of  a  distinguishable 
upper  yield  point  in  the  CF  tests  consolidated  to  higher  stresses  is  illustrated  in 

Figure  5.11,  which  compares  Q-e^  curves  for  the  two  types  of  tests  at  varying 

confining  pressures.  The  conventional  (numbered  FRS)  tests  all  show  either  a 
slightly  increasing  or  a  constant  slope  after  yielding,  whereas  the  CF  tests  with 

a'c  >  2  MPa  show  a  continuous  decrease  in  the  slope  of  the  stress-strain  curve.  The 

higher  stress  CF  tests  also  exhibit  a  much  stiffer  response  beyond  a  few  tenths 
percent  axial  strain.  This  suggests  that  the  frictional  resistance  of  the  sand  skeleton 
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becomes  significant  within  the  small  strain  region  (i.e.,  before  yielding  occurs  for 
conventional  tests)  for  specimens  consolidated  to  very  high  stresses  prior  to 
freezing.  The  shapes  of  the  stress-strain  curves  from  CIUC  tests  on  unfrozen  MFS  in 
Figure  5.2  support  this  hypothesis. 

Large  Strain  Results 

Sheet  D9  plots  the  Q  -  -  Eg  response  of  all  ten  consolidate-freeze  tests; 

Figure  5.12  compares  the  "best  quality"  CF  test  at  each  confining  stress  with 
corresponding  data  from  conventional  tests.  Comparison  of  the  two  data  sets,  along 
with  that  for  unfrozen  sand,  leads  to  the  following  observations: 

•  Increasing  confinement  in  the  conventional  tests  causes  a  change  from 
curve  Type  C  to  curve  Type  D,  whereas  all  CF  tests  have  Type  C  curves.  CF 

tests  with  o'c  ^  5  MPa  also  have  similar  rates  of  post-peak  strain  softening, 

whereas  a'^.  =  10  MPa  seems  to  produce  somewhat  less  strain  softening 
(also  see  Sheet  D9).  It  is  interesting  to  note  that  the  rate  of  strain 
hardening  (dQ/dEg)  during  undrained  shear  of  unfrozen  sand  is  also  very 

similar  for  values  of  <  5  MPa,  as  shown  in  Figure  5.2. 

•  Increasing  confinement  causes  similar  increases  in  peak  strength  up  to  5 
MPa  for  both  types  of  testing  (also  see  Sheet  DIO),  in  spite  of  very  different 
values  of  the  axial  strain  at  failure.  However,  further  confinement 

increases  Qp  for  the  conventional  tests,  but  not  for  the  CF  tests. 

The  volumetric  strain  (e^)  data  from  the  consolidate-freeze  tests  are 
disappointing,  due  to  the  lack  of  results  at  high  strains  caused,  in  part,  by  the  poor  to 
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very  poor  stabilities  exhibited  by  nearly  all  of  the  tests  (see  Sheet  D3).  However,  the 
available  results  do  show  a  decrease  in  the  amount  of  dilation  with  increasing 
confinement,  as  should  be  expected  (Figure  5.12  and  Sheet  D9).  In  addition,  the  CF 
tests  may  cause  a  slight  increase  in  the  amount  of  dilation  relative  to  conventional 
tests. 

5.4.3  Summary  and  Discussion 

The  consolidate-freeze  (CF)  test  program  presents  the  first  known 
attempt  to  measure  the  influence  on  frozen  sand  behavior  of  the  magnitude  of  the 
effective  stresses  acting  on  the  sand  skeleton  prior  to  (and  during)  freezing.  On  the 
other  hand,  the  results  from  the  program  have  definite  limitations,  due  to  possible 
expansion  of  the  sand  specimens  during  freezing  (hence,  the  possibility  of  non- 
uniform,  ice-rich  samples),  the  lack  of  data  at  very  large  strains  and  generally  poor  to 
very  poor  large  strain  specimen  stability,  and  restriction  to  testing  at  one  strain  rate 
and  temperature.  With  these  limitations  in  mind,  increases  in  the  pre-freezing 

effective  stress  level  from  =  0.2  to  2,  5  and  10  MPa  appear  to  cause  the  following 
behavior  trends,  compared  to  conventional  tests  with  similar  levels  of  confinement 
run  on  specimens  prepared  by  multiple  sieve  pluviation  (and,  hence,  having 
essentially  zero  effective  stress  during  freezing  and  also  prior  to  shear  independent 
of  the  confining  pressure): 

•  no  change  in  stress-strain  behavior  at  very  small  strains,  i.e.,  constant 
values  of  Young's  modulus  (E)  and  yield  offset  stress  (Qyo)  (Sheet  D8),  as 
should  be  expected; 
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•  perhaps  a  slight  increase  in  the  amount  of  dilation  (Figure  5.12); 

•  a  disappearance  of  the  upper  yield  point  (except  for  the  two  CF  tests  with 

a’c  <  0.3  MPa),  which  is  thought  to  present  the  onset  of  significant  cracking 
of  the  pore  ice  matrix  from  detailed  analysis  of  data  from  conventional 
tests  on  frozen  MFS  (and  from  the  literature); 

•  substantial  increases  in  the  stiffness  within  the  upper  yield  region  (Figure 
5.11),  which  apparently  reflects  contributions  from  the  frictional  resistance 
of  the  sand  skeleton  based  on  the  measured  undrained  shear  behavior  of 
unfrozen  sand  (e.g..  Figure  5.2);  and 

•  the  same  increase  in  peak  strength  (Qp)  up  to  5  MPa,  but  with  very 
different  strains  to  failure  and  post-peak  rates  of  strain  softening  (Figure 
5.12  and  Sheet  DIO). 

In  retrospect,  the  CF  test  program  should  have  included  more  testing  at 
consolidation  stresses  less  than  2  MPa,  in  order  to  better  define  the  disappearance  of 
the  upper  yield  point  and  the  presumed  influence  of  sand  skeleton  frictional 
resistance  within  the  upper  yield  region. 

The  consolidate-freeze  test  program  was  expected  to  affect  frozen  sand 
behavior  to  a  much  greater  extent  than  observed  from  comparisons  of  conventional 
tests  run  on  specimens  prepared  by  multiple  sieve  pluviation  (MSP)  versus  wet 
tamping  (WT),  wherein  the  latter  preparation  technique  presumably  causes  locked- 
in  pre-freezing  effective  stresses  due  to  compaction  (e.g..  Figure  5.10).  Figure  5.13 
compares  CF  tests  run  on  MSP  specimens  with  conventional  tests  run  on  WT 

specimens  of  dense  MFS  having  confining  stresses  of  0.2  ±  0.1,  2  and  5  MPa.  This 
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comparison  suggests  that  compaction  to  a  dense  state  produces  locked-in  stresses  of  a 
magnitude  equivalent  to  about  2  MPa.  Nevertheless,  all  compacted  specimens 
exhibited  a  distinct  yield  stress  with  values  of  Qyy  that  closely  approximated  those 
for  conventional  tests  run  on  pluviated  specimens! 

5.5  APPLICATION  OF  DILATANCY-HARDENING  MODEL  TO  FROZEN  MFS 
5.5.1  Selection  of  Model  Parameters 

From  Section  5.2,  Ladanyi’s  dilatancy-hardening  model  predicts  that 

the  peak  compressive  strength  of  frozen  sand  (Qps)  has  two  components: 

Qps^Qfs'^'Qi  . [5-^] 

where  Qfs  is  the  peak  compressive  strength  of  the  sand  skeleton,  and  Q,  is  the  peak 
compressive  strength  of  the  ice  matrix.  In  turn,  the  frictional  resistance  of  the  sand 
skeleton  at  high  confinement  leading  to  the  non-cavitation  regime  is  given  by 

Qfs  =  CJ‘cr(Rf-l) . [5.7a] 

where  a'^r  is  the  critical  confining  pressure  (=  a’3f  for  zero  excess  pore  pressure  at 
failure  for  unfrozen  sand)  and  Rf  =  (o'i/a'3)f  =  (1  +  sin  <{)')/ (1  -  sin  ([)').  At  low 

confinement  leading  to  cavitation  (i.e.,  for  +  Tj  <  o'er  producing  tensile  failure 
within  the  ice  matrix),  the  sand  skeleton  compressive  strength  is  given  by 

Qfs  =  (<^c+Ti)(Rf-l) . [5.7b] 

where  Oc  =  03  is  the  applied  confining  pressure  and  Tj  is  the  tensile  strength  of  the 
ice  matrix.  This  relationship  assumes  that  a'3f  =  a<;  +  Tj  within  the  sand  skeleton  at 
failure  of  the  frozen  sand. 
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Unfrozen  Sand  Parameters 


Figure  5.7  showed  the  effective  stress  envelope  from  consolidated- 
undrained  triaxial  compression  tests  at  maximum  obliquity  for  a  wide  range  of  sand 
densities  and  consolidation  stresses,  leading  to  the  selection  of  (])'  =  34.9°.  This 

friction  angle  gives  Rf  - 1  =  2.675.  Changing  (()'  by  10%  changes  the  calculated  by 
about  20%. 

Figure  5.8  showed  the  critical  confining  pressure  line  that  was  obtained 
by  determining  the  state  parameter  leading  to  zero  excess  pore  pressure  at  failure 
(i.e.,  Af  =  0  in  Figure  5.9).  Values  of  o’er  were  then  calculated  for  each  test  using 

log  o'er  (MPa)  =  (0-699  -  e)/0.254 . [5.8] 

This  relationship  and  Rf  -  1  =  2.675  predict  the  following  values  of  Qfs  for  the  non- 
cavitating  case: 


Density 

D,  (%) 

o’er  (MPa) 

Qfc  (MPa) 

Loose 

35 

0.423 

1.13 

Medium 

65 

1.035 

2.77 

Dense 

95 

2.52 

6.74 

Note  that  these  values  of  Qf^  apply  to  all  confining  pressures,  strain  rates,  and 

temperatures,  as  long  as  o'^r  is  less  than  o^  +  Tj. 

Pore  Ice  Parameters 

Sheets  El  to  E3  in  Appendix  E  indicate  that  the  tensile  strength  of 
polycrystalline  ice  exhibits  relatively  little  change  with  strain  rate  and  temperature 
over  the  range  of  conditions  used  for  the  frozen  MFS  test  program,  i.e., 

e  =  3  X  10^/sec  to  5  X  10^/ sec  and  T  =  -10°C  to  -25°C.  The  data  in  Sheet  E3  show 
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values  of  Tj  typically  ranging  between  about  1.7  and  2.7  MPa.  The  authors  selected 

Tj  =  2.15  MPa  for  analysis  of  frozen  MFS.  It  is  emphasized  that  variations  in  Tj  of  ± 
0.5  MPa  will  not  affect  the  predictions,  except  for  dense  sand  at  low  confinement 
where  may  exceed  +  Tj  (i.e.,  leading  to  the  cavitation  regime). 

The  literature  does  not  provide  data  on  the  peak  compressive  strength 
(Qi)  of  polycrystalline  ice  for  the  range  of  confining  pressures,  strain  rates  and 
temperatures  used  for  the  frozen  MFS  test  program,  which  also  corresponds  to  the 
ductile-to-brittle  transition  regime  for  granular  ice.  Consequently,  values  of  Qj  were 
predicted  as  follows: 

•  From  the  extensive  test  program  by  Jones  (1982)  at  an  average  temperature 

of  T  « -11.5  ±  1°C,  interpolate  values  of  Q,  at  Oc  =  0  10  MPa  for  strain 

-6  -2 

rates  varying  from  e  =  10  /sec  to  10’  /sec.  These  values  are  plotted  on 
Sheet  E4  and  lead  to  n  =  5.43  and  4.26  at  low  and  high  confinement, 
respectively.  Note:  Sheet  E5  indicates  that  Jones'  data  at  Oj.  =  0  agree 

reasonably  well  with  results  from  a  variety  of  unconfined  creep  and 
constant  strain  rate  test  programs  on  isotropic  polycrystalline  ice. 

•  Use  the  Arrhenius  equation  combined  with  Glen's  (1955)  power  law 
equation  (shown  at  the  top  of  Sheet  E4)  to  fit  Jones'  data  for  an  assumed 
activation  energy  equal  to  68  kj/mole.  This  leads  to  the  following 
relationships: 


125 


predicted  peak  strength,  Qps-  Note:  linear  interpolation  was  used  to 
obtain  values  of  Q,  for  confining  pressures  between  0.1  and  10  MPa. 

Effect  of  Confinement  for  Dense  Sand 
Both  conventional  and  consolidate-freeze  tests  were  run  at  varying 
levels  of  confinement  on  dense  MFS  at  the  moderate  strain  rate  and  T  =  -10°C. 
Figure  5.14  compares  these  data  with  predicted  peak  strengths.  Note  that  the  values 
of  Qp  for  the  conventional  tests  were  obtained  from  linear  regression  on  Qp  versus 

Df  data  at  each  confining  pressure.  At  low  confinement  =  o’c  <  0.25  MPa)  there  is 
moderate  agreement,  since  Qps  is  only  about  15%  to  20%  less  than  the  mean 

measured  strengths.  For  the  assumed  tensile  strength  of  ice  (Tj  =  2.15  MPa),  Oc  +  Tj 
usually  is  greater  than  the  critical  confining  pressure;  hence,  the  non-cavitation 
Eq.  5.7a  usually  applies.  This  leads  to  Qfg  ~  6  ±  1  MPa,  compared  to  Qj  =  4  MPa. 
However,  the  agreement  may  be  fortuitous,  due  to  severe  problems  with  strain 
compatibility  since  the  ice  matrix  presumably  reaches  its  peak  strength  at  around  1% 
axial  strain;  the  sand  skeleton  requires  over  20%  strain  based  on  the  unfrozen  sand 
data  in  Figure  5.2;  and  the  frozen  specimens  failed  at  about  6%  axial  strain. 

Figure  5.14  shows  that  increasing  confinement  produces  a  much  larger 
increase  in  measured  peak  strengths  than  predicted  by  the  dilatancy-hardening 
model,  e.g.,  predicted  strengths  about  two  thirds  of  the  measured  strengths  at  5  MPa 

confinement.  The  predicted  strengthening  occurs  solely  from  increases  in  Qj  due  to 
added  confinement,  since  Qfs  remains  constant  within  the  non-cavitation  regime. 
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Qi  (MPa)  =  0.0897  for  Cc  =  0 . [5.9a] 

Qi  (MPa)  =  0.0415  for  =  10  MPa . [5.9b] 

where  T  is  in  degrees  Kelvin. 

•  Use  Eq.  5.9  to  calculate  values  of  Qj  as  a  function  of  e  and  T  at  low  and 
high  confinements.  These  values  are  plotted  on  Sheet  E6,  which  also  plots 
data  from  Murrell  et  al.  (1989).  There  is  reasonable  agreement  between 
predicted  and  measured  ice  strengths  at  T  =  -20°C  at  low  and  high 
confinement. 

5.5.2  Results  of  Analyses 

The  analyses  proceeded  as  follows: 

1)  Determine  a'^r  for  each  conventional  frozen  and  consolidate-freeze 
MFS  test,  based  on  the  specimen  (void  ratio). 

2)  Based  on  the  values  of  a'^r,  Tj  =  2.15  MPa  and  the  actual  confining 
pressure  (Oc «  0.1  to  10  MPa),  determine  if  the  specimen  falls  within 

the  non-cavitation  or  cavitation  regimes  for  computation  of  Qfg 
using  Eq.  5.7.  The  non-cavitation  case  applied  for  all  tests  except 
dense  sand  (D^  >  90%)  at  =  0.1  MPa. 

3)  Calculate  Qj  for  each  test  using  Eq.  5.9  based  on  the  nominal  and 
T  and  the  actual  strain  rate;  add  this  value  to  Qfg  to  obtain  the 
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Strain  compatibility  also  remains  problematic,  given  the  following  trends  in  the 
axial  strain  at  the  peak  strength  (Ep): 

•  for  the  conventional  tests,  Ep  increases  to  more  than  20%,  whereas  Ep  for 
the  CF  tests  decreases  slightly  (see  Figure  5.12);  and 

•  for  unfrozen  sand,  Ep  undergoes  a  large  decrease  over  the  range  of  a'c  =  0.1 

to  5  MPa  during  undrained  shear  (Figure  5.2). 

In  any  case,  it  would  appear  that  added  confinement  causes  a  much  greater 
contribution  from  the  frictional  resistance  of  the  sand  skeleton  than  predicted  by  the 
model.  Confinement  may  also  produce  greater  strengthening  of  the  ice  matrix  than 
assumed  by  the  authors.  However,  uncertainty  in  the  tensile  strength  of  ice  should 
not  be  a  factor,  except  possibly  within  the  region  of  very  low  confinement. 

Note:  Sheet  Ell  plots  predicted  and  measured  peak  strengths  versus  relative 
density  from  conventional  tests  having  confining  pressures  of  0.1,  2,  5,  and  10  MPa 
for  the  moderate  strain  rate  and  T  =  -10°C.  In  contrast  to  dense  sand,  the  relative 
difference  between  predicted  and  measured  strengths  for  loose  sand  does  not 
increase  with  increasing  levels  of  confinement.  The  subsequent  text  discusses  loose 
versus  dense  sand  comparisons  in  detail. 

Effect  of  Strain  Rate  and  Temperature  for  Dense  and  Loose  Sand 

Figures  5.15  and  5.16  compare  measured  and  predicted  peak  strengths 

as  a  function  of  strain  rate  and  temperatvure  for  dense  MFS  at  a^.  =  0.1  and  10  MPa, 

respectively.  The  measmed  strengths  are  from  actual  tests  (not  from  linear 
regression);  hence,  they  reflect  scatter  due  to  variations  in  the  testing  conditions 
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(especially  relative  density).  The  predictions  are  for  values  of  D^,  e  and  T 
corresponding  to  each  test.  At  low  confinement,  Figure  5.15  shows  the  following 
trends: 

•  At  T  =  -10°C,  there  is  excellent  agreement  (within  10%)  between  predicted 
and  measvued  strengths  at  the  slow  strain  rate.  However,  the  predictions 
become  less  satisfactory  with  increasing  strain  rate.  This  results  from  the 
fact  that  the  model  predicts  much  higher  values  of  the  power  law 
coefficient  (n)  than  measured  and,  hence,  much  less  rate  sensitivity. 

•  At  T  =  -20°C,  the  predicted  strengths  are  much  less  than  measured  (about 
one  half)  and  the  model  again  underestimates  the  rate  sensitivity. 

•  At  T  =  -25°C,  the  discrepancy  in  strengths  becomes  even  worse,  but  with 
similar  rate  sensitivities  (but  based  on  only  two  tests). 

Collectively,  the  model  predictions  become  less  satisfactory  at  higher  strain  rates  and 
at  lower  temperatures. 

For  dense  MFS  at  high  confinement,  the  results  in  Figure  5.16  at 
T  =  -10  and  -15°C  show  the  following  trends: 

•  The  predicted  strengths  are  again  too  low  and  also  become  less  satisfactory 
at  the  lower  temperature. 

•  However,  increases  in  the  strain  rate  now  improve  the  predictions,  which 
is  opposite  to  the  results  at  low  confinement 

Figures  5.17  and  5.18  provide  similar  comparisons  for  loose  MFS  at  low 
and  high  confinement,  respectively;  they  show  the  following  trends: 
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•  At  Cc  =  0.1  MPa,  the  predicted  strengths  are  much  too  low  and  again 
become  less  satisfactory  at  the  lower  temperatures,  but  with  little  influence 
of  strain  rate  (i.e.,  similar  values  of  n) 

•  At  Oc  =  10  MPa,  the  model  greatly  overpredicts  the  rate  sensitivity  (i.e., 
values  of  n  that  are  too  low);  hence,  increases  in  strain  rate  lead  to  better 
agreement. 

It  should  be  emphasized  that  the  predicted  contribution  of  the  sand  skeleton  friction 
for  loose  sand  (Dj  =  35%)  is  only  Qfs  =  1  MPa  at  both  levels  of  confinement  (since  the 
sand  is  well  within  the  non-cavitation  regime).  Hence,  the  ice  matrix  strength  (Qj) 

dominates  the  predicted  values  of  Qps- 

5.5.3  Collective  Comparison  and  Revised  Model  Predictions 
Ratio  of  Predicted  to  Measured  Peak  Strengths,  Qps/Qp 
Sheet  E8  in  Appendix  E  plots  Qps/Qp  versus  temperature  for  dense 
(Dr  =  95%)  and  loose  (D^  =  35%)  MFS.  Qps  equals  Qfs  +  Qj  from  Ladanyi's  model,  and 
Qp  is  the  peak  strength  from  conventional  frozen  tests.  The  Qps/ Qp  ratios  were 
scaled  from  Figures  5.15  to  5.18  at  the  slow,  moderate  and  fast  strain  rates.  For  dense 
frozen  MFS,  the  Qp  values  range  between  about  10  and  25  MPa,  compared  to  a 
predicted  Qfg  of  about  6  ±  1  MPa,  and  Qps/Qp  decreases  from  roughly  0.7  ±  0.1  at 
T  =  -10°C  down  to  0.5  at  T  =  -25°C.  For  loose  frozen  MFS,  the  Qp  values  range 
between  about  7  and  30  MPa,  compared  to  a  predicted  Qfs  of  only  about  1.0  to  1.5 
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MPa,  and  Qfs/Qp  decreases  from  roughly  0.6  ±  1  at  T  =  -10°C  down  to  0.45  ±  0.1  at 
T  =  -25°C.  These  results  will  be  discussed  further  after  evaluating  the  Q,  (ice  matrix 
strength)  component  of  the  model  predictions. 

Ratio  of  Predicted  Ice  Strength  to  Measured  Upper  Yield  Stress,  Qi/Quy 

As  demonstrated  in  Section  3.4,  the  upper  yield  stress  (Quy)  of  frozen 
MFS  occurs  at  low  strains  (Eg  =  0.3  to  1.1%)  and  is  not  affected  by  sand  density  (D^) 
and  confining  pressure  {of).  If  one  assumes  zero  contribution  from  sand  skeleton 

friction  to  Quy,  which  appears  reasonable  for  conventional  frozen  tests  (i.e.,  very  low 
effective  stresses  from  pluviation  and  no  increase  due  to  application  of  the  cell 
pressure),  then  Quy  may  represent  the  operational  or  effective  strength  of  the  ice 

matrix  prior  to  the  onset  of  significant  cracking.  Assuming  that  the  effective  Qj  is 
equal  to  the  measured  upper  yield  stress,  Qj  (eff)  =  Quy,  how  do  these  values  compare 
with  the  values  of  Qi  used  for  the  model  predictioris?  Sheet  E9  in  Appendix  E  plots 

the  ratio  of  the  predicted  Qj  to  the  measured  Quy  versus  temperature  as  a  function  of 
strain  rate  for  zero  and  high  (10  MPa)  confinement.  One  observes  the  following 
trends: 

•  At  Oc  *=  0,  Qj/Quy  decreases  from  about  0.5  at  T  =  -10°C  to  about  0.3  at 
T  =  -25°C  with  relatively  little  effect  of  strain  rate. 

•  At  Oc  =  10  MPa,  Qi/Quy  decreases  from  about  0.65  to  roughly 

0.45  ±  0.1  over  the  same  temperature  range.  In  addition,  the  ratio  is 
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significantly  higher  in  going  from  the  slow  to  the  fast  strain  rate  at  the 
three  lower  temperatures. 

These  trends  suggest  that  data  in  the  literature  from  tests  on  polycrystalline  ice  could 
significantly  underestimate  the  strength  of  the  ice  matrix  in  frozen  sand,  and  that 
the  error  becomes  progressively  worse  with  decreasing  temperature. 

Comparison  of  Trends  in  Qps/Qp  Q/Quy 
The  comparison  is  made  separately  at  low  and  high  levels  of  confining 


pressure. 
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If  the  effective  strength  of  the  ice  matrix  can  be  approximated  by  Q^y,  then  the  above 
comparison  offers  a  reasonable  (but  still  tentative)  hypothesis  to  help  explain  the 
following  key  aspects  of  the  model  predictions: 

•  the  fact  that  Qps/Qp  decreases  with  decreasing  temperature,  as  a  general 
rule,  independent  of  sand  density  and  coiafining  pressure;  and 
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•  the  fact  that  Qps/Qp  increases  with  increasing  strain  rate  at  high 

confinement  for  both  dense  and  loose  sand  (as  per  Qi/Quy)/  whereas  the 
opposite  occurs  at  low  confinement.  (Qi/Quy  ^iso  shows  the  opposite,  but 
only  at  -10°C.) 

Revised  Model  Predictions 

In  order  to  further  check  the  hypothesis  that  the  effective  strength  of 
the  ice  matrix  can  be  approximated  by  Qi(eff.)  =  Quy/  revised  model  predictions  were 
made  using  the  relationship 

QFS*  =  Qfs  +  Quy . [5-10] 

where  Qps*  denotes  the  new  predicted  strength,  Qfg  is  calculated  as  before,  and  Quy 
replaces  Qi-  Values  of  Qfg  were  calculated  for  loose  and  dense  sand,  giving  1.15  MPa 
at  Dj  =  35%  and  6.75  MPa  at  Dj.  =  95%.  Note:  the  latter  strength  becomes  6.0  MPa  at 
CTc  =  0.1  for  the  previously  assumed  Ti  =  2.15  MPa,  but  this  difference  is  ignored, 
given  the  scatter  in  Ti  as  per  Sheet  E3.  Mean  values  of  Quy  were  used  at  each  strain 
rate  and  temperature,  and  values  of  Qp  were  obtained  for  comparison  from  linear 

regression  on  the  measured  Qp  versus  data  for  tests  having  Type  C  and  D  curves. 

Sheet  ElO  compares  the  difference  in  predicted  and  measured  strengths 
for  the  revised  model.  It  plots  Qps*  -  Qp  versus  temperature  for  dense  and  loose 

MFS  at  Oc  =  0.1  and  10  MPa.  As  a  general  rule,  the  predicted  Qps*  is  either  quite 
accurate  or  slightly  too  high  at  low  confinement;  whereas  at  high  confinement,  the 
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predicted  strength  is  too  low.  For  added  perspective,  refer  to  Figure  4.6  for  measured 
values  of  Qp  -  Quy  plotted  versus  relative  density.  In  essence,  Qps*  -  Qp  is  the  same  as 
comparing  Qfg  -  (Qp  -  Quy)-  Hence,  errors  in  Qps*  -  Qp  will  approach  Qfg  as  Qp 
approaches  Quy. 

The  following  comparisons  focus  on  test  conditions  with  Qp 
significantly  larger  than  Quy/  i-G-/  having  "solid"  Type  C  and  D  curves. 

Comparison  for  Dense  MFS  at  -  0.1  MPa 

•  Qp  -  Qyy  =  4  MPa  for  slow  and  moderate  ^  at  all  temperatures;  £p  =  6  ±  1%; 
fast  e  produces  Type  A  curves  at  T  <  -10°C. 

•  Predicted  Qps*  is  too  high  by  about  2.5  MPa,  which  represents  about  one 
third  of  the  predicted  Qfg  =  6.75  MPa. 

•  Commentary:  the  Ep  is  much  too  small  to  expect  full  mobilization  of  the 
sand  friction  (compared  to  undrained  shear  of  unfrozen  sand;  see  figure 
5.2),  which  is  consistent  with  the  observed  overprediction.  However,  bulk 

ice  would  undergo  significant  strain  softening  before  reaching  Ep  (e.g.. 
Figure  4.22).  Hence,  the  overall  agreement  is  better  than  expected,  based 
on  separate  evaluation  of  the  behavior  of  unfrozen  sand  and  bulk  ice. 

Comparison  for  Dense  MFS  at  =  10  MPa 

•  Qp  -  Qyy  =  10  ±  1  MPa  for  slow  and  moderate  e  at  T  >  -20°C,  Ep  =  22%  at 
slow  6,  -10°C  and  decreases  to  =  10%  at  higher  6  and  lower  T. 
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•  Predicted  Qps*  is  too  low  by  slightly  over  2  MPa,  or  about  double  the 
predicted  Qfg  =  1.15  MPa. 

•  Commentary:  the  Ep  level  should  be  sufficient  to  fully  mobilize  the  sand 

friction  (although  no  direct  comparison  can  be  made  to  unfrozen  loose 
sand  at  high  confinement).  One  would  expect  some  weakening  of  the  ice 
matrix  at  these  large  strain  levels.  Hence,  the  measured  strengths  are 
larger  than  expected  from  consideration  of  the  separate  behaviors. 

5.5.4  Summary  and  Conclusions 

Ladanyi's  dilatancy  hardening  model  underpredicts  the  measured  peak 
strength  of  frozen  Manchester  Fine  Sand  (MFS),  as  measured  in  conventional  tests 

on  loose  and  dense  sand  (D^  =  35  and  95%)  at  low  and  high  confinement 

(0(,  =  0.1  and  10  MPa),  as  illustrated  in  Figures  5.15  to  5.18  and  Sheet  E8  (Qpg  is  the 
strength  predicted  by  his  model).  The  magnitude  of  the  error  increases  with 
decreasing  temperature  (a  pervasive  trend).  For  dense  sand  sheared  at  the  moderate 
strain  rate  and  T  =  -10°C,  the  error  also  increases  with  confinement  for  both 
conventional  and  consolidate-freeze  tests  (Figure  5.14). 

The  above  predictions  used  values  of  strength  for  the  ice  matrix  (Qj) 
based  on  data  in  the  literature  from  testing  polycrystalline  ice  having  very  large 
grain  sizes  (>  1  mm),  compared  to  the  size  of  the  pores  in  MFS.  It  is  hypothesized 
that  the  operational  or  effective  peak  strength  of  the  ice  matrix  can  be  approximated 
by  the  measured  upper  yield  stress  (Quy)  of  frozen  MFS,  since 
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•  Predicted  Qps"^  is  too  low  by  3  ±  1  MPa,  which  represents  about  one  half  of 
the  predicted  Qfg  =  6.75  MPa. 

•  Commentary:  for  ep  =  20  ±  5%,  one  might  expect  near-full  mobilization  of 
the  sand  friction  (again  see  Figure  5.2).  However,  bulk  ice  would  still 
undergo  some  strain  softening,  in  spite  of  the  high  confinement  (e.g.. 
Figure  4.22  plus  Sheet  E7).  Hence,  the  large  strain  interaction  of  the  sand- 
ice  composite  produces  strengths  greater  than  expected  from  separate 
behaviors,  even  after  incorporating  a  much  higher  peak  strength  for  the 
ice  matrix. 

Comparison  for  Loose  MFS  at  =  0.1  MPa 

•  Qp  -  Qyy  « 1  MPa  for  slow  e  at  all  temperatures  and  moderate  e  at  -10°C; 
ep«6±2%. 

•  Predicted  Qps*  is  in  excellent  agreement  (for  slow  e  tests). 

•  Commentary:  the  Ep  level  is  too  small  to  expect  full  mobilization  of  the 
sand  friction  (compared  to  undrained  shear  of  unfrozen  sand;  see  Figure 
5.5).  Moreover,  bulk  ice  would  again  undergo  significant  weakening  at 

these  Ep  levels.  Hence,  there  is  better  agreement  than  expected  from 
evaluation  of  the  separate  behaviors. 

Comparison  for  Loose  MFS  at  Cc  =  10  MPa 

•  Qp  -  Qyy  ==  3.5  for  slow  at  T  >  -20°C  and  moderate  ^  at  -10°C; 

Ep  »  25-20%. 
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•  Quy  occurs  at  strain  levels  (Eg  ^  1%)  consistant  with  the  peak  strength  of 
bulk  polycrystalline  ice;  and 

•  because  Q„v  is  not  affected  by  sand  density  and  confining  pressure  (i.e.,  by 

uy 

factors  known  to  affect  the  frictional  resistance  of  unfrozen  sand),  it 
probably  involves  very  little  frictional  resistance  of  the  sand  skeleton. 

If  the  above  hypothesis  is  true,  then  Sheet  E9  shows  that  Q^/ Q^y,  where 
Quy  =  Qi  (effective),  has  the  following  trends: 

•  Qi/Quv  decreases  substantially  with  decreasing  temperature.  The  same 

1  uy 

trend  is  observed  for  Qps/ Qp/  where  Qpg  Qp  model  predicted 

and  measured  strengths,  respectively  (see  Sheet  E8). 

•  Qi/Quv  increases  with  increasing  strain  rate  at  high  confinement;  the  same 

1  uy 

trend  is  observed  for  Qps/Qp-  Note:  at  low  confinement,  Qj/Quy  shows 
little  effect  of  strain  rate,  whereas  Qps/Qp  decreases  with  increasing  strain 
rate. 

The  authors  can  readily  envision  a  higher  effective  strength  for  the  ice  matrix  than 
measured  from  tests  on  bulk  ice  due  to  the  higher  (effective)  strain  rate,  added 
deformational  constraints,  etc.,  as  hypothesized  by  Ting  et  al  (1983).  However,  they 

can  not  offer  any  explanation  for  why  Qj/ Quy  should  undergo  such  a  large  decrease 
with  decreasing  temperature. 

Revised  model  predictions  were  made  by  replacing  Qj  with  Qj(effective) 
=  Quy,  leading  to  Qps*  =  Qf*  +  Quy,  where  Qfg  is  the  predicted  sand  skeleton  friction. 
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Sheet  ElO  plots  the  predicted  Qfs*  minus  the  measured  peak  strength  (Qp)  versus 
temperature  for  loose  (Dj.  =  35%)  and  dense  (Dj.  =  95%)  MFS  at  low  and  high 
confinement  (Oj,  =  0.1  and  10  MPa)  for  tests  having  Type  C  and  D  curves  (i.e., 

Qp  >  Quy).  This  comparison  generally  shows  that  the  predicted  is  either  quite 
accurate  or  slightly  too  high  at  low  confinement,  whereas  the  predicted  strength  is 
too  low  at  high  confinement  (for  strain  rates  and  temperatures  producing 
maximum  values  of  Qp  -  Quy)-  After  considerations  of  strain  compatibility  (e.g., 

insufficient  strains  to  mobilize  the  full  sand  friction  (Qfs)/  but  excessive  strains  that 
produce  weakening  of  the  ice  matrix),  it  is  tentatively  concluded  that  the  measured 
strengths  are  larger  than  would  be  expected,  based  on  the  separate  behaviors  of  an 
unfrozen  sand  skeleton  and  an  enhanced  ice  matrix  strength  that  follows  trends 
observed  for  bulk  ice. 

In  spite  of  the  above  limitations  of  the  revised  model  predictions 
(which  use  Qyy  as  a  base  strength),  this  model  does  provide  reasonable  estimates  of 

the  peak  strength  of  frozen  MFS  over  a  very  wide  range  of  sand  densities,  confining 
pressures,  strain  rates  and  temperatures.  However,  larger  issues  still  remain,  such 
as 

•  an  explanation  for  why  the  ice  matrix  in  frozen  sands  appears  to  be  much 
larger  than  predicted  from  tests  on  bulk  ice  (if,  indeed,  the  effective  Qj 

equals  Quy), 

•  a  methodology  for  predicting  Quy;  and 
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a  methodology  for  predicting  the  absence  of  sand  friction,  i.e.,  when  the 
peak  strength  (Qp)  approaches  the  upper  yield  stress  (Quy)- 
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Fig.  5.3  Effect  of  Consolidation  Stress  on  Effective  Stress  Paths  from  CU 
Triaxial  Compression  Tests  on  Unfrozen  Dense  MFS 
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Fig.  5.4  Effect  of  Consolidation  Stress  on  Drained  Stress  — Strain  and 

Volumetric  Strain  Behavior  from  CD  Triaxial  Compression  Tests 
on  Unfrozen  Dense  MFS 
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Fig.  5.5  Effect  of  Relative  Density  on  Undrained  Stress-Strain  and 
Pore  Pressure  Behavior  from  CU  Triaxial  Compression  Tests 
at  Low  Confinement  on  Unfrozen  MFS 
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Fig.  5.6  Effect  of  Relative  Density  on  Undrained  Stress-Strain  and 
Pore  Pressure  Behavior  from  CU  Triaxial  Compression  Tests 
at  Moderate  Confinement  on  Unfrozen  MFS 
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Triaxial  Compression  Tests  on  Unfrozen  MFS 
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Fig.  5.10  Effect  of  Specimen  Preparation  Technique  on  Stress-Strain 
Curves  from  Conventional  Tests  on  Frozen  Dense  MFS 
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Moderate 


Consolidate-Freeze  and  Conventional  Tests  on  Dense  MFS  at 
Moderate  Strain  Rate  and  T  =  —  10°C 
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Figure  5.12  Comparison  of  Stress-Strain  Curves  for  Consolidate 
Freeze  and  Conventional  Frozen  Tests  on  Dense 
MFS  at  Moderate  Strain  Rate  and  T  =  -10°C 
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Figure  5.13  Comparison  of  Stress-Strain  Curves  for  Consolidate- 

Freeze  Tests  on  MSP  Specimens  and  Conventional  Frozen 
Tests  on  Wet  Tamped  Specimens  at  Moderate  Strain  Rate 
and  T  =  -10°C 
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6.  SUMMARY  AND  CONCLUSIONS 


6.1  RESEARCH  OBJECTIVES  AND  SCOPE 

Even  the  simplest  form  of  frozen  soil  (sand  with  ice-filled  pores)  entails  a 
highly  complex  interaction  between  the  pore  matrix  (composed  of  ice  and  unfrozen 
water)  and  the  skeleton  of  sand  particles  that  changes  with  time  as  a  function  of 
temperature  and  stress-strain  level.  Consequently,  the  development  of  constitutive 
relationships  to  model  the  stress-strain-time-temperature  behavior  of  frozen  sands 
lags  far  behind  the  modeling  capabilities  of  the  two  principal  components,  ice  and 
sand.  Moreover,  progress  in  modeling  first  requires  a  better  understanding  of  the 
physical  mechanisms  that  control  frozen  sand  behavior.  The  approach  taken  by 
MIT  to  achieve  this  latter  objective  has  three  components: 

1)  conduct  a  comprehensive  experimental  program  to  precisely  measure  the 
behavior  of  a  frozen  natural  sand  over  a  wide  range  of  testing  conditions; 

2)  from  a  knowledge  of  the  general  behavior  of  polycrystalline  ice  and  of  the 
measured  strain-strain  behavior  of  the  same  sand  in  an  unfrozen  state, 
attempt  to  deduce  the  relative  importance  of  the 

•  ice  matrix  per  se 

•  frictional  resistance  of  the  sand  skeleton  per  se,  and 

•  interaction  between  the  ice  matrix  and  the  sand  skeleton; 

3)  where  possible,  employ  models  for  composite  materials  to  quantify  the 
physical  mechanisms  responsible  for  the  macro-behavior  of  frozen  sands. 
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Manchester  Fine  Sand  (MFS)  was  selected  for  the  test  material.  This  river 
deposit  is  composed  mainly  of  subangular  quartz  and  feldspar  particles  having  the 
following  characteristics: 

•  gradation — 99%  smaller  than  0.35  mm,  mean  diameter  =  0.18  mm,  and  7% 
smaller  than  0.074  mm  (No.  200  sieve). 

•  density— minimum  dry  density  =  1408  kg/m^  (e^ax  =  0-909)  and  maximum 

dry  density  =  1701  kg/m^  (en,in=  0.580).  The  report  presents  density  values 

in  terms  of  relative  density  defined  asDr=  (Cmax"  s)/  (®max"®min)'  where  e  is 
the  void  ratio  (volume  of  voids  divided  by  volume  of  solids). 

•  unfrozen  water  content — ^less  than  0.14%  at  T  =  -9.5°C. 

All  test  specimens  were  prepared  by  multiple  sieve  pluviation  which  gives  a  very 
uniform  density,  in  contrast  to  many  prior  programs  that  employed  various  types  of 
mechanical  compaction. 

The  experimental  program  consisted  of  triaxial  compression  tests  performed 
in  a  high-pressure  triaxial  apparatus  on  three  types  of  specimens:  1)  specimens  that 
were  frozen  in  a  mold  before  being  placed  in  the  triaxial  cell,  called  conventional 
frozen  tests;  2)  specimens  that  were  frozen  after  being  consolidated  to  varying 
effective  stress  levels  in  the  triaxial  cell,  called  consolidate-freeze  tests;  and  3) 
specimens  that  were  not  frozen,  i.e.,  triaxial  tests  on  unfrozen  sand.  The 
experimental  procedures  used  for  testing  frozen  sand  had  three  special  features, 
which  collectively  make  the  results  unique  compared  to  prior  data.  These  features 
are 
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•  measurement  of  on-specimen  axial  strains  (Ea)  using  special  yokes  that  can 

resolve  Ea  to  0.001%  and  which  eliminate  the  initially  S-shaped  stress- 
strain  curves  commonly  reported  in  the  literature; 

•  measurement  of  volumetric  strains  (Ey),  although  these  are  subject  to  an 

error  inEyOf  up  to  ±  0.2%; 

•  the  use  of  enlarged,  lubricated  (via  ice  caps  for  the  conventional  tests)  end 
platens  to  facilitate  lateral  straining  of  the  specimen  and  enable 

measurement  of  stress-strain  behavior  to  very  large  strains  (Eg  =  20-25%). 
The  experimental  procedures  also  had  two  less  desirable  features.  First,  compliance 
in  the  triaxial  apparatus  required  about  one  percent  axial  strain  before  the  test 

specimens  experienced  a  constant  strain  rate,  dEg/ dt  =  e;  however,  analysis  of  the 

data  did  use  actual  values  of  6.  Second,  there  was  a  temperature  gradient  of  about 
0.4°C  between  the  top  and  bottom  of  the  test  specimens  during  shear  that  precluded 
testing  at  higher  temperatures. 

Table  2.1  summarizes  the  testing  matrix  for  the  99  conventional  frozen  sand 
tests.  The  four  variables  in  the  parametric  study  were 

1)  the  nominal  testing  temperature,  which  varied  in  5°C  increments  from 
T  =  -10°C  to  -25°C. 

2)  the  nominal  axial  strain  rate,  which  varied  from  ^  =  3  x  10"^/ s  (slow)  to 
6  =  3.5  X  lO’^/s  (moderate)  to  6  =  5  x  10"^/ s  (fast). 
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3)  the  relative  density  of  the  sand  specimens,  which  varied  from 

Dj.  =  20  to  100%.  Except  at  T  =  -10°C,  most  specimens  were  either  loose  (D^ 

=  35%)  or  dense  (Dj  =95%). 

4)  the  confining  pressure  during  shear,  which  varied  fromo^  =  0.1  to  10  MPa. 

Testing  at  T  =  -10°C  also  included  Oc  =  2  and  5  MPa. 

This  is  the  first  set  of  data  from  conventional  frozen  tests  to  fully  characterize 
behavior  as  a  fimction  of  sand  density  (Dj.),  confining  pressure  (o^),  strain  rate  (&) 
and  temperature  (T).  In  fact,  only  one  other  program  (Baker  and  Konrad,  1985)  has 
systematically  varied  sand  density  over  all  strain  levels,  and  that  was  done  at  only 

one  value  of  o^.,  ^  and  T. 

The  program  of  consolidate-freeze  tests  consisted  of  ten  tests  run  on  dense 
sand  consolidated  to  effective  consolidation  stresses  of  o  ^  =  0.2, 2, 5  and  10  MPa, 
frozen  at  T  =  -10°C  and  then  sheared  at  the  moderate  strain  rate.  This  is  the  first 
research  to  perform  this  type  of  test  on  frozen  soil,  although  some  aspects  of  the 
freezing  process  still  require  further  verification.  The  program  of  triaxial 
compression  tests  on  unfrozen  sand  consisted  of  49  consolidated-undrained  and 

consolidated-drained  tests  covering  a  wide  variation  in  sand  density  30  to 
100%)  and  consolidation  stresses  (ct'c=  0.1  to  over  10  MPa).  This  is  the  first  research 
to  perform  tests  on  unfrozen  sand  for  direct  comparison  with  results  of  tests  on 
frozen  sand. 
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6.2  OVERVIEW  AND  DEFINITION  OF  STRESS-STRAIN  BEHAVIOR 
6.2.1  Overview  of  Stress-Strain  Behavior 

Figure  2.2  plots  deviator  stress  (Q  =Oi  -  03)  versus  engineering  axial 

strain  (£3)/  in  order  to  illustrate  the  range  in  types  of  stress-strain  curves  measured 
during  the  conventional  testing  program  on  frozen  MFS.  Two  stress-strain  curves 
are  plotted  at  each  strain  rate  and  temperature.  One  corresponds  to  spedmeiis 
having  a  low  relative  density  (Dj=  35%)  and  low  confining  pressure  (0^=  0.1  MPa) 
and,  hence,  presumably  reflecting  the  minimum  frictional  resistance  of  the  sand 
skeleton.  The  other  corresponds  to  specimens  having  a  high  relative  deirsity 
(Dr  =  95%)  and  high  confining  pressure  (a^  =  10  MPa)  and  presumably  reflecting  the 
maximum  frictional  resistance  of  the  sand  skeleton.  Figure  2.2a  shows  the  effect  of 
increasing  the  strain  rate  (e)  by  two  orders  of  magnitude  from  slow  to  fast  for 
specimens  sheared  at  T  =  -10°C;  Figure  2.2b  shows  the  effect  of  decreasing  the 
temperature  by  10°C  for  specimens  sheared  at  the  moderate  strain  rate. 

All  of  the  stress-strain  curves  in  Figure  2.2  exhibit  a  very  distinctive 
yield  point  (the  knee  of  the  curve)  at  levels  of  axial  strain  ranging  from  0.3  to  1.0%, 
as  shown  by  the  circles  in  the  figure.  This  point  represents  the  onset  of  highly  non¬ 
linear  behavior  and  the  development  of  very  significant  plastic  deformations;  it  is 
termed  the  upper  yield  point. 

A  qualitative  assessment  of  the  data  in  Figure  2.2  shows  that  the 
behavior  in  the  small-strain  region  is  largely  unaffected  by  changes  in  relative 
density  and  confining  pressure,  both  of  which  affect  the  frictional  resistance  of  the 
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sand  skeleton.  The  results  in  Figure  2.2  also  clearly  demonstrate  that  strain  rate  and 
temperature,  which  affect  ice  behavior,  have  a  tremendous  influence  on  the 
magnitude  of  the  deviator  stress  at  the  upper  yield  point,  termed  the  upper  yield 

stress  (Quy)- 

At  larger  strains  beyond  the  upper  yield  region,  the  data  in  Figure  2.2 
show  that  changes  in  relative  density  and  confining  pressure  obviously  have  a 
tremendous  influence  on  stress-strain  behavior.  These  two  variables  also  affect  the 
frictional  resistance  of  the  sand  skeleton.  Moreover,  at  large  strains,  all  of  the  test 
specimens  exhibited  varying  degrees  of  dilation  (volumetric  expansion)  that 
depended  mainly  on  sand  density  and  confiiung  pressure,  as  illustrated  by  the 
results  presented  in  Figures  4.1  to  4.3. 

6.2.2  Definition  of  Stress-Strain  Parameters 

Small-Strain  Behavior 

Figure  2.3  presents  actual  stress-strain  curves  drawn  to  two  different 
strain  scales  in  terms  of  the  deviator  stress  (Q  =0^  -03)  and  the  axial  strain  (£3).  The 
parameters  that  have  been  selected  for  the  evaluation  of  the  small-strain  behavior 
are  the  Young's  modulus,  the  yield  offset  at  0.01%  strain,  and  the  upper  yield  stress. 
The  graphical  construction  techniques  used  to  obtain  these  parameters  are 
summarized  on  the  figure  and  explained  as  follows: 

•  the  Young's  modulus  (E).  The  initial  slope  of  the  stress-strain  curve  is 
determined  visually  at  levels  of  axial  strain  less  than  0.01%. 
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•  the  yield  offset  stress  at  0.01%  strain  (Qyo).  The  intersection  point  of  the 
stress-strain  curve  and  a  line  with  a  slope  of  the  Young's  modulus 
translated  by  0.01%  strain. 

•  the  upper  yield  stress  (Quy)-  The  first  point  on  the  stress-strain  curve 

where  the  slope  (dQ/deg)  either  becomes  zero  or  reaches  its  minimum 

positive  value  before  significant  strain  hardening. 

Large-Strain  Behavior 

The  deviator  stress  versus  axial  strain  curves  in  Figure  2.2  have 
different  characteristics  regarding  the  degree  of  post-upper  yield  strain  softening  or 

strain  hardening,  i.e.,  whether  Q  decreases  or  increases  after  reaching  Q^y;  and 
whether  or  not  strain  softening  occurs  after  the  specimen  reaches  its  peak  strength 
(Qp).  These  different  characteristics,  which  depend  mainly  on  sand  density  and 
confining  pressure,  are  illustrated  in  the  upper  portion  of  Figure  2.4  for  tests 
conducted  at  low  and  high  confining  pressures  =  0.1  and  10  MPa).  Note  that  the 
deviator  stress  is  normalized  by  the  upper  yield  stress,  since  Quy  does  not  vary  with 
Dj.  andOc.  Figure  2.4  also  summarizes  characteristic  volumetric  behavior  measured 

at  low  and  high  confining  pressures. 

Figure  2.4  labels  four  curve  types,  and  Table  2.2  describes  their 
characteristics  and  the  typical  testing  conditions  leading  to  each  type.  In  way  of 
summary.  Type  C  curves  predominate  at  low  confining  pressures  (<T(.  =  0.1  MPa)  and 

are  characterized  by  post-yield  strain  hardening  to  give  a  peak  strength  (Qp)  at 
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moderate  strains,  followed  by  significant  strain  softening  and  significant  dilation. 

For  loose  sand  sheared  at  higher  strain  rates  and/or  lower  temperatures,  this 
behavior  switches  to  Type  A  curves,  i.e.,Qp  =Quy  followed  by  significant  strain 
softening  and  volumetric  expansion.  In  contrast.  Type  D  curves  predominate  at 
high  confining  pressures  10  MPa)  and  are  characterized  by  significant  post-yield 

strain  hardening  to  give  a  peak  strength  at  large  strains,  all  with  minimal  dilation. 
For  loose  sand  sheared  at  higher  strain  rates  and/ or  lower  temperatures,  this 
behavior  switches  to  Type  B  curves,  i.e.,  initial  strain  softening  followed  by  strain 

hardening  to  produce  Qp  «Quy/  again  with  minimal  dilation.  * 

6.3  Small  Strain  Behavior  of  Frozen  Sand 
6.3.1  Young's  Modulus 

Figure  3.1  demonstrates  that  the  Young's  modulus  (E)  is  nearly 
independent  of  all  testing  variables  for  conventional  tests  on  frozen  MFS.  There  is  a 

slight  (10-15%)  decrease  in  E  withOc  increasing  from  0:1  to  10  MPa  (considered  a 
second-order  effect),  and  a  slight  increase  (15%)  withDj  increasing  from  20  to  100%. 
The  measured  modulus  equals  26.5  ±  4.0  GPa  over  all  the  testing  variables.  This 
value  agrees  well  with  reliable  measurements  (i.e.,  using  on-specimen  extensometer 
or  wave  propagation)  on  three  other  frozen  sands  and  from  the  consolidate-freeze 
tests  on  MFS. 

The  mineralogy  of  the  sand  grains,  the  sand  density,  and  the  ice/ silicate 
bond  strength  are  considered  to  be  the  most  significant  factors  that  contribute  to  the 
initial  stiffness  properties  of  frozen  MFS  in  the  ranges  of  temperature  and  strain  rate 
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investigated  in  this  program.  A  composite  materials  model  developed  by  Counto 
(1964)  for  concrete  has  been  extended  to  frozen  sands.  The  excellent  agreement 
between  predicted  and  measured  data  in  Figure  3.1  indicates  that  composite  action 
between  the  pore  ice  and  silicate  particles  is  the  mechanism  responsible  for  a 
Young's  modulus  that  is  several  times  greater  than  that  of  polycrystalline  ice  and  up 
to  two  orders  of  magnitude  greater  than  that  of  the  unfrozen  sand  skeleton. 

6.3.2  Yield  Offset  Stress 

The  yield  offset  stress  (Qyo)  is  presented  in  order  to  bridge  the  gap 

between  the  initial  modulus  (E)  and  the  upper  yield  stress  (Quy)-  The  results  from 
conventional  tests  on  frozen  MFS  in  Figure  3.2  show  the  following  trends:  relative 
density  has  no  effect;  there  is  a  modest  decrease  in  Qyo  with  increasing  confining 
pressure;  and  increases  in  strain  rate  and  decreases  in  temperature  both  cause  a 
significant  increase  in  Qyo,  but  to  a  lesser  degree  than  for  the  upper  yield  stress  (see 
Figure  3.5  for  comparison).  This  indicates  that  the  behavior  of  frozen  sand  in  the 
small  strain  region  becomes  more  strain  rate  and  temperature  sensitive  as  the  level 
of  strain  increases.  In  addition,  consolidation  to  very  high  stresses  prior  to  freezing 

did  not  change  Qyo  (see  Sheet  D8). 

6.3.3  Upper  Yield  Stress 

The  use  of  on-specimen  axial  strain  measurements  enabled  reliable 
determination  of  the  upper  yield  stress  (Quy)  for  all  stress-strain  curves  from  the 

program  of  conventional  tests  on  frozen  MFS.  The  value  of  Quy  (which  occurs 
between  0.3  to  1.0%  axial  strain)  does  not  vary  with  relative  density  and  is  affected 
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only  slightly  by  the  magnitude  of  the  confining  pressure,  as  illustrated  in  Figure  3.3. 
This  indicates  that  the  frictional  characteristics  of  the  sand  skeleton  play  a  minor 
role,  since  density  and  confining  (consolidation)  pressure  have  a  major  influence  on 
the  small  strain  behavior  of  unfrozen  MFS  (e.g..  Figures  5.2  and  5.6).  However,  this 
important  conclusion  needs  further  study  because  of  the  following  observations: 

•  The  consolidation-freeze  tests  on  MFS  show  that  consolidation  to  high 
stresses  (o'c  ^  2  MPa)  prior  to  freezing  causes  a  stiffer  small-strain  response 
and  the  disappearance  of  a  distinct  upper  yield  point  (Figures  5.11  and 
5.12).  This  suggests  that  the  frictional  resistance  of  the  sand  skeleton  can 
become  significant  within  the  small  strain  region  for  dense  sand  having 
very  high  pre-freezing  consolidation  stresses. 

•  Some  prior  data  from  conventional  test  programs  show  a  very  large 

pressure  sensitivity  inQuy  (i.e..  Figure  3.6).  Whether  this  reflects  a  vastly 
different  behavior,  compared  to  frozen  MFS,  or  shortcomings  in  the 
experimental  procedures  (e.g.,  S-shaped  initial  stress-strain  curves) 
remains  unknown.  In  addition,  preparation  of  specimens  by  compaction 
rather  than  pluviation  alters  the  measured  stress-strain  response  (e.g.. 
Figure  5.10  and  Baker  and  Konrad,  1985).  This  fact  further  complicates 

direct  comparison  of  Qyy  data  for  different  frozen  sands. 

Figures  3.4  and  3.5  show  that  the  upper  yield  stress  of  frozen  MFS  is 
very  dependent  on  the  strain  rate  and  the  temperature.  Figure  3.7  compares  Qyy  of 
frozen  MFS  with  the  peak  strength  of  polycrystalline  (granular)  ice  as  a  function  of 
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strain  rate.  The  power  law  coefficients  (n)  for  the  ice  range  from  4.3  to  5.4  and 
compare  well  with  those  for  frozen  MFS  (n  =  5.6  ±  1.0).  This  comparison,  typical 

stress-strain  behavior  for  bulk  ice  (e.g..  Figure  4.22),  and  the  independence  of  Qyy  on 

sand  density  and  confining  pressure,  suggest  that  the  principal  physical  mechanisms 
occurring  within  the  upper  yield  region  are  probably  similar  to  those  controlling  the 
behavior  of  granular  ice  at  its  peak  strength.  Chamberlain  et  al.  (1972),  Bragg  and 
Andersland  (1980),  Parameswaran  and  Jones  (1981),  and  Sayles  (1989),  among  others, 
have  also  concluded  that  the  upper  yield  behavior  of  frozen  sand  is  dominated  by 
the  ice  matrix. 

The  authors  hypothesize  that  the  upper  yield  stress  of  frozen  MFS  as 
measured  in  conventional  tests  is  caused  almost  entirely  by  the  strength  of  the  ice 
matrix.  However,  Figure  3.7  shows  that  the  peak  strength  of  ice  (Qj)  is  about  one- 

half  of  Quy  for  testing  at  T  =  -10°C.  Moreover,  analyses  in  Section  5.5  suggest  that 
Qi/Quy  decreases  substantially  at  lower  temperatures  (e.g..  Sheet  E9).  Consequently, 

if  Quy  is  controlled  by  the  cohesive  strength  of  the  ice  matrix,  then  the  presence  of 
the  sand  grains  must  cause  a  very  significant  strengthening  of  the  ice  matrix, 
compared  to  bulk  polycrystalline  ice.  The  authors  agree  with  Ting  et  al.  (1983)  who 
attribute  this  enhanced  strength  of  the  intact  ice  matrix  to  a  difference  in  the  ice 
structure  (e.g.,  an  order  of  magnitude  or  more  decrease  in  the  grain  size  of  ice  in 
frozen  sand),  an  increase  in  strain  rate,  and  changes  in  the  state  of  stress  and  the 
added  deformational  restraints  due  to  the  presence  of  the  sand  grains.  However, 
these  potential  strengthening  mechanisms  are  highly  speculative  and  difficult  to 
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quantify.  In  addition,  the  validity  of  the  assumed  Qj  relationships  for  bulk  ice  as  a 
function  of  confinement,  strain  rate  and  temperature  requires  further  verification. 
More  research  is  also  needed  to  determine  possible  trends  inQ^y  =  fia^,  e  and  T),  in 
terms  of  sand  characteristics  such  as  grain  composition,  shape,  surface  roughness 
and  size  distribution. 

6.4  LARGE  STRAIN  BEHAVIOR  OF  FROZEN  SAND 
6.4.1  Volumetric  Behavior 

The  current  program  is  one  of  the  very  few  to  measure  volumetric 
strains  (Ey)  and  only  Shibatad  al.  (1985)  report  data  at  varying  confining  pressures, 

strain  rates  and  temperatures  (but  for  dense  sand  only).  Both  programs  show 
essentially  zero  volumetric  strain  below  the  upper  yield  stress  and  then  varying 
amounts  of  volumetric  expansion  (dilation)  that  continues  to  increase  with  further 
straining.  Figures  4.1  to  4.3  illustrate  typical  volumetric  behavior  for  frozen  MFS. 

Since  an  ice-saturated  frozen  sand  continuum  subjected  to  increasing 
axial  stress  can  undergo  only  volumetric  compression  (for  Poisson's  ratio  less  than 
0.5),  any  volumetric  expansion  (dilation)  must  reflect  fracturing  (cracking)  of  the  ice 
matrix  and  loss  of  bonding  between  the  ice  matrix  and  the  sand  particles.  In  other 
words,  the  onset  of  dilation  suggests  a  decrease  in  the  cohesive  strength  of  the 
frozen  sand  system  caused  by  damage  of  the  ice  matrix. 

Two  parameters  are  used  to  quantify  the  dilation  behavior  of  frozen 

MFS  at  large  strains:  the  maximum  rate  of  dilation  (MRD)  =  (de^/ dea)inax 
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volumetric  expansion  at  20%  axial  strain,  ev20-  illustrated  in  Figures  4.14  and 


4.15,  both  parameters  follow  similar  trends: 

•  increasing  relative  density  causes  an  increasing  dilation  at  low 

confinement  =  0.1  MPa),  with  minimal-to-no  influence  of  strain  rate 
and  temperature.  It  is  interesting  that  frozen  loose  MFS  undergoes  about 
twice  as  much  dilation  at  20%  axial  strain  as  unfrozen  MFS. 

•  increasing  confining  pressure  causes  a  very  significant  decrease  in  the  rate 
and  amount  of  dilation. 

•  at  high  confinement  (Oc  =  10  MPa)/  sand  density  plays  a  very  minor  role, 
whereas  higher  strain  rates  or  lower  temperatures  may  cause  a  slight 
increase  in  dilation. 

These  trends  are  used  to  help  explain  the  post-upper  yield  stress-strain  behavior  of 
frozen  sand. 

6.4.2  Post-Upper  Yield  Behavior:  Curve  Types  and  Conceptual  Model 

As  previously  emphasized,  the  post-upper  yield  behavior  of  frozen 
sand  is  highly  dependent  on  sand  density  and  confinement,  in  addition  to  strain 
rate  and  temperature.  The  four  curve  types  depicted  in  Figure  2.4  and  further 

described  in  Table  2.2  are  used  to  illustrate  this  dependency  onD^,  ^  and  T. 
Section  4.6  presents  a  conceptual  model  to  explain  the  post-upper  yield  strain 
hardening /strain  softening  and  volumetric  behavior  associated  with  these  four 
curve  types  in  terms  of  weakening  mechanisms  associated  with  damage  to  the  ice 
matrix  due  to  cracking  that  is  reflected  by  dilation  of  the  frozen  sand,  and 
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strengthening  mechanisms  associated  with  the  frictional  resistance  of  the  frozen 
sand  system. 

Weakening  and  Strengthening  Mechanisms 

Figure  4.22  illustrates  typical  stress-strain  curves  for  polycrystalline  ice 
for  the  range  of  testing  conditions  used  for  frozen  MFS.  According  to  Mellor  and 
Cole  (1982)  and  Murrell  et  al  (1989),  the  onset  of  internal  cracking  begins  at  the  first 
yield  point.  The  continued  formation  and  growth  of  aacks  result  in  the  peak 
strength  (Qj)  occurring  at  about  0.5  to  1%  axial  strain,  followed  by  strain  softening 

that  increases  at  higher  strain  rates  and  lower  temperatures,  but  decreases  with 
added  confinement.  However,  assuming  that  the  ice  matrix  in  frozen  sand  behaves 
similarly  to  bulk  ice,  one  would  always  expect  a  significant  loss  in  the  strength  of  ice 

at  axial  strains  greater  than  about  0.5  to  1%  (i.e.,  after  Qj  for  ice  andQ^y  for  frozen 
sand).  This  represents  the  weakening  mechanism. 

On  the  other  hand,  continued  straining  of  frozen  sand  beyond  Quy  also 
will  produce  a  strengthening  mechanism  due  to  the  frictional  resistance  of  the  ice- 
sand  composite  that  progressively  behaves  more  like  a  particulate  system,  due  to 
increasing  damage  of  the  ice  matrix  and  increasing  interaction  between  the  particles 
of  the  sand  skeleton.  Section  4.6.3  contains  a  detailed  discussion  of  the  potential 
strengtherung  due  to  sand  skeleton  friction  if  this  frictional  resistance  exhibited  the 
same  behavior  as  observed  for  unfrozen  sand  (which  is  not  true,  but  still  may 
provide  general  guidance).  In  addition,  one  can  envision  strengthening  due  to 
interference  between  the  sand  particles  and  the  damaged  ice  matrix,  and  perhaps 
some  frictional  resistance  of  the  fractured  ice  matrix.  In  any  case,  these 
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strengthening  mechanisms  should  increase  with  increasing  relative  density  and 
increasing  confining  pressure. 

Based  on  the  observed  volumetric  behavior  of  frozen  MFS,  one  could 
expect  the  ice-sand  system  to  behave  in  an  "undrained"  manner  at  low-to-moderate 

strains  (say  up  to  £3  *  5%),  since  the  amount  of  dilation  is  very  small,  even  for  tests 

at  low  confinement  (Oj,  =  0.1  MPa).  For  a  truly  constant  volume  condition,  the 
effective  stresses  producing  frictional  resistance  of  the  sand  skeleton  are  generated 
within  the  ice-sand  system.  On  the  other  hand,  as  dilatancy  becomes  more 
pronounced  at  larger  strains,  the  ice-sand  system  could  be  expected  to  approach  a 
"drained"  condition.  For  a  truly  drained  shear  condition,  the  effective  stresses 
causing  frictional  resistance  are  generated  by  the  applied  stresses,  i.e.,  by  the  axial 
stress  and  the  confining  pressure. 

Behavior  at  Low  Confinement  (<Tc  =  0.1  MPa) 

Type  C  curves  represent  the  most  prevalent  type  of  behavior.  After 
reaching  the  upper  yield  stress,  the  frozen  sand  exhibits  strain  hardening,  since  the 
strengthening  mechanisms  are  stronger  than  the  weakening  mechanisms.  The 

degree  of  strengthening  and  the  resultant  peak  strength  (Qp)  increase  with 
increasing  relative  density  (Dj),  as  illustrated  in  Figures  4.6  and  4.17,  probably  similar 

to  imdrained  shear.  With  continued  straining,  dilatancy  reduces  the  internal 
frictional  resistance,  and  the  weakening  mechanisms  now  predominate  to  cause 
strain  softening.  Moreover,  the  normalized  rate  of  strain  softening  is  uniquely 
related  to  the  amount  of  dilation,  as  per  Figure  4.21  and  Section  6.4.4. 
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Faster  strain  strains  and  lower  temperatures  produce  Type  A  curves, 
e.g.,  test  number  36  in  Figure  2.2a  and  test  number  97  in  Figure  4.1b,  especially  at  the 
lower  sand  densities.  The  weakening  mechanisms  dominate,  due  to  both  a  more 
rapid  rate  of  ice  cracking  and  the  relatively  low  post-yield  frictional  resistance.  It  is 
important  to  note  that  the  normalized  rate  of  strain  softening  versus  volumetric 
expansion  relationship  for  Type  A  curves  is  the  same  as  for  Type  C  curves  (Fig.  4.21). 
However,  Type  A  curves  exhibit  a  more  rapid  initial  strength  loss  immediately  after 

yielding  at  Qp  =  Q^y  than  occurs  after  the  peak  strength  for  Type  C  curves.  But  both 
curve  types  show  an  almost  linear  decrease  in  strength  and  increase  in  dilation  at 
the  larger  strain  levels  (e.g..  Fig.  4.1). 

Behavior  at  High  Confinement  (Cc  =  10  MPa  ) 

Type  D  curves  represent  the  most  prevalent  type  of  behavior.  After 
reaching  the  upper  yield  stress,  the  frozen  sand  exhibits  continued  strain  hardening 
until  reaching  the  peak  strength  at  large  strains.  The  amoimt  of  strength  increase 
with  increasing  relative  density  is  much  larger  than  for  tests  at  low  confinement 
(Fig.  4.6),  since  high  confinement  increases  the  frictional  resistance  (more 
strengthening)  and  reduces  the  degree  of  ice  cracking  (less  weakening). 

Faster  strain  rates  and  lower  temperatures  produce  Type  B  curves,  e.g., 
test  number  129  in  Figure  4.3b,  especially  at  the  lower  sand  densities.  The 
weakening  mechanisms  dominate  immediately  after  yielding  for  the  same  basic 
reasons  as  for  Type  A  curves.  But  thereafter,  the  increased  cor\finement  produces 
sufficient  frictional  resistance  to  offset  damage  of  the  ice  matrix. 
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6.4.3  Peak  Strength:  Background 

Table  B  in  Appendix  B  summarizes  the  scope  of  17  prior  experimental 
programs  that  have  investigated  the  peak  strength  of  saturated  (or  nearly  so)  frozen 
sands.  Of  these  programs,  only  one  measured  on-specimen  axial  strains  (and,  hence, 
could  identify  the  upper  yield  stress),  only  one  included  a  detailed  study  of  relative 
density  (D^)  (although  two  included  ice-sand  mixtures),  and  only  one  included 

variations  in  confining  pressure  (o^),  strain  rate  (^)  and  temperature  (T).  In 
addition,  many  programs  used  compacted  rather  than  pluviated  specimens,  which 
leads  to  non-uniform  densities  and  more  scattered  results.  Hence,  the  scope  of  the 
testing  program  shown  in  Table  2.1  for  frozen  MFS  is  unique. 

Figure  4.4  compares  the  peak  strength  (Qp)  of  frozen  MFS  with  data  on 

other  frozen  sands  as  a  function  of  strain  rate  (e)  at  temperatures  ranging  from 
T  =  -10°  to  -30°C.  The  comparison  is  restricted  to  dense  sands,  since  all  prior  studies 
of  e  and  T  used  only  dense  specimens.  For  testing  at  low  confinement,  Figure  4.4a 
shows  that  the  location  of  MFS  and  its  values  of  the  power  law  coefficient  (n)  are 
representative  of  frozen  sand  data,  except  for  the  strain  rate  insensitivity  observed  in 
Program  A  above  ^  =  10  sec  for  a  30-to-140-mesh  silica  sand  (sunilar  trends  also 

occurred  at  T  =  -6  and  -15°C).  But  also  note  the  two-to-three-fold  variations  in  Qp  at 
the  same  strain  rate  for  sands  tested  at  T  =  -10°C.  This  large  difference  presumably 
reflects  differences  in  relative  density  (and  maybe  method  of  preparation),  degree  of 
ice  saturation  (all  used  pure  water)  and  the  grain  size  distribution.  For  testing  at 
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high  confinement.  Figure  4.4b  shows  close  agreement  between  MFS  and  the 
uniform  fine  sand  used  in  Program  1. 

Figure  4.12  compares  the  upper  yield  stress  (Quy)  and  peak  strength  (Qp) 

of  frozen  MFS  as  a  function  of  strain  rate  and  temperature  at  two  extreme 
conditions:  low  density  and  low  confinement  (hence,  minimal  effect  of  sand 
skeleton  frictional  resistance);  and  high  density  and  high  confinement  (hence, 
maximum  effect  of  sand  skeleton  frictional  resistance).  This  comparison  illustrates 

the  extreme  importance  of  comparing  these  two  stresses,  since  Quy  is  thought  to  be 

dominated  by  the  strength  of  the  ice  matrix.  Therefore,  Quy  is  an  important  baseline 

strength.  The  comparison  also  demonstrates  that  inspection  of  Qp  data  alone  (as 

done  in  Figure  4.4)  provides  little  insight  regarding  the  physical  mechanisms 
controlling  the  peak  strength  of  frozen  sands.  Consequently,  evaluation  of  the 
variables  that  affect  the  peak  strength  of  frozen  MFS  is  made  in  terms  of  the  strength 
increase  relative  to  the  upper  yield  stress.  In  other  words,  the  results  are  analyzed  by 

using  AQ  =  Qp  -  Quy  or  Qp/Quy,  in  order  to  quantify  the  net  effect  of  the  post-upper 
yield  strengthening  and  weakening  mechanisms.  Note  that  these  analyses  are 
restricted  to  curve  Types  C  and  D  (Figure  2.4),  since  Qp  =  Quy  for  curve  Types  A  and 
B.  The  conclusions  from  these  analyses  are  summarized  in  Section  6.4.5;  but  for 
perspective,  that  summary  will  be  introduced  with  a  description  of  a  methodology 
for  predicting  the  maximum  strengthening  that  can  occur  due  to  the  frictional 
resistance  of  the  sand  skeleton  based  on  Ladanyi  s  dilatancy-hardening  model. 
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6.4.4  Peak  Strength:  Dilatancy-Hardening  Model 

Assumptions 

Section  5.2.1  describes  the  model  presented  by  Ladanyi  (1985)  and 
Ladanyi  and  Morel  (1990)  for  predicting  the  peak  strength  of  frozen  dense  sand. 

Basically,  the  model  assumes  that  the  predicted  strength  (Qps)  is  given  by 

Qfs  =  Qi  +  Qfs  . 

Qj  is  the  strength  of  the  ice  matrix;  it  is  equated  to  the  peak  strength  measured  in 

triaxial  compression  tests  run  on  bulk  polycrystalline  ice.  Qfg  is  the  frictional 

resistance  of  the  sand  skeleton;  it  is  equated  to  the  peak  strength  computed  for 
imdrained  shear  in  triaxial  compression.  This  frictional  resistance  is  basically 
computed  as 

Qfe  =  a’3f(Rf-l)  . [6.2] 

where  <y'3f  is  the  minor  principal  stress  at  failure  for  undrained  shear  of  the 
unfrozen  sand  skeleton,  and  Rf  is  the  obliquity  of  the  principal  stresses  at  failure, 
(a’i/a’3)  =  (1  +  sin  <1)')/(1  -  sin 

The  prediction  of  in  the  dilatancy-hardening  model  depends  upon 

whether  or  not  cavitation  of  the  pore  ice  occurs  during  constant  volume  shear  of 
the  frozen  sand.  For  combinations  of  sand  density  and  confining  pressure  that  do 

not  cause  cavitation  (i.e.,  tests  with  low  density  or  high  confinement),  a’3f  is  given 
by  the  critical  confining  pressure  (a'c),  which  is  uniquely  related  to  the  preshear 
density  of  the  (frozen)  sand  specimen.  For  combinations  of  sand  density  and 
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confining  pressure  that  do  cause  cavitation  of  the  pore  ice,  (i.e.,  high  density  and  low 
confinement),  a'sf  equals  plus  Tj,  where  is  the  applied  confining  pressure,  and 

Tj  is  the  tensile  strength  of  the  pore  ice  (equated  to  the  tensile  strength  of  bulk 
polycrystalline  ice).  Consequently,  the  predicted  peak  strength  (Qps)  foJ^  the  non¬ 
cavitation  case  is  give  by 

QFS  =  Qi  +  [Qfs  =  <^’cr(Rf-l)]  . [6.3a] 

and  for  the  cavitation  case,  it  is  given  by 

Qfs  =  Qi  +  [Qfs  =  +  Ti)  (Rf  - 1)] . [6.3b] 

Model  Parameters  for  Frozen  MFS 
Section  5.5.1  describes  the  approach  used  to  select  parameters  for  the 

pore  ice.  Values  for  the  peak  strength  of  ice  (Qj)  were  obtained  by  using  a  combined 
Arrhenius-Glen  (1955)  power  law  equation  to  fit  the  Jones  (1982)  data  at  T  =  -11.5°C 
for  Cc  =  0  arid  10  MPa  and  then  extrapolated  to  other  temperatures,  assuming  a 
constant  activation  energy  (Sheets  E4  and  E6).  This  approach  was  adopted  due  to  the 
lack  of  experimental  data  at  lower  temperatures  for  the  strain  rates  of  interest;  it 

involves  significant  uncertainty  in  the  Qj  values.  The  tensile  strength  of  pore  ice 
was  taken  as  Tj  =  2.15  MPa,  based  on  the  data  summarized  in  Sheets  El  to  E3. 

The  extensive  program  of  consolidated-undrained  (CU)  and 
consolidated-drained  (CD)  triaxial  compression  tests  on  unfrozen  MFS  (Section  5.3) 
provided  two  types  of  vital  information:  1)  how  changes  in  sand  density  and 
confining  (consolidation)  stress  cause  vastly  different  stress-strain  curves  and 
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effective  stress  paths  (e.g..  Figures  5.2  to  5.6);  and  2)  reliable  predictions  of  the 
effective  stress  friction  angle  (())',  Figure  5.7)  and  of  the  variation  in  critical  confining 

pressure  (o’er)  ^^ith  sand  density.  Determination  of  a’^r  involved  the  following 
steps: 

•  definition  of  the  steady  state  line  (SSL),  which  represents  a  state  of 
continued  shearing  at  constant  void  ratio  (density),  effective  stress  and 
shear  stress  (Sheet  Cll).  This  step  required  very  extensive  experimental 
data. 

•  use  of  the  state  parameter  \y  proposed  by  Been  and  Jefferies  (1985)  to  predict 
how  changes  in  density  and  confinement  affect  Skempton  s  pore  pressure 
parameter  A  at  failure  [Af  =  Au  -  A03)/  ( Aoj  -  A<73)],  i.e.,  the  value  of 

\|r  =  -0.061  that  produced  Af  =  0  (Figure  5.9). 

•  definition  of  the  critical  confining  pressure  line  (Figure  5.8). 

Sheet  E12  shows  how  relative  density  (D^)  affects  the  critical  confining 

pressure  (o'er)  predicted  frictional  resistance  of  the  sand  skeleton  (Qfs)  for 

frozen  MFS  having  Tj  =  2.15  MPa  and  tested  at  Oc  =  0.1  MPa  confinement.  Several 

important  conclusions  can  be  drawn  from  these  model  predictions: 

•  The  non-cavitation  case  applies  for  all  relative  densities  less  than 

Dr  =  90%,  i.e.,  for  <  Oc  +  Tj  =  2.25  MPa.  Within  this  region,  Qfs  varies 
from  0.7  MPa  at  D^  =  20%  to  Qfs  =  5.8  MPa  at  D^  =  90%.  Moreover,  increases 
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in  confinement  do  not  affect  the  predicted  undrained  frictional  resistance 
of  the  sand  skeleton. 

•  The  cavitation  case  occurs  only  for  relative  densities  between 

Dj  =  90  and  100%.  Within  this  region,  the  predicted  frictional  resistance 

remains  constant  at  Qfg  =  6.0  MPa  independent  of  D,  for  =  0.1  MPa. 
However,  increases  in  will  cause  a  large  increase  in  Qf^  (corresponding 
to  (|)  =  4)’  =  34.9°),  as  long  as  Cf.  +  Tj  remains  less  than  a'cr 

•  Uncertainty  in  the  tensile  strength  of  ice  (Tj)  will  not  affect  the  predicted 

Qfs  of  frozen  MFS,  except  for  dense  sand  at  low  confinement. 

Comparison  of  Predicted  and  Measured  Strengths 
Section  5.5.2  presents  a  detailed  comparison  of  predicted  (Qps)  and 
measured  (Qp)  strengths  for  conventional  frozen  tests  on  dense  (D,  =  95%)  and  loose 

(Dj  =  35%)  MFS  at  low  (Oc  =  0.1  MPa)  and  high  (Oc  =  10  MPa)  confinement,  e.g.. 

Figures  5.15  to  5.18.  The  overall  agreement  is  very  poor,  except  for  relatively  slow 
shearing  at  T  =  -10°C  of  dense  sand  at  low  confinement.  Perhaps  by  coincidence, 
Ladanyi  and  Morel  (1990)  obtained  excellent  agreement  for  a  dense,  20-30  Ottawa 
sand  sheared  at  =  0  to  0.3  MPa,  ^  =  2.7  x  10“*/ sec  and  T  =  -5°C.  (Their  prediction 
used  the  cavitation  equations;  see  Section  5.2.2). 

Sheet  E8  quantifies  the  model  error  by  plotting  Qps/ Qp  versus 

temperature  for  dense  and  loose  frozen  MFS  at  =  0.1  and  10  MPa.  For  dense  sand. 
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the  Qp  values  range  from  10  to  25  MPa  (compared  to  Qfg «  6  ±  1  MPa),  and  Qps/Qp 
decreases  from  roughly  0.7  ±  0.1  at  T  =  -10°C  down  to  about  0.5  at  T  =  -25°C.  For  loose 
sand,  the  Qp  values  range  from  7  to  30  MPa  (compared  to  Qjg  of  less  than  1.5  MPa), 
and  Qps/Qp  decreases  from  roughly  0.6  ±  0.1  at  T  =  -10°C,  down  to  0.45  ±  0.1  at 
T  =  -25°C.  The  fact  that  the  predicted  Qfs  represents  a  small  fraction  of  the  measured 
Qp  (at  least  for  loose  sand)  suggests  that  the  main  problem  with  the  model 
predictions  lies  with  the  predicted  values  for  the  strength  of  the  ice  matrix  (Qj). 

Section  6.3.3  hypothesized  that  the  upper  yield  stress  (Quy)  may  represent  the 
operational  or  effective  strength  of  the  ice  matrix  prior  to  the  onset  of  significant 
cracking.  Sheet  E9  plots  the  predicted  Q,  divided  by  the  measured  Quy  versus 
temperature  for  testing  at  Oc  =  0  and  10  MPa.  Since  the  magnitudes  and  trends  in 
Qps/Qp  and  Qj/Quy  are  amazingly  similar,  the  authors  conclude  that  Qj  should  be 
replaced  by  Quy  when  applying  the  dilatancy-hardening  model  to  frozen  MFS.  The 
revised  predictions  are  denoted  by  Qps*  =  Quy  Qfs/  where  Qfs  is  again  given  by 
Eq.  6.3.  Note  that  the  predicted  Qps*  minus  the  measured  Qp  should  equal  Qfs  for 
perfect  predictions.  The  next  section  includes  such  comparisons. 

6.4.5  Peak  Strength:  Effect  of  Testing  Variables 

As  noted  earlier,  analysis  of  the  influence  of  testing  variables  on  the 

peak  strength  is  made  primarily  in  terms  of  AQ  =  Qp  -  Quy  ^or  the  combination  of  D,, 
Oc,  ^  and  T  producing  curve  Types  C  and  D  that  strain  harden  after  reaching  the 
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upper  yield  stress  (Quy).  Note:  dQ  will  be  used  to  denote  changes  in  both  dAQ  and 
dQp  as  a  function  of  changes  in  D,,  Oc/ 

Effect  of  Relative  Density  (D^) 

Figure  4.6  shows  a  linear  increase  in  AQ  with  increasing  D,  for  frozen 

MFS.  The  rate  of  increase  (dQ/dD,)  is  roughly  independent  of  e  and  T,  but  it  is 
highly  dependent  on  confinement,  as  further  illustrated  in  Figure  4.7.  At 
Oc  =  0.1  MPa,  dQ/dDr  =  5.6  ±  1.0  MPa  and  almost  doubles  to  11.9  ±  0.6  MPa  at 

Oc  =  10  MPa.  The  only  prior  program  with  a  detailed  study  of  relative  density 
(program  R  in  Table  B)  showed  a  much  higher  dQp/ dD^  =  13  MPa  from  unconfined 

compression  tests  on  a  pluviated  16-100  Ottawa  sand  at  Dj>  65%.  Both  the  original 
and  revised  versions  of  the  dilatancy-hardening  model  predict  a  non-linear  increase 
in  AQ  =  Qfs  with  D^,  i.e..  Sheet  E12  predicts  a  linear  -  log  Qfg  relationship.  In  any 
case,  further  research  is  needed  to  determine  how  relative  density  affects  post-upper 
yield  behavior  for  further  comparison  with  the  frozen  MFS  data. 

Figure  4.11  shows  results  of  analyses  to  obtain  values  of  D,  producing 

Qp  =  Qyy  (i.e.,  AQ  =  0  for  Type  A  curves)  for  =  0.1  MPa  tests  on  frozen  MFS  as  a 
function  of  ^  and  T.  Although  the  single  data  point  from  program  R  agrees  with 
the  frozen  MFS  data,  added  research  is  needed  to  predict  conditions  leading  to 
failure  at  the  upper  yield  stress. 
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Ejfect  of  Confining  Pressure  (afi 


Figure  4.8  shows  the  effect  of  confinement  for  frozen  MFS  tested  at  the 
moderate  strain  rate  and  T  =  -10°C.  The  pressure  sensitivity  (in  terms  of  dQ/ dOc  and 

the  corresponding  total  stress  friction  angle,  0)  is  much  higher  for  dense  than  for 
loose  frozen  MFS  and  tends  to  decrease  at  very  high  confining  pressures  (especially 
for  loose  sand).  Mean  values  of  dQ/ d  a,,  (from  data  at  =  0.1  and  10  MPa)  for  slow 
and  moderate  shearing  range  from  about  0.2  to  0.7  (corresponding  (j)  =  5  to  15°)  and 

vary  in  a  complex  fashion  with  D^,  e  and  T,  as  summarized  in  Figure  4.9. 

Several  prior  programs  have  studied  the  effect  of  confinement  on  the 
peak  strength  of  medium-dense  to  dense  frozen  sands.  Figure  4.10  summarizes  data 

for  tests  conducted  at  T  =  -11  ±  1°C.  The  results  from  the  four  test  series  on  medium 
to  fine  sands  (programs  G,  I  and  L;  also  see  Table  B)  give  dQp/ dOc  *  0.56  to  0.9 
(<j>  =  12.5°  to  18°)  and  are  consistent  with  data  for  dense  MFS.  However,  tests  on  the 
coarse  20-30  Ottawa  sand  used  for  program  J  produced  a  much  higher  pressure 
sensitivity  (([)  =  28°).  In  addition,  Shibata  et  al  (1985)  did  not  observe  a  decrease  in 

pressure  sensitivity  at  faster  strain  rates  and  lower  temperatures,  as  observed  for 
frozen  MFS.  Hence,  the  pressure  sensitivity  of  frozen  sands  does  not  appear  to 
follow  a  consistent  trend. 

As  shown  in  Figure  5.14,  consolidation  of  dense  MFS  to  o'c  =  0.1  to  10 

MPa  prior  to  freezing  (i.e.,  results  from  the  consolidate-freeze  test  program 
presented  in  Section  5.4)  did  not  produce  peak  strengths  very  different  from  those 
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measured  in  conventional  frozen  tests.  However,  consolidation  to  high  stresses 
(Oc  >  2  MPa)  did  cause  a  very  significant  increase  in  the  deviator  stress  (Q)  at  axial 

strain  levels  from  about  =  0.5  to  5%,  as  illustrated  in  Figures  5.11  and  5.12.  This 

pronounced  stiffening  is  consistent  with  the  behavior  observed  in  consolidated- 
undrained  triaxial  compression  tests  on  unfrozen  MFS,  as  per  Figure  5.2.  The  stress- 
strain  curve  in  Figure  5.2  for  test  C-29  also  suggests  that  anisotropic  consolidation 
(i.e.,  higher  vertical  than  horizontal  stresses)  will  greatly  enhance  this  stiffening  at 
small  strains. 

Figure  5.14  also  compares  dilatancy-hardening  model  predicted 
strengths  (Qfs)  with  the  measured  peak  strengths  (Qp)  for  dense  MFS  at  the 
moderate  strain  rate  and  T  =  -10°C.  The  model  vastly  underestimates  the  measured 
pressure  sensitivity,  because  it  predicts  no  increase  in  the  frictional  resistance  (Qfg) 
of  the  sand  skeleton  with  added  confinement  (the  modest  increase  in  Qps  occurs 
only  from  the  predicted  increase  in  Qj  with  confinement).  The  stress-strain  data  in 

Figure  5.2  also  show  little  change  in  the  peak  undrained  strength  of  unfrozen  sand. 
However,  increases  in  corxsolidation  stress  do  cause  a  much  stiffer  response  at 
smaller  strains.  Hence,  strain  compatibility  between  the  cohesive  strength  of  the  ice 
matrix  and  the  frictional  resistance  of  the  sand  skeleton  probably  plays  an  important 
role  in  the  behavior  of  the  ice-sand  composite  material. 

Effects  of  Strain  Rate  and  Temperature  (6  and  T) 

Figure  4.6  plots  the  post-upper  yield  strength  gain,  AQ  =  Qp  -  Quy,  versus 
relative  density  for  frozen  MFS  at  low  and  high  confinement  (Oj.  =  0.1  and  10  MPa). 


183 


In  general,  shearing  at  the  faster  strain  rates  and  lower  temperatures  tends  to  cause  a 
downward  shift  in  the  AQ  versus  Dj,  relationships.  Figure  4.13  provides  a  more 

detailed  picture  of  the  influence  of  e  and  T  by  plotting  AQ  for  loose  (D,  =  30%)  and 

dense  (D^  =  90%)  specimens  (interpolated  from  Figure  4.6)  versus  log  e  for  the 
different  temperatures.  Inspection  of  this  figure  indicates  the  following  complex 
trends,  using  the  slow  strain  rate  at  -10°C  as  the  reference: 

•  For  loose  sand  at  low  confinement,  the  minimal  strength  gain  rapidly 
vanishes  at  higher  strain  rates. 

•  For  loose  sand  at  high  confinement  and  for  dense  sand  at  low 
confinement,  the  moderate  strength  gain  disappears  more  rapidly  with 
faster  shearing  at  the  lower  temperatures. 

•  For  dense  sand  at  high  confinement,  the  large  strength  gain  is  temperature 
independent  at  the  slow  rate;  for  faster  shearing,  the  amount  of 
strengthening  reduces  at  the  lower  temperatures. 

Revised  Model  Predictions 

Sheet  ElO  plots  the  predicted  Ops*  =  Quy  +  Qfs  minus  the  measured  peak 
strength  (Qp)  versus  temperature  for  loose  (D^  =  35%)  and  dense  (D^  =  95%)  MFS  at 
low  and  high  confinement  =  0.1  and  10  MPa)  for  tests  having  curve  Types  C  and 

D.  This  comparison  generally  shows  that  the  predicted  Qps*  is  either  quite  accurate 
or  slightly  too  high  at  low  confinement  for  combinations  of  6  and  T  producing  the 
maximum  values  of  AQ  =  Qp  -  Quy.  However,  at  high  confinement  (again  for 
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conditions  producing  the  maximum  AQ),  the  predicted  Qps’^  is  too  low.  After 
considerations  of  strain  compatibility,  it  is  tentatively  concluded  that  the  measured 
strengths  are  larger  than  would  be  expected,  based  on  the  separate  behaviors  of  an 
unfrozen  sand  skeleton  and  an  enhanced  ice  matrix  strength  that  follows  trends 
observed  for  bulk  ice.  That  is,  the  strain  levels  at  the  peak  strength  of  frozen  MFS 
are  generally  too  small  to  mobilize  the  full  frictional  resistance  of  the  sand  skeleton 

(Qfs);  and  the  strains  are  sufficiently  large  to  cause  significant  weakening  of  the  ice 
matrix. 

6.4.6  Post-Peak  Strain  Softening  Behavior 

The  use  of  lubricated  end  platens  and  measurement  of  volumetric 
strains  have  made  possible  a  unique  evaluation  of  the  post-peak  strain  softening 
behavior  of  frozen  MFS.  For  tests  at  low  confinement  that  produce  curve  Types  A 
and  C  (i.e.,  post-peak  strain  softening),  there  appears  to  be  a  unique  relationship 
between  the  rate  of  strain  softening  (after  normalization  to  the  upper  yield  stress) 
and  the  amount  of  volumetric  expansion  (dilation).  Figure  4.21  presents  this 

relationship  in  terms  of  the  normalized  rate  of  strain  softening,  NRSS  =  [(Qp  - 
Q2o)/Quy]/(  ^20  ■  ^)'  versus  the  amount  of  dilation  at  20%  axial  strain,  ev2o-  The 
collective  data,  which  cover  a  wide  range  of  relative  densities  (D,  =  35  to  95%),  have 

relatively  little  scatter  about  the  mean  trend  and  do  not  show  a  consistent  deviation 
as  a  function  of  strain  rate  or  temperature.  However,  the  individual  stress 
components  of  the  NRSS  do  vary  sigiuficantly  with  the  testing  variables,  as 
summarized  below  (for  curve  Types  A  and  C): 
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•  Qp/Quy  increases  linearly  with  relative  density;  these  lines  move 
downward  and  become  flatter  with  increasing  strain  rate  and  decreasing 
temperature.  (See  Fig.  4.17  for  Type  C  curves.) 

•  Q2o/Quy  is  independent  of  relative  density,  but  also  has  a  smaller 
magnitude  with  increasing  strain  rate  and  decreasing  temperature.  (See 
Fig.  4.18  for  both  Type  A  and  C  curves.) 

Analysis  of  test  results  at  all  confining  pressures  (i.e.,  up  to  =  10  MPa) 

produced  the  results  presented  in  Figures  4.19  and  4.20.  For  all  curve  types  (i.e., 
including  Type  B  and  D  with  little  or  no  post-peak  strain  softening),  the  rate  of 
strain  softening  increases  with  the  maximum  rate  of  dilation,  and  the  stress  ratio 

t 

Qio/Qp  decreases  with  the  amount  of  dilation.  Although  both  plots  have  significant 

scatter,  they  provide  additional  evidence  that  the  degree  of  post-peak  strain 
softerung  or  strength  loss  is  strongly  related  to  the  rate  or  amount  of  dilation. 
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APPENDIX  A:  Data  on  Small-Strain  Behavior  of  Frozen  MFS 


Sheet  No.  Description 

Tables  from  Andersen  et  al  (1994)  giving  specimen  properties, 
testing  conditions  and  measured  small  strain  parameters 


Al&cA2 

A3 

A4 

A5 

A6 


Table  5.  Results  at  T  =  -10°C 

Table  6.  Results  at  T  =  -15°C 

Table  7.  Results  at  T  =  -20°C 

Table  8.  Results  at  T  =  -25°C 

Results  of  Linear  Regression  Analyses  on 
Upper  Yield  Stress  of  Frozen  MFS 
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TABLE  5.  Specimen  properties,  testing  conditions  and  measured  small  strain  parameters: 

Results  at  T  =  -10°C  (Page  1/2) 
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TABLE  5.  Specimen  properties,  testing  conditions  and  measured  small  strain  parameters 

Results  at  T  =  -10°C  (Page  2/2) 
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TABLE  6.  Specimen  properties,  testing  conditions  and  measured  small  strain  parameters: 
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TABLE  T.  Specimen  properties,  testing  conditions  and  measured  small  strain  parameters: 

Results  at  T  =  -20°C 


Q  rt  Q  o  O  O  O  O  O  O  21  21  O  O 

OOOOctfcdOOOcdOOp  P  C«  ed  O 

OJ)  OO  &i)  U-j  60  DO  W)  O*  DO  DO^  ^  Cli 

^^0'^vooooo^^o^^oc^vo^^o^c^o^ooc^^^^^ 

0^00000000000000000 

inoo-H*— *vn»n\o»nino\322  22S!552iSf:: 


csvo'noo^^r^^o  ^  oo^ 

0'S  1  00  \o  O ^ 2: 00  21  00  z VO  2:  00  fr> 


wa  C\  CO  ON  K  vn  CS  ^  Ov  <S  CN  oo  ^  oo  «o 

(5J  #«•••••••••  •  ♦•  •• 

t— '  tu  cu  to  ^  ^ ON  <s  o  rsi  ^  vn  ■<  CO  ^  00  CO -<  VO  o 

O  cs  CO  CO  CN  CS  CS  CO  <S  CN  CO  CN  21  CO  IZ  »-H  cs  Z  CS  CO 


O  5  ON  VO  VO  00  00  VO  VO  ^  00  <S  ON  VO  00  ^  CO  CO  CO  00 

c/a  S  00 •  •  VO  VO  NO  VO  •  CO  CO  VO  •  vo  vo  •  o  vo 

^  o^rtvocococococor^pN;5;»cococo^coco 

r*  X^COCO  COVOVOVOCO  COVO 

I  ed 

o  s 


ovo^cococovo25sO^S2S*^r^5SSSS 

^COONOCOCOCOCOCSpoOCSONCSCSCO  "^CS 
CO  <SCO  <S  ^  ^  CN 


o 

cd 

oc 

r*  _ 

<s 

9 

o 

Cd 

o 

1-1 

O' 

CO 

b  Cd  S 


csvovor^ONr^-^vo 
^  ^  ^  ^ 

*-H  ^  CS  ON 


^  vn  vo  o 
^  *-H  vo  ^ 

<s 


CSVONOvo^COCSVOCSCOONCO^^COOn*-«»«hON 

O  O  O  O  O  O  O  O  O  OONQOOOONOpON 

c>ic>j<sr'irQ(Scscscscs^<scscscs^cscst-H 

I  I  I  I 


cscscscscscscscscscs 

I  I  I  I  t  I  I  I  I  I 


OONOON^*-HCMVOr^^^^r;j^O;;:JCOC^CO 

o  o  o  ON  ooooooooooooogo 

•  •OOn  •  •  *0  •  *  *00  •  •  *0  *0 

OO’^  000’^000^’^000*^0»^ 


oouooooouooooouoooo 


o  o  a< 


CUOuOnO-iOuOxOuC^CUO^OuDmOLiOmOuD-iCuCUCU 

COSSC/OC/DC/DC/DC^C/DCOCOC/DCOC^COCOCOCOCOCO 


vor^or^csONONOr-vovocsvooocsvo’^r^oo 

^-lOOr^OOONONONOOOOOvOONOOONOONOsONQ 

OOnOnOnOnOnOnOONOnQnOnONOnOOnOnOnO 


'^cor^cNJOcovoo^^^vo'^cscor^cococN 

rfo^coooor^t^ovor^vnr^covovo-rfvc5c4 

OOONONCO^^t^COVOCOCOCOr^ONOOCOCOCOOO 


€d  o  t 

(X  9 

CO  CU 

O  JO  ”■ 
Vi  ed  ed 

p<*N  cd 

O  cd  ’►-* 

o 

C 


cs  CO  vo  CO  vo^^  csoN’^ovor^voovocsp 
^^r-ONco^fMO^ONOOcsr^'^oa\o^ 
NO^SVO’^VOVOVOVOVO’^VOVOVOVOVOVO'^VOVO 


o^(sj^u-^vovovor^oooN<sco'^vor^obONO 

t^f^t^l-^t^t^OOOOOOOOOOONONONONONONONO 


TABLE  8.  Specimen  properties,  testing  conditions  and  measured  small  strain  parameters: 

Results  at  T  =  -25“C 
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APPENDIX  B:  Programs  Reporting  Large-Strain  Behavior  of  Saturated 
Frozen  Sands  and  Data  on  Large-Strain  Behavior 


Sheet  No.  Description 


B1  to  B4  Table  B  Summary  of  Programs  Reporting  Large-Strain  Behavior  of 
Saturated  Frozen  Sands  (Modified  from  Andersen  1991) 

Tables  from  Andersen  (1991)  and  Swan  (1994)  Giving  Specimen 
Properties,  Testing  Conditions  and  Measured  Parameters 

B5  &  B6  Results  at  T  =  -10°C 

B7  Results  at  T  =  -15°C 

B8  Results  at  T  =  -20°C 

B9  Results  at  T  =  -25°C 

BIO  &  Bll  Additional  Large  Strain  Parameters  to  Supplement  Sheets  B5  to  B9 

B12  Parametric  Analysis  Illustrating  the  Relative  Importance  of  Testing 

Variables  on  the  Peak  Strength  of  Frozen  MFS 

B13  Normalized  Strength  Increase  vs.  Temperature  at  Varying  Strain  Rates 

for  Loose  and  Dense  Frozen  MFS 

B14  &  B15  Axial  Strain  at  Peak  Strength  vs.  Relative  Density  and  Confining  Pressure 
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TqLIi.  B  Lot^^  SifQin  BahQ'jiot  of  Sofu^crhd  Frozen  Sonds  (p3of4) 
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MEASURED  CONVENTIONAL  FROZEN  TESTS  RESULTS 
Test  Results  at  T  =  -10  "C  (from  Andersen  1991) 
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MEASURED  CONVENTIIINAL  FROZEN  TESTS  RESULTS 
Test  Results  at  T  =  -10  *C  (from  Andersen  1991) 


zn  > 


Bi  CD 

6-1  l-H  - 

W  rt 

ss  •  'U  CS 
s>  =:  X  o 

j  ua  rt  '^t-H 
C3  c=  s  >  X 


bOC-hOftbOOitiOCutiObOCi.  CuHh  WJ  O,  G,  bO  C^Mh  m-iHh  hOgLhOhObObO 
& 

Mh  bO  hO  bOHn  bOMn  bO  >  &Hh  Hh  d.  bOM-i  Sh  bO  bOMH  bO  hO  bO^n  bO 

uDoooococaoooo  oomt^;ou3C^J'^cacoTj<ooi-i'^^c^'»-<t*oa 

OCaOOOPOCOt^  Z^OOr-iOCOOOCOOOOCSOOOOCacO^ 


to  CDCOP^mZIZ^  i-H  CNCOnCD  CS  O^T-l 


E-1  Z  VJ  5< 

b«a  w  a 

-<  Z  >— I - 

:  W  C==  ^ 

,  (=j  aee 

erf  cyc- 


srf  ca  a  w 

w  i-a  >-H 

ffL.  W  C5 - — 

CU  HH  W  >»rf 

>-  erf  so- 


lOOOOffsacacsocaoQOOOoooocsoooo 


trtVC»c»^©QC5CSO'^'^^C<J'^'^C-0500O00t0CCOC0'^00^rH^ 

-^CO  CC‘^^t-ht:^COOOOO^;;^COM  £2 

^  CS  CS  ^  C^>  _ y-l  g^l  C^l  CS  CS  1~< 

^  coirt'rH'-^ooOino^'»-iocovr5'^OooT:<ooc30THC^ 

ir^i>^ao®^cj;oooc5«^c^t-^c500«cr:cio-r-^^cacs;^^0£^o 

w  t-H _ ^  ^ 

oo'*atG'ci'«r'®<ir'J«^cooo^ 

00ga<a500C>g50<S>0000'^0000 

OOCSOOCSO'r-tOOCi'COOCSJtra'^Oirtt^OOCS'rHC^C^dJOOT^O 

xnOOOOOOOi«l^'^^i^OOOOOOLftt^OOOuaoOCilGOOOOt*lOOOOO 


cn  Cs£3  &3  rf- 

erf  « 


e  o 

Z  O  tH 

hM  M  W  Z  X 

•<  5h  ca 


twC^.r-(»iT-iCieoc^i>-co‘>^‘>aC}oo  co'^i*':>inO'^^wococo 

*^c<?  c^CO«^<N^^®^CSCS«^®^  CO*^CCCOC^CO^COir3*<3<CO 

l>-Wao*»^t^'OO^OOO^t^COCO'»*^t^O«^C^  00  05  ^  M  CO  CO  o  . 


ub*<bir5^io®oO<5>vnT?*'*^^ir3co'C3'oiO' 

CSS0CN^CS«^’0'«^  <N  OJ  »  CO  CN  ’ 


r*ai  i-i  m  to  Id 
<c^co(rQcsg^(Mc^cocN 


ooc<icstoc5coc^csoomcsvnc^r-i^c5coc3cai-oioooo^ca 

c^ootorr^to  •  •  •  •oor-xoto  *001^00  *0000  •£2£2£2S£2£2 

'^cocococoo*^cat*-corococo<Ncocococ'acocooococococscico 

lOlO  COC^CSJC^-^  CO  CO 

OOC5tOt-tOCOOOOOt^iOT^'^(NOOt-«OOOOOOCOO»tO 

ooo^^-'-H  •  •  *  •c^T-(i— (  ‘C^CN  •i-HCScso-r-ics 

i-Hi-Hi-iC^JCO  CO  o  cs 


oot^t*-cooT:^c*-cs 

totoxrt^xoxrtxrtxrtto 


050505C505C505C505 


CQC^C0CSC0(N'-<C0'^^C^*7‘C0'^O'r-l'— 

1-(0”|-I000000000000000000000000 

o'^0^'^’^'^OOOOOtraOC^'^'^'^XOC^OXOOO**C^ 
•rHO^OOOO»-<THi-l  ^  «rt  ■<->  OOP _ 1-1  tH  ^  O  O 

TrO'^WCO^OxrtCOOOOOirsO'rHOOtOt--i-HC^T?CJ'rHb-tvO-pH 

t^tOCOOtOOOC^COtDtOOOC^'^OOCMtOCOO*^CO'«?''^^^COCO 

C005a505C0C0C5TrC5C50QOC5C5  00  C5^O_C5^C0_^n_C0j£_XX£Q£C£0O^2^ 

0^?aCO'^Vf5COt^0005^C>5CO“^UOtOt^C30C50^CSCO^^CDC^ 

i^J-ij^TrTrTr’^’^rr'^'^voioioiGxoxoxoxrDuototototototototo 


205 


Sheet  B6 


TABLE  5.3  (page  3  of  5) 

MEASURED  CONVENTIONAL  FROZEN  TESTS  RESULTS 
Tests  Results  at  T  =  -15  "C 
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MEASURED  CONVENTIONAL  FROZEN  TESTS  RESULTS 
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strain  at  Q 


APPENDIX  C:  Data  from  Triaxial  Compression  Tests  on  Unfrozen  MFS 


Sheet  No. 

C1&C2 

C3&C4 

C5 

C6 

C7 

C8&C9 

CIO 

Cll 

C12 


Description 

Shear  Characteristics  at  or  near  Steady  State  from  CIUC  and  CAUC  Tests 

Steady  State  Parameters  from  CIUC  Tests 

Steady  State  Parameters  from  CIDC  Tests 

Void  Ratio  vs.  Consolidation  Stress 

Effect  of  Consolidation  Stress  on  Upper  Yield  Stress 

Effect  of  Relative  Density  on  Effective  Stress  Paths  from 
CIUC  tests  with  =  0.1  and  2  MPa 

Effect  of  Relative  Density  on  Drained  Behavior  from 
CIDC  Tests  with  o’c  =  0.1  MPa 

Steady  State  Line  from  CIUC  and  CAUC  Tests 

Undrained  Strength  Ratio  vs.  State  Parameter  from  CIUC  Tests 
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SHEAR  CHARACTERISTICS  -  SERIES  C  ^  ^  ^ 

Steady  State  Parameters 
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Sheet  C5 


0.1  MPa 


Sheet  C7 


Figure  4.38  Upper  Yield  Stress  versus  Effective  Confining  Stress  for 

Undrained  Triaxial  Compression  Tests  (Swan  1994) 


4.44  Effective  Stress  Paths  for  Unfrozen  CIUC  Tests 

with  Similar  Confinement:  0.1  MPa  (Swan  1994) 


2  MPa 


Effective  Stress  Paths  for  Unfrozen  CIUC  Tests 
■with  Similar  Confinement:  2  MPa  (Swan  1994) 


Figure  4.49  Stress-Strain  and  Volumetric  Strain  Responses  for 

Unfrozen  CIDC  Tests  at  Similar  Confinement:  0.1  MPa 
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APPENDIX  D: 


Data  from  Consolidate-Freeze  Triaxial  Compression  Tests  on  MFS 


Sheet  No. 


D1  to  D3 


D4 

D5 

D6 

D7 

D8 

D9 


DIO 


Description 

Tabulated  Data  from  Consolidate-Freeze  Tests: 

Consolidation,  Freezing  and  Shearing  Results 

Compression  Curves  during  Consolidation 

Cell  Modifications  for  Consolidate-Freeze  Tests 

Pore  Water  Volume  and  Temperature  Changes  during  Freezing  of  CF02 
Efficiency  during  Freezing 

Comparison  of  E  and  Qy^,  from  Conventional  and  Consolidate-Freeze  Tests 

Stress-Axial  Strain-Volumetric  Strain  Curves  from 
Ten  Consolidate-Freeze  Tests 

Qp  vs.  Confining  Stress  from  Conventional  and  Consolidate-Freeze  Tests 
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lAHLCi  5.7  /  ^  lQoj\ 

CONSOLIDATION  CHARACTERISTICS  ln^4) 


TABLE  5.8 

SPECIMEN  FREEZING  CHARACTERISTICS 
CONSOLIDATE- FREEZE  TESTS 
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TABLE  5.9 

CONSOLIDATE- FREEZE  TESTS  RESULTS 
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Sheet  D3 


Isotropically  Consolidated  Consolidate— Freeze  Tests  (Swan  1994) 


Figure  3.4 


Specimen 


Details  of  Internal  Modifications  for  Consolidate-Freeze  Testing  Program 

(Sw^n  Ii94) 
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Sheet  D5 


Efficiency  of  Specimen  Freezing  versus  Elapsed  Freezing  Time 
for  Consolidate— Freeze  Tests  on  Dense  Specimens  (Swan  1994^ 


Volumetric  Strain,  e  (%)  Deviator  Stress,  a.  -  (MPa) 


Deviator  Stress  vs  Axial  Strain 


Axial  Strain  (%) 


Volumetric  Strain  vs  Axial  Strain 


Axial  Strain  (%) 

Figure  5.58  Comparison  of  Stress-Strain  Curves  for  Consolidate- 
Freeze  Tests  on  Dense  Specimens  at  Moderate  Strain 
Rate  and  T  =  -lO^C  (Swan  1994) 
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APPENDIX  E: 


Application  of  Dilatancy-Hardening  Model  to  Frozen  MFS 


Sheet  No. 

EltoE3 

E4 

E5 

E6 

E7 

E8 

E9 

ElO 

Ell 

E12 


Description 

Data  on  the  Tensile  Strength  of  Polycrystalline  Ice 

Compressive  Strength  of  Polycrystalline  Ice  at 
Oc  =  0  and  10  MPa  from  Jones  (1982) 

Unconfined  Compressive  Strength  of  Isotropic  Polycrystalline  Ice  from 
Constant  Stress  and  Constant  Strain  Rate  Test  Programs 

Compressive  Strength  of  Polycrystalline  Ice  Used  in  Dilatancy  Hardening 
Model  Compared  to  Results  of  Murrell  et  at  (1989) 

Effect  of  Confinement  on  Stress-Strain  Behavior  of  Ice 

Ratio  of  Dilatancy-Hardening  Model  Predicted  to  Measured 
Peak  Strength  for  Frozen  MFS 

Ratio  of  Predicted  Ice  Strength  to  Measured  Upper  Yield  Stress  of 
Frozen  MFS 

Modified  D-H  Model  Predicted  Strength  minus  Measured 
Peak  Strength  for  Frozen  MFS 

Comparison  of  D-H  Model  Predicted  and  Measured  Peak  Strength  vs. 
Relative  Density  at  Varying  Confining  Pressures 

Relative  Density  vs.  Critical  Confining  Pressure  and 
Undrained  Strength  for  MFS 
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Figure  2.46  Stress— strain  Curves  From  Uniaxial  Tension  Tests  on 

Polycrystalline  Ice  (from  Hawkes  and  MeUor  1972) 
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Sheet  El 


Granular  Ice  at  —7*  C  (from  Mellor  1980) 


urre 


CCL  Ij Zol9S__A^0^ 


T— 26*C  e-lO'*s**- 


Fig.  6.  Stress-strain  curves 
for  pure  ice  specimens  under 
triaxial  compression  at  -20 
deg  c  and  a  strain-rate  of 
lO'Vs . 


T— 20’C  e-io'^s'* 


Fig.  7.  Stress-strain  curves 
for  pure  ice  in  triaxial 
compression  at  -20  deg  c  and 
a  strain-rate  of  lO'^/s.  p  » 
confining  pressure. 
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Comparison  of  Peak  Strength  versus  Relative  Density  Estimated  from 
the  Dilatancy-Hardening  Model  and  Measured  in  Conventional  Frozen  MFS 
Tests  at  Moderate  Strain  Rate.  T=-10°  and  Various  Confining  Pressures 
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